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1 Basic relations 

A cosmological model represents the Universe at a particular scale. We will assume that on large scales, 
space-time geometry is described by Einstein's general theory of relativity (see, e.g., d'Inverno [IJ, Wald 
lIU, Hawking and Ellis 1,3], or Stephani [4J). Then a cosmological model is defined by specifying 15J-[|7J: 

* the space-time geometry represented on some specific averaging scale and determined by the metric 
9ab{x^), which — because of the requirement of compatibility with observations — must either have some 
expanding Robertson-Walker ('RW') geometries as a regular limit (see |l8l), or else be demonstrated to 
have observational properties compatible with the major features of current astronomical observations of the 
Universe; 

* the matter present, represented on the same averaging scale, and its physical behaviour (the energy- 
momentum tensor of each matter component, the equations governing the behaviour of each such component, 
and the interaction terms between them), which must represent physically plausible matter (ranging from 
early enough times to the present day, this will include most of the interactions described by present-day 
physics); and 

* the interaction of the geometry and the matter — how matter determines the geometry, which in 
turn determines the motion of the matter (see e.g. [9]). We assume this is through Einstein's relativistic 
gravitational field equations ('EFE') given b}Q 

Gab = Rab ~ \ Rgab = Tab — ^ dab , (1) 

which, because of the twice-contracted Bianchi identities, guarantee the conservation of total energy- 
momentum 

VbG"'^==0 ^ VbT''''^0, (2) 
provided the cosmological constant A satisfies the relation VqA — 0, i.e., it is constant in time and space. 

Together, these determine the combined dynamical evolution of the model and the matter in it. The 
description must be sufficiently complete to determine 

* the observational relations predicted by the model for both discrete sources and background radiation, 
implying a well-developed theory of structure growth for very small and for very large physical scales (i.e., 
for light atomic nuclei and for galaxies and clusters of galaxies), and of radiation absorbtion and emission. 

To be useful in an explanatory role, a cosmological model must be easy to describe — that means they 
have symmetries or special properties of some kind or other The usual choices for the matter description will 
be some combination of 

* a fluid with a physically well-motivated equation of state, for example a perfect fluid with specified 
equation of state (beware of imperfect fluids, unless they have well-defined and motivated physical proper- 
ties); 

* a mixture of fluids, usually with different 4-velocities; 

* a set of particles represented by a kinetic theory description; 

* a scalar field (p, with a given potential ¥{(/)) (at early times); 

* an electromagnetic field described by Maxwell's field equations. 

As intimated above, the observational relations implied by cosmological models must be compared with 
astronomical observations. This determines those solutions that can usefully be considered as viable cos- 
mological models of the real Universe. A major aim of the present lectures is to point out that this class 

'Throughout this review we employ geometrised units characterised by c = 1 = 8itG/c^. Consequently, all geometrical variables 
occuning have physical dimensions that are integer powers of the dimension [ length ] . The index convention is such that space-time and 
spatial indices with respect to a general basis are denoted by a, 6, ■ • • = 0, 1, 2, 3 and o, /3, ■ • ■ = 1, 2, 3, respectively, while space-time 
and spatial indices in acoordinate basis are fj,,u, ■■■ = 0,1,2,3 and ■•• = 1,2, 3,, respectively. 
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is wider than just the standard Friedmann-Lemaitre-Robertson- Walker ('FLRW') cosmologies; indeed 
those models cannot be reahstic on all scales of description, but they may also be inaccurate on large scales, 
or at very early and very late times. To examine this, we need to consider the subspace of the space of all 
cosmological solutions that contains models with observational properties like those of the real Universe at 
some stage of their histories. Thus we are interested in the full state space of solutions, allowing us to see how 
realistic models are related to each other and to higher symmetry models, including particularly the FLRW 
models. 

These lectures develop general techniques for examining this, and describe some specific models of in- 
terest. The first part looks at exact general relations valid in all cosmological models, the second part at exact 
cosmological solutions of the EFE, and the third part at approximate equations and solutions: specifically, 
'almost-FLRW models, linearised about a FLRW geometry. 

2 1 + 3 covariant description 

Space-times can be described via 

(a) the metric gij{x^) described in a particular set of local coordinates, with its differential properties, 
as embodied by the connection, given through the Christoffel symbols; 

(b) the metric described by means of particular tetrads, with its connection given through the Ricci 
rotation coefficients; 

(c) 1 + 3 covariantly defined variables. In anisotropic cases, tetrad vectors can be uniquely defined in a 
1 + 3 covariant way and this approach merges into (b). 

Here we will concentrate on the 1 + 3 covariant approach, based on ll5llToll6ll7i rTT1l . but dealing also 
with the tetrad approach which serves as a completion to the 1 + 3 covariant approach. The basic point here 
is that because we have complete coordinate freedom in General Relativity, it is preferable where possible to 
describe physics and geometry by tensor relations and quantities; these then remain valid whatever coordinate 
system is chosen. 

2.1 Variables 

2.1.1 Average 4-velocity of matter 

In a cosmological space-time (Al, g), at late times there will be a family of preferred worldlines representing 
the average motion of matter at each poin|l (notionally, these represent the histories of clusters of galaxies, 
with associated 'fundamental observers'); at early times there will be uniquely defined notions of the average 
velocity of matter (at that time, interacting gas and radiation), and corresponding preferred worldlines. In 
each case their 4-velocity i^ 

"° = , UaU'^ = ~ 1 , (3) 

dr 

where r is proper time measured along the fundamental worldlines. We assume this 4-velocity is unique: that 
is, there is a well-defined preferred motion of matter at each space-time event. At recent times this is taken 
to be the 4-velocity defined by the vanishing of the dipole of the cosmic microwave background radiation 
('CBR'): for there is precisely one 4-velocity which will set this dipole to zero. It is usually assumed that this 
is the same as the average 4-velocity of matter in a suitably sized volume ||6l; indeed this assumption is what 
underlies studies of large scale motions and the 'Great Attractor' . 

^We are here assuming a fluid description can be used on a large enough scale [5 6|. The alternative is that the matter distribution is 
hierarchically structured at all levels or fractal (see e.g. 1 12 1 and refrences there), so that a fluid description does not apply. The success 
of the FLRW models encourages us to use the approach taken here. 

'Merging from the one concept to the other as structure formation takes place. 
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Given u"-, there are defined unique projection tensors 

U\ ^-u'^Ub ^ U\U'b = U\ , C/"a = 1 , Uabu'' = Ua , (4) 

hab = 5afc + Ua Mfc ^ /l"c /I'^b = /l^b , /l^a = 3 , hab u'' ^ . (5) 

The first projects parallel to the 4-velocity vector u", and the second determines the (orthogonal) metric 
properties of the instantaneous rest-spaces of observers moving with 4-velocity it". There is also defined a 
volume element for the rest-spaces: 

Vabc = U'^ Vdabc ^ Vabc = »7[afcc] , rjabc = , (6) 

where r]abcd is the 4-dimensional volume element (r]abcd. = Vlabcd]^ '7oi23 = y/\ det gab |)- 

Moreover, two derivatives are defined: the covariant time derivative ' along tlie fundamental world- 
lines, where for any tensor T^^cd 

f'^'cd = u'^V.T-'.d , (7) 
and the fully orthogonally projected covariant derivative V, where for any tensor T°-''cd 

S/eT^'cd = h'^f h\ hP, h^d T^'pg , (8) 

with total projection on all free indices. The tilde serves as a reminder that if has non-zero vorticity, V 
is not a proper 3-dimensional covariant derivative (see Eq. dZTl i below). Finally, following fTT] (and see 
also 1 13 |), we use angle brackets to denote orthogonal projections of vectors and the orthogonally projected 
symmetric trace-free part of tensors: 

/^^^;^ T<"^> = [ /^(^ /i'')^ - i /i''^ /i^d ] T^^ (9) 

for convenience the angle brackets are also used to denote othogonal projections of covariant time derivatives 
along M° (Fermi derivatives): 

v'^'''' =h''bv\ f^''''^ ^[h^\h''^d-^h''''h^d]T'"^ ■ (10) 
Exercise: Show that the projected time and space derivatives of Uab, hab and rjabc all vanish. 

2.1.2 Kinematical quantities 

We split the first covariant derivative of Ua into its irreducible parts, defined by their symmetry properties: 

Valifc ^ -UaUb + ^aUb = - lia Ufc + | 6 hab + (^ab + ^ab , (H) 

where the trace 9 = Vau" is the rate of volume expansion scalar of the fluid (with H = 6/3 the Hubble 

scalar); aab — ^{a^b) is the trace-free symmetric rate of shear tensor (aab — <y(ab)^ <^ab — 0, <7°-a — 0), 
describing the rate of distortion of the matter flow; and LOab — '^[a^b] is the skew-symmetric vorticity 
tensor {ujab = ^^[ab]^ ^abU^ — 0)Q describing the rotation of the matter relative to a non-rotating (Fermi- 
propagated) frame. The stated meaning for these quantities follows from the evolution equation for a relative 
position vector 77^ ~ h°-b'rf', where 77° is a deviation vector for the family of fundamental worldlines, i.e., 
u^Vb'n'^ — rf'Vbu"^ ■ Writing 77" = 51 e"", Cae" = 1, we find the relative distance 5t obeys the propagation 
equation 

^ = ie + K,e'^e^), (12) 
(the generalised Hubble law), and the relative direction vector the propagation equation 



' = (a b- [cTcde e ) h b - cj f,) e , 

'The vorticity here is defined with respect to a right-handedly oriented spatial basis. 



(13) 
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giving the observed rate of change of position in the sky of distant galaxies. Finally = u^Vbu"" is the 
relativistic acceleration vector, representing the degree to which the matter moves under forces other than 
gravity plus inertia (which cannot be covariantly separated from each other in General Relativity: they are 
different aspects of the same effect). The acceleration vanishes for matter in free fall (i.e., moving under 
gravity plus inertia alone). 



2.1.3 Matter tensor 

The matter energy-momentum tensor Tat can be decomposed relative to u"- in the forrrH 

Tab = HUaUb + qaUb+ Ua Qb + P Kb + T^ab , (14) 
QaU" = , TT^Q = , TTab = TT(ab) , TTab = , 

where ^ = (Tabu'^u'') is the relativistic energy density relative to u^, q°- = — Tbc u'' is the relativistic 
momentum density, which is also the energy flux relative to u"^, p — ^ {Tabh"''') is the isotropic pressure, 
and TTab = Ted h'^ {a h'^b) is the trace-free anisotropic pressure (stress). 

The physics of the situation is in the equations of state relating these quantities; for example, the com- 
monly imposed restrictions 

= TTafc = Tab = "a "6 + P ^ab (15) 

characterise a perfect fluid with, in general, equation of state p — p{fi, s). If in addition we assume that 
p = 0, we have the simplest case: pressure-free matter ('dust' or 'Cold Dark Matter'). Otherwise, we must 
specify an equation of state determining p from /i and possibly other thermodynamical variables. Whatever 
these relations may be, we usually require that various energy conditions hold: one or all of 

H>0, (Ai + p) > , (/^ + 3p) > , (16) 

(the latter, however, being violated by scalar fields in inflationary universe models), and additionally demand 
the isentropic speed of sound — {dp/dfi)s=const obeys 

< < 1 ^ < ( 1^ ) < 1 , (17) 

\ / s=const 

as required for local stability of matter (lower bound) and causality (upper bound), respectively. 



2.1.4 Maxwell field strength tensor 

The MaxweU field strength tensor Fab of an electromagnetic field is split relative to w° into electric and 
magnetic field parts by the relations (see Q) 

Ea=Fabu'' ^ -BaW"=0, (18) 
Ha = ^VabcF'" ^ HaU-' = Q. (19) 



2.1.5 Weyl curvature tensor 

In analogy to Fab, the Weyl conformal curvature tensor Cabcd is split relative to into 'electric' and 
'magnetic' Weyl curvature parts according to 

Eab = Cacbd U^U'' => E\^0, Eab= E^ab) , EabU^'^O, (20) 
Hab = \ Vade C'^'^bc u'^ +> H'^ a = , Hab = F[(^ab) i Hab u'^ — ■ (21) 
^We should really write fi = fi{u"'), etc; but usually assume this dependence is understood. 



2 1 + 3 COVARIANT DESCRIPTION 



9 



These represent the 'free gravitational field', enabling gravitational action at a distance (tidal forces, gravi- 
tational waves), and influence the motion of matter and radiation through the geodesic deviation equation 
for timelike and null congruences, respectively lfT4l - |fT8l . Together with the Ricci curvature tensor Rat 
(determined locally at each point by the matter tensor through Einstein's field equations ([T])), these quantities 
completely represent the space-time Riemann curvature tensor Rabcd, which in fully 1 + 3-decomposed 
form become^ 

T->ab rtab , j^ab , nab , nab 

n cd — cd + cd + cd + i^H cd , 

Rpcd = |(M + 3p-2A)7.["u[,;j''l<i] + |(/i + A)/i'^[,/i^], 

Rfcd = -2u['^/i^l[,grf]-2u[,/i[><z''l-2w['^U[,7r'']d]+2/i['^[,7r^ld], (22) 
R'icd = ^u^''u^cE''^d]+^h^\cE''^d], 

R'k'cd = 2 7?'^''^M[,i/d],+2 7?edeM["-ff'l" . 



2.1.6 Auxiliary quantities 

It is useful to define some associated kinematical quantities: the vorticity vector 

c^" = i77"^^Wte ^ w„u" = 0, w„fcc^'' = 0, (23) 

the magnitudes 

= i Kfco;'^'') > , = i Kba"') > ' (24) 
and the average length scale S determined by 

|-ie, (25) 

so the volume of a fluid element varies as S^. Further it is helpful to define particular spatial gradients 
orthogonal to characterising the inhomogeneity of space-time: 

Xa = Vaf^ , Za = V aQ ■ (26) 

The energy density /i (and also 8) satisfies the important commutation relation for the V-derivative ||T9l 

V[aVf,]/X = T^ahc^^ A ■ (27) 

This shows that if uj"- (j. ^ in an open set, then Xa ^ there, so non-zero vorticity implies anisotropic 
number counts in an expanding universe f20l (this is because there are then no 3-surfaces orthogonal to the 
fluid flow; see [6j|). 



2.2 1 + 3 covariant propagation and constraint equations 

There are three sets of equations to be considered, resulting from Einstein's field equations ^ and their 
associated integrability conditions. 

*Here P is the perfect fluid part, / tlie imperfect fluid part, E that due to the electric Weyl curvature, and H that due to the magnetic 
Weyl curvature. This obscures the similarities in these equations between E and tt, and between H and q; however, this partial symmetry 
is broken by the field equations, so the splitting given here (due to M Shedden) is conceptually useful. 
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2.2.1 Ricci identities 

The first set arise from the Ricci identities for the vector field u", i.e., 

2V[,Vb]?/"-i?a6"rfU^ (28) 

On substituting in from (fTTT i. using ([T]), and separating out the orthogonally projected part into trace, symmet- 
ric trace-free, and skew symmetric parts, and the parallel part similarly, we obtain three propagation equations 
and three constraint equations. The propagation equations are, 

1. The Raychaudhuri equation |21 1 

e - Wau"" = - i + {iiau") - 2a^ + 2uj^ - ^i^i + 3p)+A, (29) 

which is the basic equation of gravitational attraction im-Q, showing the repulsive nature of a positive 
cosmological constant, leading to identification of (/i + 3p) as the active gravitational mass density, and 
underlying the basic singularity theorem (see below). 

2. The vorticity propagation equation 

together with (|38T l below, showing how vorticity conservation follows if there is a perfect fluid with acceler- 
ation potential $ |5, 7 1, since then, on using (|27]l, ■q"-'"' VbUc = "q"-^" VbVc$ = 2 w° <i>, 

3. The shear propagation equation 

a{o.b) _ ^{a^b) ^_2Q^ab^ ^{a -^b) _ ^{a^ ^b)c _ ^{a ^b) _ (^ab _ 1 ^ab^ 

the anisotropic pressure source term iVab vanishing for a perfect fluid; this shows how the tidal gravitational 
field, the electric Weyl curvature Eab, directly induces shear (which then feeds into the Raychaudhuri and 
vorticity propagation equations, thereby changing the nature of the fluid flow). 

The constraint equations are, 

1. The (Oa) -equation 

= {CiT = Vhd"^ - I V'^e + 77"^^^ [ VbUJ, + 2 ?ib CO, ] + , (32) 

showing how the momentum flux (zero for a perfect fluid) relates to the spatial inhomogeneity of the expan- 
sion; 

2. The vorticity divergence identity 

= (C2) = VaUj'' - {uaLo'') ; (33) 

3. The iJab-equation 

= (C3)"'' = H°-^ + 2u<"tj''> + V^'^cj''^ - (curlCT)"'' , (34) 

characterising the magnetic Weyl curvature as being constructed from the 'distortion' of the vorticity and the 
'curl' of the shear, (curlo-)"^ = if'^^''Vca^'>d- 

2.2.2 Twice-contracted Bianclii identities 

The second set of equations arise from the twice-contracted Bianchi identities which, by Einstein's field 
equations ([TJ, imply the conservation equations Projecting parallel and orthogonal to u", we obtain the 
propagation equations 

A + Vaq" =-eifl+p)-2 [iLaq") - (aabTT''') (35) 
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and 

q{a) + v> + Vb^''^ = -leq''^a%q'' -i^i+p)u''^ iit n'^' - 77"''^ q, , (36) 

which constitute the energy conservation equation and the momentum conservation equation, respec- 
tively. For perfect fluids, characterised by Eq. ( fTSl l. these reduce to 

A = -e(M + p), (37) 

and the constraint equation 

O^VaP+i^i+p)Ua. (38) 

This shows that (/i + p) is the inertial mass density, and also governs the conservation of energy. It is clear 
that if this quantity is zero (an effective cosmological constant) or negative, the behaviour of matter will be 
anomalous. 

Exercise: Examine what happens in the two cases (i) {fJ. + p) = 0, (ii) (/.t + p) < 0. 
2.2.3 Other Bianchi identities 

The third set of equations arise from the Bianchi identities 

V[ai?bc]de = . (39) 

Double contraction gives Eq. (|2|l, akeady considered. On using the splitting of Rated into Rab and Cabcd, 
the above 1 + 3 splitting of those quantities, and Einstein's field equations, the once-contracted Bianchi 
identities give two further propagation equations and two further constraint equations, which are similar in 
form to Maxwell's field equations in an expanding universe (see 1221 171). 
The propagation equations are, 

+ i7r<"''>) - (curli/)"'' + iV^^'g^^ = - ^ {p + p) a"'' - Q {E"'' + ^ t:"'') (40) 

+ Ij'''^- [ 2 + LO, + i ) ] , 

the i?-equation, and 

ij^'^''^ + (curl i;)''^- i (curl tt)'''' = - 9 + 3 ct^'^c i?''^" + | t^^N''^ (41) 



the ij-equation, where we have defined the 'curls' 



(curli?)'''' = Vci^/^ (42) 

(curl^;)''^ = Vc£;d''^ (43) 

(curlTr)"^ = 'n'^^^^'S/cnd''^ . (44) 

These equations show how gravitational radiation arises: taking the time derivative of the £^-equation gives 
a term of the form (curl H)' ; commuting the derivatives and substituting from the i/-equation eliminates H, 
and results in a term in E and a term of the form (curl curl E), which together give the wave operator acting 
on E II22I I23I: similarly the time derivative of the iJ-equation gives a wave equation for H. 
The constraint equations are 

= = Wb{E'''' + ^TT'''')-lW''fi+lQq''-^(J%q''-3uJbH''^ 

-7j''''^[abdH/-lcjbqc], (45) 
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the (div i?)-equation with source the spatial gradient of the energy density, which can be regarded as a vector 
analogue of the Newtonian Poisson equation |24|, enabling tidal action at a distance, and 

= {C5T = ^bH''' + {fi+p)L0'' + 'SLUb{E''''-lTT'''') 

+ rj-'^ [ i Vb9c + <yM {E/ + i tt/) ] , (46) 

the (div i?) -equation, with source the fluid vorticity. These equations show respectively that scalar modes 
will result in a non-zero divergence of Eat (and hence a non-zero _E-field), and vector modes in a non-zero 
divergence of Hab (and hence a non-zero 7J-field). 

2.2.4 Maxwell's field equations 

Finally, we turn for completeness to the 1 + 3 decomposition of Maxwell's field equations 

Vf,^^"^=j:, V[aFbc]=0. (47) 

As shown in Q, the propagation equations can be written as 

_ ^-fc- Vbi/c = -j^'''^ -lQE'' + a\E'' + Tj''''[ubH,+uJbE,], (48) 

^(a> ^ ^ahc ^ -^OH'' +a\H' -7f'''[ubE,-L0bH,], (49) 

while the constraint equations assume the form 

- (Ce) - VaE"" ~ 2 [LOaH'') ~ Pc , (50) 
= {Ch) = VaH'' +2{uJaE'') , (51) 

where po = (-joau")- 

2.3 Pressure-free matter ('dust') 

A particularly useful dynamical restriction is 

= p = q^^nab ^ iia^O, (52) 

so the matter (often described as 'baryonic') is represented only by its 4-velocity and its energy density 
/i > 0. The implication follows from momentum conservation: ([38]l shows that the matter moves geodesically 
(as expected from the equivalence principle). This is the case of pure gravitation: it separates out the (non- 
linear) gravitational effects from all the fluid dynamical effects. The vanishing of the acceleration greatly 
simplifies the above set of equations. 

2.4 Irrotational flow 

If we have a barotropic perfect fluid: 

= q'' = nab, P = Pifj-) , v"''" VbUe = , (53) 

then cj" = is involutive: i.e., 

tj" = initially = w'' = at all later times 

foHows from the vorticity conservation equation dSOl l (and is true also in the special case p = 0). When the 
vorticity vanishes: 

1. The fluid flow is hypersurface-orthogonal, and there exists a cosmic time function t such that Ua = 
— g{x'') Vat; if additionally the acceleration vanishes, we can set g to unity; 



2 1 + 3 COVARIANT DESCRIPTION 



13 



2. The metric of the orthogonal 3-spaces is hat, 

3. From the Gauss embedding equation and the Ricci identities for the Ricci tensor of these 3-spaces 
is given by |5 , 6| 

^Rab = -C7(a6) -6 0-ab + V(„Mb) +M(aWb) +7rafc+ i/lah [2 /X- §92 + 2 0-2 + 2 A] , (54) 

which relates '^Rab to Eab via ( l3Tl i. and their Ricci scalar is given by 

^i? = 2^-| 6^ + 2 + 2A, (55) 

which is a generalised Friedmann equation, showing how the matter tensor determines the 3-space average 
curvature. These equations fully determine the curvature tensor ^Rabcd of the orthogonal 3-spaces Q. 

2.5 Implications 

Altogether, in general we have six propagation equations and six constraint equations; considered as a set of 
evolution equations for the 1 + 3 covariant variables, they are a first-order system of equations. This set is 
determinate once the fluid equations of state are given; together they then form a complete set of equations 
(the system closes up, but is essentially infinite dimensional because of the spatial derivatives that occur). 
The total set is normal hyperbolic at least in the case of a perfect fluid, although this is not obvious from the 
above form; it is shown by completing the equations to tetrad form (see the next section) and then taking 
combinations of the equations to give a symmetric hyperbolic normal form (see |25, 26|). We can determine 
many of the properties of specific solutions directly from these equations, once the nature of these solutions 
has been prescribed in 1 + 3 covariant form (see for example the FLRW and Bianchi Type 1 cases considered 
below). 

The key issue that arises is consistency of the constraints with the evolution equations. It is believed 
that they are generally consistent for physically reasonable and well-defined equations of state, i.e., they are 
consistent if no restrictions are placed on their evolution other than implied by the constraint equations and 
the equations of state (this has not been proved in general, but is very plausible; however, it has been shown 
for irrotational dust [11, 27J). It is this that makes consistent the overall hyperbolic nature of the equations 
with the 'instantaneous' action at a distance implicit in the Gauss-like equations (specifically, the (div E)- 
equation), the point being that the 'action at a distance' nature of the solutions to these equations is built into 
the initial data, which ensures that the constraints are satisfied initially, and are conserved thereafter because 
the time evolution preserves these constraints (cf. |28 |). A particular aspect of this is that when cj" — 0, the 
generalised Friedmann equation ( |55] ) is an integral of the Raychaudhuri equation ( |29] l and energy equation 

m- 

One must be very cautious with imposing simplifying assumptions (such as, e.g., vanishing shear) in order 
to obtain solutions: this can lead to major restrictions on the possible flows, and one can be badly misled 
if their consistency is not investigated carefully 1291 l24ll . Cases of particular interest are shear- free fluid 
motion (see ll30ll - ll32ll ) and various restrictions on the Weyl curvature tensor, including the 'silent universes', 
characterised by Hab = (and p = 0) t33,i34J, or models with Vhi/"'' = ll35l . 

2.5.1 Energy equation 

It is worth commenting here that, because of the equivalence principle, there is no agreed energy conservation 
equation for the gravitational field itself, nor is there a definition of its entropy (indeed some people — Free- 
man Dyson, for example ||36l — claim it has no entropy). Thus the above set of equations does not contain 
expressions for gravitational energjQ or entropy, and the concept of energy conservation does not play the 



^ There ai'e some proposed 'super-energy' tensors, e.g., the Bel-Robinson tensor 1371 . but they do not play a significant role in the 
theory. 
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major role for gravitation that it does in the rest of physics, neither is there any agreed view on the growth 
of entropy of the gravitational field0 However, energy conservation of the matter content of space-time, ex- 
pressed by the divergence equation V^T"^^ — 0, is of course of major importance. 

If we assume a perfect fluid with a (linear) 7-law equation of state, then ( l37l ) shows that 

p = (7 -1)^,7 = ^ M = M/S-^T , M = . (56) 

One can approximate ordinary fluids in this way with 1 < 7 < 2 in order that the causality and energy 
conditions are valid, with 'dust' and Cold Dark Matter ('CDM') corresponding to 7 = 1 /i = M/ S^, and 
radiation to 7 = | =4' /i = M/ S^. 

Exercise: Show how to generalise this to more realistic equations of state, taking account of entropy and 
of matter pressure (see e.g. H-Q). 

In the case of a mixture of non-interacting matter, radiation and CDM having the same 4-velocity, rep- 
resented as a single perfect fluid, the total energy density is simply the sum of these components: ^ = 
Mdust + McDM + A*radn- (NB: This is only possible in universes with spatially homogeneous radiation energy 
density, because the matter will move on geodesies which by the momentum conservation equation implies 
VaPiadn = <^ Va/^radn = 0. This will not be truc for a general inhomogeneous or perturbed FLRW 
model, but will be true in exact FLRW and orthogonal Bianchi models.) 

Exercise: The pressure can still be related to the energy density by a 7-law as in (|56] l in this case of non- 
interacting matter and radiation, but 7 will no longer be constant. What is the equation giving the variation 
of (i) 7, (ii) the speed of sound, with respect to the scale factor in this case? (See (381.) 

A scalar field has a perfect fluid energy-momentum tensor if the surfaces {4> ^ const} are spacelike and 
we choose normal to these surfaces. Then it approximates the equation of state (|56] | in the 'slow -rolling' 
regime, with 7 w 0, and in the velocity-dominated regime, with 7 « 2. In the former case the energy condi- 
tions are no longer valid, so 'inflationary' behaviour is possible, which changes the nature of the attractors in 
the space of space-times in an important way. 

Exercise: Derive expressions for /i, p, [ji + p), (/i + ip) in this case. Under what conditions can a scalar 
field have (a) {p. + p) = 0, (b) (/x + 3p) = 0, (c) (/i + 3p) < 0? 

2.5.2 Basic singularity theorem 

Using the definition (IZST i of S, the Raychaudhuri equation can be rewritten in the form (cf. El]) 

3 I = - 2 (fT^ _ + v^-^ + {uau") - i (a* + 3p) + A , (57) 

showing how the curvature of the curve S{t) along each worldline (in terms of proper time r along that 
worldline) is determined by the kinematical quantities, the total energy density and pressur^ in the combi- 
nation [ji + 3p), and the cosmological constant A. This gives the basic 

Singularity Theorem 11211 151161171: In a universe where + 3p) > 0, A < 0, and ii"- = 

tj" = at all times, at any instant when Hq = i > 0, there must have been a time < 1/ ^^o 
ago such that S* ^ as < ^ to; a space-time singularity occurs there, where fi ^ 00 and p ^ 00 
for ordinary matter (with (/i + p) > 0). 



Entropy is well understood in the case of black holes, but not for gravitational fields in the expanding Universe. 
'This form of the equation is valid for impeii'ect fluids also: the quantities and nab do not directly enter this equation. 
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The further singularity theorems of Hawking and Penrose (see ||39l |3] |40l) utilize this result (and its null 
version) as an essential part of their proofs. 

Closely related to this are two other results: the statements that (a) a static universe model containing 
ordinary matter requires A > (Einstein's discovery of 1917), and (b) the Einstein static universe is unstable 
(Eddington's discovery of 1930). Proofs are left to the reader; they follow directly from ( |57] ). 



2.5.3 Relations between important parameters 

Given the definitions 



1 5*0 Mo po A 



for the present-day values of the Hubble scalar ('constant'), deceleration parameter, density parameter, pres- 
sure to density ratio, and cosmological constant parameter, respectively, then from ( |57] i we obtain 

90 = ^ 772 TTri 1- 2 i'o (,1 + owoj - "A ■ (59) 

Now CBR anisotropics let us deduce that the first two terms on the right are very small today, certainly less 
than lO^'^, and we can reasonably estimate from the nature of the matter that po ^ piQ and the third term on 
the right is also very small, so we estimate that in realistic Universe models, at the present time 

go « i f^o - f^A • (60) 

(Note we can estimate the magnitudes of the terms which have been neglected in this approximation.) This 
shows that a cosmological constant can cause an acceleration (negative go); if it vanishes, as commonly 
assumed, the expression simplifies: 

A = ^ qo~k^o, (61) 

expressing how the matter density present causes a deceleration of the Universe. If we assume no vorticity 
(cj" = 0), then from (1551 ) we can estimate 

3i?o~6i7o'(r!o-l + f^A), (62) 

where we have dropped a term (ao/Ho)'^. If A = 0, then •^i?o ~ 6 Hq {Hq ~ 1), showing that flo = 1 is 
the critical value separating irrotational universes with positive spatial curvature (ilo > 1 ^Rq > 0) from 
those with negative spatial curvature (ilo < 1 ^Ro < 0). 

Present day values of these parameters are almost certainly in the ranges f4T]: baryon density: 0.01 < 
Qbaryons ^ 0.03, total matter density: 0.1 < flo < 0.3 to 1 (implying that much matter may not be baryonic), 
Hubble constant: 45 km/sec/Mpc < Hf) < 80 km/sec/Mpc, deceleration parameter: —0.5 < Qq < 0.5, 
cosmological constant: < fl^ < 1- 



2.6 Newtonian case 

Newtonian equations can be developed completely in parallel PSl l43l |6l and are very similar, but simpler; 
e.g., the Newtonian version of the Raychaudhuri equation is 

e + ie2 + 2(cr2_t^2^_Dc<a" + ip-A = 0, (63) 

where p is the matter density and Ua — Va + Dq$ is the Newtonian analogue of the relativistic 'acceleration 
vector', with ' the convective derivative and $ the Newtonian potential (with suitably generalised boundary 
conditions 1,44. ,45 J ). The Newtonian analogue of Eab is 

E^p = D„D^$ - i (D^DT$) /i„/3 , (64) 
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where ha/s denotes the metric, and Dq, the covariant derivative, of EucHdean space. For the latter Dahp~^ = 
and [Da, ] =0. There is no analogue of Hab in Newtonian theory |6|, as shown by a strict limit process 
leading from relativistic to Newtonian solutions ||46l . 

Exercise: Under what conditions will a relativistic cosmological solution allow a representation (|64| | for 
the electric part of the Weyl curvature tensor? Will the potential $ occurring here necessarily also relate to 
the acceleration of the timelike reference worldlines? 

2.7 Solutions 

Useful solutions are defined by considering appropriate restrictions on the kinematical quantities, Weyl cur- 
vature tensor, or space-time geometry, for a specified plausible matter content. Given such restrictions, 

(a) we need to understand the dynamical evolution that results, particularly fixed points, attractors, etc., 
in terms of suitable variables, 

(b) we particularly seek to determine and characterise involutive subsets of the space of space-times: 
these are subspaces mapped into themselves by the dynamical evolution of the system, and so are left invariant 
by that evolution. The constraint and evolution equations must be consistent with each other on such subsets. 
A characterisation of these subspaces goes a long way to characterising the nature of self-consistent solutions 
of the full non-linear EFE. 

As far as possible we aim to do this for the exact equations. We are also concerned with 

(c) linearisation of the equations about known simple solutions, and determination of properties of the 
resulting linearised solutions, in particular considering whether they accurately represent the behaviour of the 
full non-linear theory in a neighbourhood of the background solution (the issue of linearisation stability), 

(d) derivation of the Newtonian limit and its properties from the General Relativity equations, and under- 
standing how accurately this represents the properties of the full relativistic equations (and of its linearised 
solutions); see 1,24 J for a discussion. 

3 Tetrad description 

The 1 + 3 covariant equations are immediately transparent in terms of representing relations between 1 + 3 
covariantly defined quantities with clear geometrical and/or physical significance. However, they do not form 
a complete set of equations guaranteeing the existence of a corresponding metric and connection. For that we 
need to use a tetrad description. The equations determined will then form a complete set, which will contain 
as a subset all the 1 + 3 covariant equations just derived (albeit presented in a slightly different form). For 
completeness we will give these equations for a general dissipative relativistic fluid (recent presentations, giv- 
ing the following form of the equations, are B7l l26i ). First we briefly summarize a generic tetrad formalism, 
and then its application to cosmological models (cf. Il30ll48l ). 

3.1 General tetrad formalism 

A tetrad is a set of four mutually orthogonal unit basis vector fields { Cq }a=o.i,2.3, which can be written in 
terms of a local coordinate basis by means of the tetrad components ea^{x^): 



(the latter stating that the i-th component of the a-th tetrad vector is just the directional derivative of the i-th 
coordinate in the direction Cq). This can be thought of as just a general change of vector basis, leading to 
a change of tensor components of the standard tensorial form: T^^^d = e'^i s^'j ej" T^^ ki with obvious 
inverse, where the inverse components e°'i{x^) (note the placing of the indices!) are defined by 






(65) 



6 j — j vv" Ca G- i — Oa 



(66) 
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However, it is a change from an integrable basis to a non-integrable one, so non-tensorial relations (specifi- 
cally: the form of the metric and connection components) are a bit different than when coordinate bases are 
used. A change of one tetrad basis to another will also lead to transformations of the standard tensor form for 
all tensorial quantities: if Cq = Aq" (x*) Cq' is a change of tetrad basis with inverse e^/ = Aa'°'{x'') (each 
of these matrices representing a Lorentz transformation), then T°-''cd — Aa'° A;,''' Ac^ A^;^ T° ^'d'- Again, 
the inverse is obviousFl 

The components of the metric in the tetrad form are given by 

9ab = gij e-a = Ba- Bb = Tjab , (67) 

where r^ab = diag ( — 1, + 1, + 1, + 1 ), showing that the basis vectors are unit vectors mutually orthogonal 
to each other (because the components gab are just the scalar products of these vectors with each other). The 
inverse equation 

=?7afce"<(a;')e^(a;") (68) 

explicitly constructs the coordinate components of the metric from the (inverse) tetrad components e'^i{x^). 
We can raise and lower tetrad indices by use of the metric gab — Vab and its inverse 5"'' = ry"''. 

The commutation functions related to the tetrad are the quantities '~f°'bc{x'^) defined by the commutators 

of the basis vectors 

[ea,eb]^j'ab{x')e, ^ j\c{x') ^ - j\b{x') ■ (69) 
It follows (apply this relation to the local coordinate x') that in terms of the tetrad components, 

7%e(x') = e\ ( eb' d,e,' - e,^ d.Cb' ) = - 2 Cb' ej V[,e'^,] . (70) 

These quantities vanish iff the basis { } is a coordinate basis: that is, there exist local coordinates such 
that ea = (5a' dldx\ iff [ e^, e;, ] = ^ 7"^^ = 0. 

The connection components V°-bc for the tetrad ('Ricci rotation coefficients') are defined by the relations 

Ve.e^ = r^afc ^ T'afc = e", Cb^' V^ Ca* , (71) 

i.e., it is the c-component of the covariant derivative in the 6-direction of the a-vector It follows that all 
covariant derivatives can be written out in tetrad components in a way completely analogous to the usual 
tensor form, for example V aThc = e.a{Thc) ~ ^"^baTdc - ^"^caThd, where for any function /, ea(/) = 
ea^ df /dx"^ is the derivative of / in the direction e^. In particular, because ea{gbc) = for gab = rjab, 
applying this to the metric gives 

Vag&c = ^ -r\a9dc-T''ca9bd = ^ Tj^b), = , (72) 

— the rotation coefficients are skew in their first two indices, when we raise and lower the first indices only. 
We obtain from this and the assumption of vanishing torsion the tetrad relations that are the analogue of the 
usual Christoffel relations: 

T^bc = — (r'^bc — r"cfc) , Tabc = I ( gad l'^ cb — 9bd l'^ ca + gcd l^\b ) ■ (73) 

This shows that the rotation coefficients and the commutation functions are each just linear combinations of 
the other 



The tetrad components of any quantity are invariant when the coordinate basis is changed (for a fixed tetrad), and coordinate 
components are invariant when a change of tetrad basis is made (for a fixed set of coordinates); however, either change will alter the 
tetrad components relative to the given coordinates. 

"Remember that the commutator of any two vectors X,Y is [X,Y] = XY — YX. 
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Any set of vectors whatever must satisfy the Jacobi identities: 

[X, [Y, Z]] + [Y, [Z, X]] + [Z, [X,Y]]=0, 

which follow from the definition of a commutator. Applying this to the basis vectors Ca, ei, and Cc gives the 
identities 

e[a(7'M)+'T'[a6 7''c]e =0, (74) 

which are the integrabihty conditions that the 7" 6c (a^* ) are the commutation functions for the set of vectors . 

If we apply the Ricci identities to the tetrad basis vectors Cq, we obtain the Riemann curvature tensor 
components in the form 

R'^bcd = ec(r"'bd) ~ e-diVbc) + r"ec T'^bd — T^ed T'^bc — ^""be l"^ cd ■ (75) 

Contracting this on a and c, one obtains Einstein's field equations in the form 

Rbd = ea (r°'M) — ed(r'^ba) + T'^ea T'^bd — T^de T'^ba ~ T'bd — \ T g^d + ^ Qbd ■ (76) 

It is not immediately obvious that this is symmetric, but this follows because (|74] | implies Raibcd] = ^ 

Rab = R(ab)- 

3.2 Tetrad formalism in cosmology 

For a cosmological model we choose eo to be the future-directed unit tangent of the matter flow, This 
fixing implies that the initial six-parameter freedom of using Lorentz transformations has been reduced to 
a three-parameter freedom of rotations of the spatial frame { }. The 24 algebraically independent frame 
components of the space-time connection r°f,c can then be split into the set (see Il30ll49ll26l ) 



TaOO = Ua (77) 

TaO/S = I 0<5q/3 + — Eq/J^W^ (78) 

Tq/SO = ^al3-i (79) 

Tap^ = 2 <5/5]^ + e-y5[Q n^]'' + ^ ■ (80) 



The first two sets contain the kinematical variables. In the third is the rate of rotation il" of the spatial frame 
{ } with respect to a Fermi-propagated basis. Finally, the quantities a" and riap = ?^(Q/3) determine the 
9 spatial rotation coefficients. The commutator equations ( [69l ) applied to any space-time scalar / take the 
form 

[eo,e„](/) = u^ea{f)-[\Q5j' + nj + eJ^{LO-' + n-')]epU) (81) 

[e„,e^](/) = 2eaf3iUj'^ eQ{f) + [2a[a5p{' + €apsn^'^ ]e^{f) ] (82) 

The full set of equations for a gravitating fluid can be written as a set of constraints and a set of evolution 
equations, which include the tetrad form of the 1 + 3 covariant equations given above, but complete them by 
giving all Ricci and Jacobi identities for the basis vectors. We now give these equations. 

3.2.1 Constraints 

The following set of relations does not contain any frame derivatives with respect to Cq. Hence, we re- 
fer to these relations as 'constraints'. From the Ricci identities for u° and the Jacobi identities we have the 
(OQ;)-equation (Ci which, in Hamiltonian treatments of the EFE, is also referred to as the 'momentum con- 
straint', the vorticity divergence identity (C2) and the i/afa-equation respectively; the once-contracted 
Bianchi identities yield the (divi?)- and (div iJ) -equations (C4)" and (C5)" ||6ll4Tl; the constraint (Cj)" 
again arises from the Jacobi identities while, finally, {Cg)°'^ and (Cg) stem from the EFE. In detail. 
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+ e"'^^ [{ef3 + 2ufi- ap) {lj^) - nps (J-,^ ] (83) 

= (C2) = (e„-u„-2a,)(L^") (84) 

- e^*'<" [ (e^ - a^) + n/> u;s ] (85) 

= (C5)" = (e^ - 3 (i?"^) + + p) + 3 c^^ (i?"'' - i ^"'') - i n"^ 

+ t'^P^ [ i - a^) (9^) + + \ n^^) - nps H^^ ] (87) 

= {CjT = (e/3-2a^)(n"'^) + f eco" + 2a";3tj'3 + e"'^^[e0(a^)-2L^^f)^] (88) 

= (Cg)"^ = *5"'' + iecr"'3-a<"-yCr'3>^-tJ<"tj'3>+2a;<"f}'3>_(£;a/3_^ (gg^ 

= (Cg) = *i?+|e2-(a„^a"'3) + 2(t^„w")-4(a;„f^")-2^-2A, (90) 
where 

*S'a0 = 6(^(0^3)) + - e'^'^^Q (e|^l - 2 a|^l) (91) 

*i? = 2(2e„-3a„)(a")-i6a" (92) 

6q/3 = 2 np^ -n^^ Uap . (93) 



If ti>" = 0, so that become the normals to a family of 3-spaces of constant time, the last two constraints in 
the set correspond to the symmetric trace-free and trace parts of the Gauss embedding equation ( |54] |. In this 
case, one also speaks of (Cg) as the generalised Friedmann equation, alias the 'Hamiltonian constraint' or 
the 'energy constraint' . 

3.2.2 Evolution of spatial commutation functions 

The 9 spatial commutation functions a" and Uap are generally evolved by equations (40) and (41) given in 
(^1; these originate from the Jacobi identities. Employing each of the constraints (Ci)" to [C^)"-^ listed 
in the previous paragraph, we can eliminate frame derivatives of the kinematical variables Q, aap and 
tj" from their right-hand sides. Thus, we obtain the following equations for the evolution of the spatial 
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connection components: 

eo(a") = -lieS"(3-^a"f3){u>^ + a^) + ^n"0Ujl^-y" 

- \ e"^^ [ {U0 + an) LO^ - nps -{e0 + U0-2 ap) {n^) ] + \ (Ci)« (94) 



e^^<" [ + a.) a/> - K + 2 1^^) n/> ] - | <5"^ (C2) + (Cg)"'^ . (95) 



3.2.3 Evolution of kinematical variables 

The evolution equations for the 9 kinematical variables 6, and are provided by the Ricci identities 
for m", i.e., 

eo(e)-ea(w") = -ie2 + (w„-2aa)w"-(a«0(j"'3) + 2(u;«w") 

- i (/i + 3p) + A (96) 

eoK)- ie"^^e;3(u^) = - § Ow" + tr"/? - i n"/3 li'' - i e"^T [ - 2 0^0;^ ] (97) 

- [E'^f^ - i tt"'') + e^-^^" [ 2 - n^^> ] . (98) 

3.2.4 Evolution of matter and Weyl curvature variables 

Finally, we have the equations for the 4 matter variables /i and (7" and the 10 Weyl curvature variables E^fi 
and -f/a/3> which are obtained from the twice-contracted and once-contracted Bianchi identities, respectively: 

eo(M)+ea(g") = - 6 (/x +p) - 2 («« - a„) - (a„;37r"'^) (99) 
eo(g") + r^e0(p) + e;3(7r"'^) = -| 6 - a"^ / - +p) - («/3 - 3a^) tt^'^ 

- e + \ tt"^) + 3 {ei^^^ - \ ■K^^'^) 

-h \ n^^ H"f' - 3 n^"^, H^^^ - i (2 + a<") g'^^ 
-Fe^''<"[(2u^-a^)7?/> (101) 
+ {uj^ + 2 n^) + i 7r/> ) + i n^'^^ qg ] 

+ eT^<« [ - i 7r/> ) - 2 (102) 

Exercise: (a) Show how most of these equations are the tetrad version of corresponding 1-1-3 covariant 
equations. For which of the tetrad equations is this not true? (b) Explain why there are no equations for 
eo (O" ) and eo (m" ) . [Hint: What freedom is there in choosing the tetrad?] 
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3.3 Complete set 

For a prescribed set of matter equations of state, this gives the complete set of tetrad relations, which can 
be used to characterise particular families of solutions in detail. It clearly contains all the 1 + 3 covariant 
equations above, plus others required to form a complete set. It can be recast into a symmetric hyperbolic 
form 1 26 1 (at least for perfect fluids), showing the hyperbolic nature of the equations and determining their 
characteristics. Detailed studies of exact solutions will need a coordinate system and vector basis, and usually 
it will be advantageous to use tetrads for this purpose, because the tetrad vectors can be chosen in physically 
preferred directions (see lf30l |5P | for the use of tetrads to study locally rotationally symmetric space-times, 
and EH ED forBianchi universes; these cases are both discussed below). 

Finally it is important to note that when tetrad vectors are chosen uniquely in an invariant way (e.g., as 
eigenvectors of a non-degenerate shear tensor or of the electric Weyl curvature tensor), then — because they 
are uniquely defined from 1+3 covariant quantities — all the rotation coefficients above are in fact covariantly 
defined scalars, so all these equations are invariant equations. The only times when it is not possible to define 
unique tetrads in this way is when the space-times are isotropic or locally rotationally symmetric, as discussed 
below. 

4 FLRW models and observational relations 

A particularly important involutive subspace of the space of cosmological space-times is that of the Friedmann- 
Lemaitre ('FL') models, based on the everywhere-isotropic Robertson-Walker ('RW') geometry. It is 
characterised by a perfect fluid matter tensor and the condition that local isotropy holds everywhere: 

^ il'' = aab = OJ^ ^ = Eab = Hab = Xa = Za=WaP, (103) 

the first conditions stating the kinematical quantities are locally isotropic, the second that these models are 
conformally flat, and the third that they are spatially homogeneous. 

Exercise: Show that the implications in this relation follow from the 1 + 3 covariant equations in the 
previous section when p = thus showing that isotropy everywhere implies spatial homogeneity in this 
case. 

4.1 Coordinates and metric 

It follows then that (see ll52l ): 

1. Matter-comoving local coordinates can be founcd so that the metric takes the form 

ds^ = -dt^ + S^{t){dr^ + f{r)dn^) , ^ S^o , (104) 

where dfl"^ = d6^ + sin'^ 6* Ua ~ — Vat, and S/S = 5 S, characterising S{t) as the scale factor for 
distances between any pair of fundamental observers. The expansion of matter depends only on one scale 
length, so it is isotropic (there is no distortion or rotation). 

2. The Ricci tensor ^Rab is isotropic, so the 3-spaces {t ~ const} are 3-spaces of constant (scalar) 
curvature 6k/ S^ where k can be normalised to ± 1, if it is non-zero. Using the geodesic deviation equation 
in these 3-spaces, one finds that (see |52, 53 1) 

/(r) = sinr , r , sinhr if k = +1,0,-1. (105) 

Thus, when fc = + 1, the surface area Att S'^{t) P{r) of a geodesic 2-sphere in these spaces, centred on the 
(arbitrary) point r — Q, increases to a maximum at r = 7r/2 and then decreases to zero again at the antipodal 

'^There are many other coordinate systems in use, for example with different definitions of the radial distance r. 
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point r — tt; hence the point at r = 27r has to be the same point as r = 0, and these 3-spaces are necessarily 
closed, with finite total volume. In the other cases the 3-spaces are usually unbounded, and the surface areas 
of these 2-spaces increase without limit; however, unusual topologies still allow the spatial sections to be 
closed 1541 . 

Exercise: Find the obvious orthonormal tetrad associated with these coordinates, and determine their 
commutators and Ricci rotation coefficients. 

4.2 Dynamical equations 

The remaining non-trivial equations are the energy equation dJTl l. the Raychaudhuri equation ( |29] |. which 
now takes the form 

3| + i(A* + 3p) = 0, (106) 
and the Friedmann equation that follows from ( |55] l: 

^R^2^,-le'^^, (107) 

where A: is a constant. Any two of these equations imply the third if 5* 7^ (the latter equation being a first 
integral of the other two). All one has to do then to determine the dynamics is to solve the Friedmann equa- 
tion. The solution depends on what form is assumed for the matter: Usually it is taken to be a perfect fluid 
with equation of state p = p{fJ,), or as a sum of such fluids, or as a scalar field with given potential V{(f>). For 
the 7-law discussed above, the energy equation integrates to give ( |56] |, which can then be used to represent ^ 
in the Friedmann equation. 

Exercise: Show that on using the tetrad found above, all the other 1 + 3 covariant and tetrad equations are 
identically true when these equations are satisfied. 

4.2.1 Basic parameters 

As well as the parameters Hq, fio; and qo, the FLRW models are characterised by the spatial curvature 
parameter Kq = k/S^ = These parameters are related by the equations ( I6OI 1 and ( |62] i. which are 

now exact rather than approximate relations. 

4.2.2 Singularity and ages 

The existence of the big bang, and age limits on the Universe, follow directly from the Raychaudhuri equation, 
together with the energy assumption (/i + 3p) > (true at least when quantum fields do not dominate), 
because the Universe is expanding today (9o > 0). That is, the singularity theorem above applies in particular 
to FLRW models. Furthermore, from the Raychaudhuri equation, in any FLRW model, the fundamental age 
relation holds (see e.g. 1521 ): 

Age Theorem: In an expanding FLRW model with vanishing cosmological constant and 
satisfying the active gravitational mass density energy condition, ages are strictly constrained by 
the Hubble expansion rate: namely, at every instant, the age to of the model (the time since the 
big bang) is less than the inverse Hubble constant at that time: 

(a* + 3p) > , A = ta< 1/Ho . (108) 

More precise ages io(^Oj ^0) can be determined for any specific cosmological model from the Friedmann 
equation ( |107| l; in particular, in a matter-dominated early universe the same result will hold with a factor 2/3 
on the right-hand side, while in a radiation dominated universe the factor will be 1 /2. Note that this relation 
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applies in the early universe when the expansion rate was much higher, and, hence, shows that the hot early 
epoch ended shortly after the initial singularity ll52l . 

The age limits are one of the central issues in modern cosmology ll55ll4Tl . Hipparchos satellite measure- 
ments suggest a lowering of the age estimates of globular clusters to about 1.2 x 10^ years, together with 
a decrease in the estimate of the Hubble constant to about Hq « 50 km/sec/Mpc. This corresponds to a 
Hubble time 1/Hq of about 1.8 x 10^ years, implying there is no problem, but red giant and cepheid mea- 
surements suggest Hq w 72 — 77 km/sec/Mpc |56|, implying the situation is very tight indeed. However, 
recent supernovae measurements 157 1 suggest a positive cosmological constant, allowing violation of the age 
constraint, and hence easing the situation. All these figures should still be treated with caution; the issue is 
fundamental to the viability of the FLRW models, and still needs resolution. 

4.3 Exact and approximate solutions 

If A = and the energy conditions are satisfied, FLRW models expand forever from a big bang if A: = — 1 
or A: = 0, and recollapse in the future if fc = + 1. A positive value of A gives a much wider choice for 
behaviours 1581 [591 . 

4.3.1 Simplest models 

a) Einstein static model: S{t) — const, fc = +l, A = + 3p) > 0, where everything is constant in space 
and time, and there is no redshift. This model is unstable (see above). 

b) de Sitter model: S{t) — S'unit exp{Ht), H = const. A; = 0, a steady state solution in a constant cur- 
vature space-time: it is empty, because (fi + p) = 0, i.e., it does not contain ordinary matter, but rather a 
cosmological constant^! or a scalar field in the strict 'no-rolling' case. It has ambiguous redshift because the 
choice of families of worldlines and space sections is not unique in this case; see ll60l . 

c) Milne model: S{t) — t, k — — I. This is flat, empty space-time in expanding coordinates (again 

(m + p) = 0). 

d) Einstein-de Sitter model: the simplest non-empty expanding model, with 

fc = = A : S{t) = a i^/^ ^ a = const if p = . 

= 1 is always identically true in this case (this is the critical density case that just manages to expand 
forever). The age of such a model is to = 2/{3Ho); if the cosmological constant vanishes, higher density 
models (Oq > 1) wifl have ages less than this, and lower density models (0 < ilo < 1) ages between this 
value and ( llOSI l. This is the present state of the Universe if the standard inflationary universe theory is correct, 
the high value of il then implying that most of the matter in the Universe is invisible (the 'dark matter' issue; 
see ll4l1 for a summary of ways of estimating the matter content of the Universe, leading to estimates that the 
detected matter in the Universe in fact corresponds to Hq w 0.2 to 0.3). It is thus difficult to reconcile this 
model with observations (the Universe could have flat space sections and a large cosmological constant; but 
then that is not the Einstein-de Sitter model). 



'A fluid with (/^ + p) = is equivalent to a cosmological constant. 
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4.3.2 Parametric solutions 

Use dimensionless conformal time t = J dt/S{t) and rescale S ^ y ^ S{t)/SQ. Then, for a non- 
interacting mixture of pressure-free matter and radiation, we find in the three cases fc = + 1, 0, — 1, 



fc = +l: y= a (1 - cost) +/3 sinr , (109) 

fc = 0: y= aT^/2 + PT, (110) 

fc = — 1: y= a (coshr — 1) + /3 sinhr , (HI) 

where a = H§ 17,„/2, (3 ^ {S^ f^r)^/^, and, on setting i = r = when S ^ 0, 

k = +l: t= Sa[a{T -sinr) + P{1--cost)] , (112) 

fc = 0: t= 5o [arV6 + /3rV2] , (113) 

k^-l: 5o [a(sinliT-T) +/3(coshT- 1) ] . (114) 

It is interesting how in this parametrization the dust and radiation terms decouple; this solution includes as 

special cases the pure dust solutions, /? = 0, and the pure radiation solution, a = 0. The general case 



represents a smooth transition from a radiation dominated early era to a matter dominated later era, and (if 
fc 7^ 0) on to a curvature dominated era, recollapsing if fc = +1. 

4.3.3 Early-time solutions 

At early times, when matter is relativistic or negligible compared with radiation, the equation of state is 
p = i /i and the curvature term can be ignored. The solution is 

/ Q \ 1/4 . 

Sit)=ct'/\ c = const, fi=^t-\T^ ( — ] (115) 

which determines the expansion time scale during nucleosynthesis and so the way the temperature T varies 
with time (and hence determines the element fractions produced), and has no adjustable parameters. Con- 
sequently, the degree of agreement attained between nucleosynthesis theory based on this time scale and 
element abundance observations li61J -[.63 1 may be taken as supporting both a FLRW geometry and the valid- 
ity of the EFE at that epoch. 

The standard thermal history of the hot early Universe (e.g. f6T\) follows; going back in time, the tempera- 
ture rises indefinitely (at least until an inflationary or quantum-dominated epoch occurs), so that the very early 
Universe is an opaque near-equilibrium mixture of elementary particles that combine to form nuclei, atoms, 
and then molecules after pair production ends and the mix cools down as the Universe expands, while various 
forms of radiation (gravitational radiation, neutrinos, electromagnetic radiation) successively decouple and 
travel freely through the Universe that has become transparent to them. This picture is very well supported by 
the detection of the extremely accurate black body spectrum of the CBR, together with the good agreement 
of nucleosynthesis observations with predictions based on the FLRW time scales (II 15l l for the early Universe. 

Exercise: The early Universe was radiation dominated but later became matter dominated (as at the 
present day). Determine at what values S'cqu of the scale factor S{t) matter-radiation equality occurs, as a 
function of fl^. For what values of does this occur before decoupling of matter and radiation? (Note that 
if the Universe is dominated by Cold Dark Matter ('CDM'), then equality of baryon and radiation density 
occurs after this time.) When does the Universe become curvature dominated? 

4.3.4 Scalar field 

The inflationary universe models use many approximations to model a FLRW universe with a scalar field (f> as 
the dominant contribution to the dynamics, so allowing accelerating models that expand quasi-exponentially 
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through many efoldings at a very early time f64",'65l, possibly leading to a very inhomogeneous structure on 
very large (super-particle-horizon) scales |66|. This then leads to important links between particle physics 
and cosmology, and there is a very large literature on this subject. If an inflationary period occurs in the very 
early Universe, the matter and radiation densities drop very close to zero while the inflaton field dominates, 
but is restored during 'reheating' at the end of inflation when the scalar field energy converts to radiation. 

This will not be pursued further here, except to make one point: because the potential V{4)) is unspecified 
(the nature of the inflaton is not known) and the initial value of the 'rolling rate' can be chosen at will, 
it is possible to specify a precise procedure whereby any desired evolutionary history S{t) is attained by 
appropriate choice of the potential V{4)) and the initial 'rolling rate' (see ll67l for details). Thus, inflationary 
models may be adjusted to give essentially any desired results in terms of expansion history. 



4.3.5 Kinetic tlieory 

While a fluid description is used most often, it is also of interest to use a kinetic theory description of the 
matter in the Universe f68l. The details of coUisionless isotropic kinetic models in a FLRW geometry are 
given by Ehlers, Geren and Sachs [69J; this is extended to collisions in |70|. Curiously, it is also possible to 
obtain exact anisotropic coUisionless solutions in FLRW geometries; details are given in IItTI . 



4.4 Phase planes 

From these equations, as well as finding simple exact solutions, one can determine evolutionary phase planes 
for this family of models; see Stabell and Refsdal |59| for {flm, Qo), Ehlers and Rindler [72] for (f2,„, ilr, Qo), 
Wainwright and Ellis |51 1 for {flo, Hq), and Madsen and Ellis |[38l for (fi, S). The latter are based on the 
phase plane equation 

^ = -(37-2)^(1-17). (116) 

This equation is valid for any 7, i.e., for arbitrary relations between /i and p, but gives a (17, S) phase plane 
flow if 7 = 7(11, S), and in particular if 7 = 7(6*) or 7 = const. Non-static solutions can be followed through 
turnaround points where 5* = (and so is infinite). This enables one to attain complete (time-symmetric) 
phase planes for models with and without inflation; see 11381 and 11731 for details. 



4.5 Observations 

Astronomical observations are based on radiation travelling to us on the geodesic null rays that generate our 
past light cone. In the case of a FLRW model, we may consider only radial null rays as these are generic 
(because of spatial homogeneity, we can choose the origin of coordinates on any light ray of interest; because 
of isotropy, light rays travelling in any direction are equivalent to those travelling in any other direction). Thus, 
we may consider geodesic null rays travelling in the FLRW metric ( 1104b such that ds^ = = d9 ~ dcj); then 
it follows that = — dt^ + S'^{t) dr^ on these geodesies. Hence, radiation emitted at E and received at O 
obeys the basic relations 

E Jt, S{t) Js, S{t) Sit) 

yielding the dimensionless matter-comoving radial coordinate distance, where the term S may be found from 
the Friedmann equation ( 1107b . once a suitable matter description has been chosen. 



4.5.1 Redshift 



The first fundamental quantity is redshift. Considering two successive pulses sent from E to O, each remain- 
ing at the same matter-comoving coordinate position, it follows from dl 17| that the cosmological redshift in a 
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FLRW model is given by 

and so directly measures the expansion of the model between when light was emitted and when it is received. 
Two comments are in order. First, redshift is essentially a time-dilation effect, and will be apparent in all 
observations of a source, not just in its spectra; this characterisation has the important consequences that (i) 
redshift is achromatic — the fractional shift in wavelength is independent of wavelength, (ii) the width of 
any emitted frequency band dvE is altered proportional to the redshift when it reaches the observer, i.e., the 
observed width of the band is di'o — {1 + z) dvE, and (iii) the observed rate of emission of radiation and the 
rate of any time variation in its intensity will both also be proportional to (1 + z). Second, there can be local 
gravitational and Doppler contributions zq at the observer, and z^; at the emitter; observations of spectra tell 
us the overall redshift z, given by 

(1 + z) = (l + Zo)(l + Zc)(l + Zi;) , (119) 

but cannot tell us what part is cosmological and what part is due to local effects at the source and the observer. 
The latter can be determined from the CBR anisotropy, but the former can only be estimated by identifying 
cluster members and subtracting off the mean cluster motion. The essential problem is in identifying which 
sources should be considered members of the same cluster. This is the source of the controversies between 
Arp et al and the rest of the observational community (see, e.g.. Field et al Ii74j ). 



4.5.2 Areas 

The second fundamental issue is apparent size. Considering light rays converging to the observer at time 
in a solid angle dO. = sin 9 dO d(f), from the metric form (1104b the coiTesponding null ray^ will be described 
by constant values of 9 and and at the time will encompass an area dA = S'^{tE)P{r)dQ, orthogonal 
to the light rays, where r is given by (II 17b . Thus, on defining the observer area distance ro by the standard 
area relation, we find 

dA = rld^ ^ rl = S^tE)fir). (120) 

Because these models are isotropic about each point, the same distance will relate the observed angle a 
corresponding to a linear length scale £ orthogonal to the light rays: 

i^roa. (121) 

One can now calculate rg from this formula together with ( II 17b and the Friedmann equation, or from the 
geodesic deviation equation (see ll53l ). to obtain for a non-interacting mixture of matter and radiation llTSl . 



ro{z) = 



(go - 1) {1 + 2qoz + qoz^il - - (go - QoPz - 1) 



Hoqoiqo + f3 - 1) {1 + zy 



(122) 



where /3 represents the matter to radiation ratio: (1 — /?) p,„o = 2 /3 pro ■ The standard Mattig relation for 
pressure-free matter is obtained for /3 — 1 1761 . and the corresponding radiation result for /? = 0. 



An important consequence of this relation is refocusing of the past light cone: the Universe as a whole 
acts as a gravitational lens, so that there is a redshift z, such that the observer area distance reaches a maxi- 
mum there and then decreases for larger z; correspondingly, the apparent size of an object of fixed size would 
reach a minimum there and then increase as the object was moved further away |77]. As a specific example, 
in the simplest (Einstein-de Sitter) case with p = A = /c = 0, we find 

/3=l,go = ^ ^ ro(z) = -|-— ^-(\/m-l), (123) 



Bounded by geodesies located at (4>o,6o), {cpo + d.<j), Sq), (ipo^ ^0 + df^), {<t>0 + dcj), Oq + d8). 
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which refocuses at z» = 5/4 llTSl : objects further away will look the same size as much closer objects. For 
example, an object at a redshift zi — 1023 (i.e., at about last scattering) will appear the same angular size 
as an object of identical size at redshift Z2 = 0.0019 (which is very close — it corresponds to a speed of 
recession of about 570 km/sec). In a low density model, refocusing takes place further out, at redshifts up to 
z K, A, depending on the density, and with apparent sizes depending on possible source size evolution 1,791 . 

The predicted (angular size, distance)-relations are difficult to test observationally because objects of 
more or less fixed size (such as spherical galaxies) do not have sharp edges that can be used for measuring 
angular size and so one has rather to measure isophotal diameters (see e.g. llSOl ), while objects with well- 
defined linear dimensions, such as double radio sources, are usually rapidly evolving and so one does not 
know their intrinsic size. Thus, these tests, while in principle clean, are in fact difficult to use in practice. 



4.5.3 Luminosity and reciprocity theorem 

There is a remarkable relation between upgoing and downgoing bundles of null geodesies connecting the 
source at and the observer at to- Define galaxy area distance re as above for observer area distance, 
but for the upgoing rather than downgoing bundle of null geodesies. The expression for this distance will be 
exactly the same as (I120l l except that the times and will be interchanged. Consequently, on using the 
redshift relation (II 18l l. 

Reciprocity Theorem: The observer area distance and galaxy area distance are identical up 
to redshift factors: 

4 - 7^ • (124) 



' G 



[l + zf 



This is true in any space-time as a consequence of the standard first integral of the geodesic deviation equa- 
tion ilsniel. 

Now from photon conservation, the flux of Ught received from a source of luminosity L at time t e will 
be measured to be 

p _ LjtE) 1 1 

47r (1 + z)2 rl ' 



with rQ = S'^{to) P{r) and r given by (II 17b . and the two factors (1 + coming from photon redshift and 
time dilation of the emission rate, respectively. On using the reciprocity result, this becomes 



L{tE) 1 1 
47r (1 + zY 



F = ^ -a - (125) 



where ro is given by (1122b . On taking logarithms, this gives the standard (luminosity, redshift)-relation of 
observational cosmology llTTl . Observations of this Hubble relation basically agree with these predictions, but 
are not accurate enough to distinguish between the various FLRW models. The hopes that this relation would 
determine go from galaxy observations have faded away because of the major problem of source evolution: 
we do not know what the source luminosity would have been at the time of emission. We lack standard 
candles of known luminosity (or equivalently, rigid objects of known linear size, from which apparent size 
measurements would give the answer). Various other distance estimators such as the Tully-Fisher relation 
have helped considerably, but not enough to give a definitive answer. Happily it now seems that Type la 
supernovae may provide the answer in the next decade, because their luminosity can be determined from 
their light curves, which should depend only on local physics rather than their evolutionary history. This is 
an extremely promising development at the present time (see e.g. [57])- 
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4.5.4 Specific intensity 

In practice, we measure (a) in a limited waveband rather than over all wavelengths, as the 'bolometric' 
calculation above suggests; and (b) real detectors measure specific intensity (radiation received per unit solid 
angle) at each point of an image, rather than total source luminosity. Putting these together, we see that if 
the source spectrum is X{ve), i.e., a fraction I{ve) di^E of the source radiation is emitted in the frequency 
range di^E, then the observed specific intensity at each image point is given bj0 

hd, J" I{v{l + z))dv , (126) 
(_i + z) 

where Be is the surface brightness of the emitting object, and the observer area distance tq has canceled 
out (because of the reciprocity theorem). This tells us the apparent intensity of radiation detected in each 
direction — which is independent of (area) distance, and dependent only on the source redshift, spectrum, 
and surface brightness. Together with the angular diameter relation (I121l i. this determines what is actually 
measured by a detector ifSOl . 

An immediate application is black body radiation: if any radiation is emitted as black body radiation at 
temperature Te, it follows (Exercise!) from the black body expression = h{v/TE) that the received 
radiation will also be black body (i.e., have this same black body form) but with a measured temperature of 

^° - (iT^ ■ ^'^'^ 

Note this is true in all cosmologies: the result does not depend on the FLRW symmetries. The importance of 
this, of course, is that it applies to the observed CBR. 



4.5.5 Number counts 



If we observe sources in a given solid angle dVl in a matter-comoving radial coordinate range (r, r + dr), 
the corresponding volume is dV = S^{tE)rQ dr dil, so if the source density isn{tE) and the probability of 
detection is p, the number of sources observed will be 



dN = pn{tE)dV = p 



{l + zf 



S^ito)f{r)drdn 



(128) 



with r given by dl 17t . This is the basic number count relation, where dr can be expressed in terms of 
observable quantities such as dz; the quantity in brackets is constant if source numbers are conserved in a 
FLRW model, that is 

nits) ^ n{to)il + zf . (129) 

The FLRW predictions agree with observations only if we allow for source number and/or luminosity evolu- 
tion (cf. the discussion of spherically symmetric models in the next section); but we have no good theory for 
source evolution. 



The additional problem is that there are many undetectable objects in the sky, including entire galaxies, 
because they lie below the detection threshold; thus we face the problem of dark matter, which is very 
difficult to detect by cosmological observations except by its lensing effects (if it is clustered) and its effects 
on the age of the Universe (if it is smoothly distributed). The current view is that there is indeed such dark 
matter, detected particularly through its dynamical effects in galaxies and clusters of galaxies (see f4T\ for a 
summary), with the present day total matter density most probably in the range 0.1 < f^o < 0.3, while the 
baryon density is of the order of 0.01 < fig'*'^''™'' < 0.03 (from nucleosynthesis arguments). Thus, most of 



Absorption effects will modify this if there is sufficient absorbing matter present; see [6] for relevant formulae. 
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the dark matter is probably non-baryonic. 

To properly deal with source statistics in general, and number counts in particular, one should have a 
reasonably good model of detection limits. It is highly misleading to represent such limits as depending on 
source apparent magnitude alone (see Disney 182]); this does not take into account the possible occurrence 
of low surface brightness galaxies. A useful model based on both source apparent size and magnitude is 
presented in Ellis, Sievers and Perry [83], summarized in [52|. One should note from this particularly that if 
there is an evolution in source size, this has a more important effect on source detectability than an evolution 
in surface brightness. 

Exercise: Explain why an observer in a FLRW model may at late times of its evolution see a situation 
that looks like an island universe (see ll84l ). 

4.6 Observational limits 

The first basic observational limit is that we cannot observe anything outside our past light cone, given 
by ( 1117b . Combined with the finite age of the Universe, this leads to a maximum matter-comoving radial 
coordinate distance from the origin for matter with which we can have had any causal connection: namely 



which converges for any ordinary matter. Matter outside is not visible to us; indeed, we cannot have had any 
causal contact with it. Consequently (see Rindler |85 |), the particles at this matter-comoving coordinate value 
define the particle horizon: they separate that matter which can have had any causal contact with us since 
the origin of the Universe from that which cannot. This is most clearly seen by using Penrose's conformal 
diagrams, obtained on using as coordinates the matter-comoving radius and conformal time; see Penrose ll86l 
and Tipler, Clarke and Ellis [40]. The present day distance to the particle horizon is 



From ( 1117b . this is a sphere corresponding to infinite measured redshift (because S{t) as t — > 0). 

Exercise: Show that once comoving matter has entered the particle horizon, it cannot leave it (i.e., once 
causal contact has been established in a FLRW universe, it cannot cease). 

Actually we cannot even see as far as the particle horizon: on our past light cone information rapidly 
fades with redshift (because of ( 11261 )); and because the early Universe is opaque, we can only see (by means 
of any kind of electromagnetic radiation) to the visual horizon (Ellis and Stoeger ||87l ), which is the sphere 
at matter-comoving radial coordinate distance 



where td is the time of decoupling of matter and radiation, when the Universe became transparent (at about 
a redshift of z = 1100). The matter we see at that time is the matter which emitted the CBR we measure 
today with a present temperature of 2.73 K; its present distance from us is 




(130) 




(131) 




(132) 



Dvh{to) = S{to)ryh ■ 



(133) 



If we evaluate these quantities in an Einstein-de Sitter model, we find an interesting paradox: (re-establishing 
the fundamental constant c,) the present day distance to the particle horizon is Dph{to) = 3cio(= 2c/i/o)- 
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The question is how can this be bigger than cio (the distance corresponding to the age of the Universe). This 
suggests that the matter comprising the particle horizon has been moving away from us at faster than the 
speed of Hght in order to reach that distance. How can this be? To investigate this (ElHs and Rothman |78l), 
note that the proper distance from the origin to a galaxy at matter-comoving radial coordinate r at time t is 
D{t, r) — S{t) r. Its velocity away from us is thus given by Hubble's law 

Q 

v = D = Sr^-D = HD. (134) 

Thus, at any time t, v — c when Dc{t) — c/H{t); this is the speed of light sphere, where galaxies are 
(at the present time) receding away from us at the speed of light; those galaxies at a larger distance will be 
(instantaneously) moving away at a speed greater than c. In the case of an Einstein-de Sitter universe, it 
occurs when Dc = 3cto/2(= c/Hq). This is precisely half the present distance to the particle horizon; the 
latter is thus not the distance where points are moving away from us at the speed of light (however, it is the 
surface of infinite redshift). 

To see that this is compatible with local causality, change from Lagrangian coordinates {t, r) to Eulerian 
coordinates {t, D), where D is the instantaneous proper distance, as above. Then we find the (non-comoving) 
metric form 

ds2 = _ [ 1 _ i^I^ 1 (dctf ^2^DdctdD + dD^ + S'^U) fMt, D)) dfl'^ . (135) 
S"^ S 

It follows that the local light cones are given by 

It is easily seen then that there is no violation of local causality. We also find from this that the past light cone 
of t = to intersects the family of speed of light spheres at its maximum distance from the origin (the place 
where the past light cone starts refocusing), i.e., at 

t* = §fto, SiU)^^Sito), z* = 1.25. (137) 

At that intersection, dD_/dt = (maximum distance!), dDj^ /dt — 2c, so there is no causality violation by 
the matter moving at speed c relative to the central worldline. That matter is presently at a distance cto from 
us. By contrast, the matter comprising the visual horizon was moving away from us at a speed v = 61c when 
it emitted the CBR, and was at a distance of about 10^ light years from our past worldline at that time. Hence, 
it is the fastest moving matter we shall ever see, but was not at the greatest proper distance to which we can 
see (which is D*, see (I137l l). For a full investigation of these matters see 1 



Finally it should be noted that an early inflationary era will move the particle horizon out to very large 
distances, thus ll64l |65]| solving the causal problem presented by the isotropy of CBR arriving here from 
causally disconnected regions (see | 871[78l for the relevant causal diagrams), but it will have no effect on the 
visual horizon. Thus, it changes the causal limitations, but does not affect the visual limits on the part of the 
Universe we can see. 



Exercise: Determine the angular size seen today for the horizon distance Dph{td) at the time of decou- 
pling. What is the physical significance of this distance? How might this relate to CBR anisotropics? 



4.6.1 Small universes 

The existence of visual horizons represent absolute limits on what we can ever know; because of them, we 
can only hope to investigate a small fraction of all the matter in the Universe. Furthermore, they imply we 



4 FLRW MODELS AND OBSERVATIONAL RELATIONS 



31 



do not in fact have the data needed to predict to the future, for at any time gravitational radiation from as 
yet unseen objects (e.g., domain walls in a chaotic inflationary universe) may cross the visual horizon and 
undermine any predictions we may have made |88|. However, there is one exceptional situation: it is possible 
we live in a small universe, with a spatially closed topology on such a length scale (say, 300 to 800 Mpc) 
that we have already seen around the universe many times, thus already having seen all the matter there is 
in the universe. The effect is like being in a room with mirrors on the floor, ceiling, and all walls; images 
from a finite number of objects seem to stretch to infinity. There are many possible topologies, whatever 
the sign of k f541; the observational result — best modelled by considering many identical copies of a basic 
cell attached to each other in an infinitely repeating patterr0 — can be very like the real Universe (Ellis and 
Schreiber 1891 ). In this case we would be able to see our own galaxy many times over, thus being able to 
observationally examine its historical evolution once we had identified which images of distant galaxies were 
in fact repeated images of our own galaxy. 

It is possible the real Universe is like this. Observational tests can be carried out by trying to identify the 
same cluster of galaxies, QSO's ll90l . or X-ray sources in different directions in the sky ||9T1 : or by detecting 
circles of identical temperature variation in the CBR sky (Cornish et al |92| ). If no such circles are detected, 
this will be a reasonably convincing proof that we do not live in such a small universe — which has various 
philosophical advantages over the more conventional models with infinite spatial sections ll88ll . Inter alia they 
give some degree of mixing of CBR modes so giving a potentially powerful explanation of the low degree of 
CBR anisotropy (but this effect is not as strong as some have claimed; see |93l). 

4.7 FLRW universes as cosmological models 

These models are very successful in explaining the major features of the observed Universe — its expansion 
from a hot big bang leading to the observed galactic redshifts and remnant black body radiation, tied in well 
with element abundance predictions and observations (Peebles et al |94J). However, these models do not 
describe the real Universe well in an essential way, in that the highly idealized degree of symmetry does 
not correspond to the lumpy real Universe. Thus, they can serve as basic models giving the largest-scale 
smoothed out features of the observable physical Universe, but one needs to perturb them to get realistic 
('almost-FLRW') Universe models that can be used to examine the inhomogeneities and anisotropics arising 
during structure formation, and that can be compared in detail with observations. This is the topic of the last 
sections. 

However, there is a major underlying issue: because of their high symmetry, these models are infinitely 
improbable in the space of all possible cosmologies. This high symmetry represents a very high degree of fine 
tuning of initial conditions, which is extraordinarily improbable, unless we can show physical reasons why it 
should develop from much more general conditions. In order to examine that question, one needs to look at 
much more general models and see if they do indeed evolve towards the FLRW models because of physical 
processes (this is the chaotic cosmology programme, initiated by Misner 1951 196| , and taken up much later by 
the inflationary universe proposal of Guth |64|). Additionally, while the FLRW models seem good models 
for the observed Universe at the present time, one can ask (a) are they the only possible models that will fit 
the observations? (b) does the Universe necessarily have the same symmetries on very large scales (outside 
the particle horizon), or at very early and/or very late times? 

To study these issues, we need to look at more general models, developing some understanding of their 
geometry and dynamics. This is the topic of the next section. We will find there is a range of models in 
addition to the FLRW models that can fulfill all present day observational requirements. Nevertheless, it is 
important to state that the family of perturbed FLRW models can meet all present observational requirements, 
provided we allow suitable evolution of source properties back in the past. They also provide a powerful the- 
oretical framework for considering the nature of and effects of cosmic evolution. Hence, they are justifiably 

'*In mathematical terms, the universal covering space. 
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the standard models of cosmology. No evidence stands solidly against them. 

Exercise: Apart from the detection of major anisotropy, there are a series of other observations which 
could, if they were ever observed, decisively disprove this family of standard models. What are these obser- 
vations? (See |97| for some.) 



4.8 General observational relations 

Before moving to that section, we briefly consider general observational relations. The present section has 
focussed on the observational relations holding in FLRW models. However, corresponding generic relations 
can be found determining observations in arbitrary cosmologies (see Kristian and Sachs |98| and Ellis |i6j). 
The essential points are as follows. 

In a general model, observations take place on our past light cone, which will develop many cusps and 
caustics at early times because of gravitational lensing, but is still locally generated by geodesic null rays. 

The information we receive comes to us along these null rays, with tangent vector 

A:" = — , kak" = , fca = V,(^ ^ fe'-V^r = , k^Wacj) = . (138) 
dv 

The phase factor <j> determines the local light cone {(f) — const}. Relative to an observer with 4-velocity 
the null vector k"^ determines a redshift factor {—kau"") and a direction Ca'. 

k" = i-kbu'') {u" + e") , eau" = , Cae" = 1 . (139) 

Considering the observed variation (p = u"-'^ a4> of the phase we see that the observed cosmological 
redshift z for comoving matte|3 is given by 

(^ I \ {kaU°')E 

{l + z) = — = . (140) 

\e [kbu'')o 

Taking the derivative of this equation along k'^, we get the fundamental equation |l5]|6l 

^ - - = [ i + i^aen + Kfee'^e") ] dl , (141) 

where dl — {—kau"-) dv. This shows directly the isotropic and anisotropic contributions to redshift from the 
expansion and shear, respectively, and the gravitational redshift contribution from the accelerationHf] In a 
FLRW model, the last two contributions will vanish. 

Area distances are defined as before, and because of the geodesic deviation equation, the reciprocity 
theorem holds unchanged 1610 Consequently, the same surface brightness results as discussed above hold 
generically; specifically, Eq. (|126t holds in any anisotropic or inhomogeneous cosmology. The major dif- 
ference from the isotropic case is that due to the effect of the electric and magnetic Weyl curvatures in the 
geodesic deviation equation, distortions occur in bundles of null geodesies which then cause focusing, re- 
sulting both in strong lensing (multiple images, Einstein rings, and arcs related to cusps and caustics in the 
past light cone) and weak lensing (systematic distortion of images in an observed area); see the lectures by 
Y Mellier and F Bernardeau, the book by Schneider, Ehlers and Falco 1,99 J . and the work by Holz and Wald 

Eool. 



'Cf. the comment on cosmological and local sources of redshift above. 

'^In a static gravitational field, this will be given by an acceleration potential: iia = Va^; see (5]. 

"Because of the first integrals of the geodesic deviation equation; this result can also be shown from use of Liouville's theorem in 
kinetic theory. 



5 SOLUTIONS WITH SYMMETRIES 



33 



Power series equations showing how the kinematical quantities and the electric and magnetic Weyl curva- 
tures affect cosmological observations have been given in a beautiful paper by Kristian and Sachs Ii98il . The 
generalisation of those relations to generic cosmologies has been investigated by Ellis et al BIOIL showing 
how in principle cosmological observations can directly determine the space-time structure on the past null 
cone, and thence off it. Needless to say, major practical observational difficulties make this a formidable task, 
but some progress in this direction is possible (see e.g. Il88l . 01 0211 and 01O3II ). 

Exercise: Apart from area distances, distortions, matter densities, and redshifts, a crucial data set needed 
to completely determine the space-time geometry from the EFE is the transverse velocities of the matter we 
observe on the past light cone IIIOII . Consider how one might try to measure these velocity components, and 
what are the best limits one might place on them by practical measurement techniques. [Hint: One possible 
route is by solar system interferometry. Another is by the Sunyaev-Zel'dovich effect II104I .1 

5 Solutions with symmetries 
5.1 Symmetries of cosmologies 

Symmetries of a space or a space-time (generically, 'space') are transformations of the space into itself 
that leave the metric tensor and all physical and geometrical properties invariant. We deal here only with 
continuous symmetries, characterised by a continuous group of transformations and associated vector fields 

una. 

5.1.1 Killingvectorfields 

A space or space-time symmetry, or isometry, is a transformation that drags the metric along a certain 
congruence of curves into itself. The generating vector field of such curves is called a Killing vector field 
(or 'KV'), and obeys Killing's equations, 

(L^5)y=0 ^ V(,Cj)=0 ^ V,;Cj=-Vje. , (142) 

where is the Lie derivative along By the Ricci identities for a KV, this implies the curvature equation: 

V,;V,Cfc =i?'VCm , (143) 

and so the infinite series of further equations that follows by taking covariant derivatives of this one, e.g., 

ViV,,V,a- - (V,i?'"yfc) + R""^Jk ViU . (144) 

The set of all KV's forms a Lie algebra with a basis { £_a }a=i,2,...,r, of dimension r < ^ n (n — 1). 
denote the components with respect to a local coordinate basis; a, b, c label the KV basis, and k the 
coordinate components. Any KV can be written in terms of this basis, with constant coefficients. Hence: if 
we take the commutator [^a, ^f,] of two of the basis KV's, this is also a KV, and so can be written in terms 
of its components relative to the KV basis, which will be constants. We can write the constants as C^ab, 
obtaining^ 

[S,a, Cb] ^ C'ab^c , C"" be ^ C"" [be] ■ (145) 

By the Jacobi identities for the basis vectors, these structure constants must satisfy 

C\ibC\d] = , (146) 

(which is just equation (|74]l specialized to the case of a set of vectors with constant commutation functions). 
These are the integrability conditions that must be satisfied in order that the Lie algebra exist in a consistent 



^°Ct equation (69). 
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way. The transformations generated by the Lie algebra form a Lie group of the same dimension (see Eisenhart 
ifT^ or Cohn tl06J). 

Arbitrariness of the basis: We can change the basis of KV's in the usual way; 

^a'^Aa'^'ia ^ Ca' ^ ^a'"" C , (147) 

where the A^'" are constants with det (Aq'") ^ 0, so unique inverse matrices A" a exist. Then the structure 
constants transform as tensors: 

C^'a'b' ^A^'cAa'^'Ab^'C^ab. (148) 

Thus the possible equivalence of two Lie algebras is not obvious, as they may be given in quite different 
bases. 

5.1.2 Groups of isometries 

The isometries of a space of dimension n must be a group, as the identity is an isometry, the inverse of an 
isometry is an isometry, and the composition of two isometries is an isometry. Continuous isometries are 
generated by the Lie algebra of KV's. The group structure is determined locally by the Lie algebra, in turn 
characterised by the structure constants 1 106|. The action of the group is characterised by the nature of its 
orbits in space; this is only partially determined by the group structure (indeed the same group can act as a 
space-time symmetry group in quite different ways). 

5.1.3 Dimensionality of groups and orbits 

Most spaces have no KV's, but special spaces (with symmetries) have some. The group action defines orbits 
in the space where it acts, and the dimensionality of these orbits determines the kind of symmetry that is 
present. 

The orbit of a point p is the set of all points into which p can be moved by the action of the isometries of a 
space. Orbits are necessarily homogeneous (all physical quantities are the same at each point). An invariant 
variety is a set of points moved into itself by the group. This will be bigger than (or equal to) all orbits it 
contains. The orbits are necessarily invariant varieties; indeed they are sometimes called minimum invariant 
varieties, because they are the smallest subspaces that are always moved into themselves by all the isometries 
in the group. Fixed points of a group of isometries are those points which are left invariant by the isometries 
(thus the orbit of such a point is just the point itself). These are the points where all KV's vanish (however, the 
derivatives of the KV's there are non-zero; the KV's generate isotropics about these points). General points 
are those where the dimension of the space spanned by the KV's (i.e., the dimension of the orbit through the 
point) takes the value it has almost everywhere; special points are those where it has a lower dimension (e.g., 
fixed points). Consequently, the dimension of the orbits through special points is lower than that of orbits 
through general points. The dimension of the orbit and isotropy group is the same at each point of an orbit, 
because of the equivalence of the group action at all points on each orbit. 

The group is transitive on a surface S (of whatever dimension) if it can move any point of S into any 
other point of S. Orbits are the largest surfaces through each point on which the group is transitive; they 
are therefore sometimes referred to as surfaces of transitivity. We define their dimension as follows, and 
determine limits from the maximal possible initial data for KV's: 

dim surface of transitivity = s, where in a space of dimension n, s < n. 

At each point we can also consider the dimension of the isotropy group (the group of isometries leaving 
that point fixed), generated by all those KV's that vanish at that point: 

dim of isotropy group = q, where in a space of dimension n, q < {n — 1). 
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The dimension r of the group of isometries of a space of dimension nis r = s + q (translations plus 
rotations). From the above limits , < r < n + {n — 1) = {n + 1) (the maximal number of transla- 
tions and of rotations). This shows the Lie algebra of KV's is finite dimensional. 

Maximal dimensions: If ?■ — i n (n + 1), we have a space(-time) of constant curvature (maximal 
symmetry for a space of dimension n). In this case, 

Rijkl = K ( gik Qji - ga g^k ) , (149) 

with K a constant; and K necessarily is a constant if this equation is true and n > 3. One cannot get 

5 = i n (n — 1) — 1 so r i n (n + 1) — 1. 

A group is simply transitive if r = s <^ q ~ (no redundancy: dimensionality of group of isometries 
is just sufficient to move each point in a surface of transitivity into each other point). There is no continuous 
isotropy group. 

A group is multiply transitive if r>s <^ q > (there is redundancy in that the dimension of the 
group of isometries is larger than is needed to move each point in an orbit into each other point). There exist 
non-trivial isotropics. 

5.2 Classification of cosmological symmetries 

We consider non-empty perfect fluid models, i.e., ( fTsT i holds with {fi + p) > 0. 

For a cosmological model, because space-time is 4-dimensional, the possibilities for the dimension of 
the surface of transitivity are s = 0, 1, 2, 3, 4. As to isotropy, we assume + p) ^ 0; then q = 3, 1, or 
because u"- is invariant and so the isotropy group at each point has to be a sub-group of the rotations acting 
orthogonally to u'^ (and there is no 2-dimensional subgroup of 0(3).) The dimension q of the isotropy group 
can vary over the space (but not over an orbit): it can be greater at special points (e.g., an axis centre of 
symmetry) where the dimension s of the orbit is less, but r (the dimension of the total symmetry group) must 
stay the same everywhere. Thus the possibilities for isotropy at a general point are: 

a) Isotropic: q = 3, the Weyl curvature tensor vanishes, kinematical quantities vanish except Q. All 
observations (at every point) are isotropic. This is the FLRW family of space-time geometries; 

b) Local Rotational Symmetry ('LRS'): q — 1, the Weyl curvature tensor is of algebraic Petrov type 
D, kinematical quantities are rotationally symmetric about a preferred spatial direction. All observations at 
every general point are rotationally symmetric about this direction. All metrics are known in the case of dust 
(30] and a perfect fluid (see [50] and also 11071 ). 

c) Anisotropic: q ~ 0; there are no rotational symmetries. Observations in each direction are different 
from observations in each other direction. 

Putting this together with the possibilities for the dimensions of the surfaces of transitivity, we have the 
following possibilities (see Figure 1): 

5.2.1 Space-time homogeneous models 

These models with s — 4 are unchanging in space and time, hence is a constant, so by the energy con- 
servation equation (|37] | they cannot expand: = 0. Thus by (1140b they cannot produce an almost isotropic 
redshift, and are not useful as models of the real Universe. Nevertheless, they are of some interest. 

The isotropic case q — 3 (=> r = 7) is the Einstein static universe, the non-expanding FLRW model 
(briefly mentioned above) that was the first relativistic cosmological model found. It is not a viable cosmol- 
ogy inter alia because it has no redshifts, but it laid the foundation for the discovery of the expanding FLRW 



5 SOLUTIONS WITH SYMMETRIES 



Dim invariant variety 

Dimension 



Isotropy 
Group 


8 = 2 

inhomogeneous 


s = 3 

spatially 
homogeneous 


s = 4 

space-time 
homogeneous 


q = 

aniso- 
tropic 


generic metric form known 
Spatially self-similar, 
Abelian G_2 on 2-d 

spacelike surfaces, 
non-Abelian G_2 


Bianchi : 
orthogonal, 
tilted 


Osvath/Kerr 


q = 1 

LRS 


Lemaitre-Tolman- 
Bondi family 


Kant owski- Sachs , 
LRS Bianchi 


G"odel 


q = 3 

isotropic 


none 

(cannot happen) 


Friedmann 


Einstein 
static 




two non-ignorable 
coordinates 


one non-ignorable 
coordinate 


algebraic EFE 
(no redshift) 





Dim invariant variety 

Dimension 

Isotropy s = s = 1 

Group 

inhomogeneous inhomogeneous /no isotropy group 



q = Szekeres-Szaf ron. General metric 

Stephani-Barnes, form independent 

Oleson type N of one coord; 

KV h . s . o . / not h . s . o . 

The real universe! 



Figure 1: Classification of cosmological models (with (n+p) > 0) by isotropy and homogeneity. 
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models. 

Exercise: What other features make this space-time problematic as a cosmological model? 

The LRS case q = 1 r = 5) is the Godel stationary rotating universe [1081, also with no redshifts. 
This model was important because of the new understanding it brought as to the nature of time in General 
Relativity (see |I3]|401l109|). Inter alia, it is a model in which causality is violated (there exist closed timelike 
curves through each space-time point) and there exists no cosmic time function whatsoever 

The anisotropic models q — (^ r = 4) are all known, 1,1 10] , but are interesting only for the light they 
shed on Mach's principle; see 111 1 111 . 

5.2.2 Spatially homogeneous universes 

These models with s — 3 are the major models of theoretical cosmology, because they express mathematically 
the idea of the 'cosmological principle' : all points of space at the same time are equivalent to each other 111 121 . 

The isotropic case q = 3 (^ r = 6) is the family of FLRW models, the standard models of cosmology 
discussed above that have the matter-comoving metric form ( 1104b . 

The LRS case q = 1 (=> r = 4) is the family of Kantowski-Sachs universes 111 1311 - 111 151 plus the LRS 
orthogonal (49] and LRS tilted 11161 Bianchi models. The simplest are the Kantowski-Sachs family, with 
matter-comoving metric form 

ds^ ^-de + A^{t) dr^ + B\t){de^ + f{e) d(f>^ ) , (150) 

where f{6) is given by ( 11051 ). 

The anisotropic case q = r = 3) is the family of orthogonal and tilted Bianchi models with a 
group of isometrics G3 acting simply transitively on spacelike surfaces. The simplest class is the Bianchi 
Type I family, discussed later in this section. The family as a whole has quite complex properties; these 
models are discussed in the following section. 

5.2.3 Spatially inhomogeneous universes 

These models have s < 2. 

The LRS cases ((7 = 1 s = 2, r = 3) are the spherically symmetric family with matter-comoving 
metric form 

ds^ = - C^{t, r) di^ + A^{t, r) dr^ + B^{t, r) ( dO^ + f{d) dcj)^ ) , (151) 

where f{9) is given by ( I105l l. In the dust case, we can set C(t, r) = 1 and can integrate Einstein's field 
equations analytically; for /c = + 1, these are the spherically symmetric Lemaitre-Tolman-Bondi mod- 
els ('LTB') 11 171 - 11 1911 . discussed later in this section. They may have a centre of symmetry (a timelike 
worldline), and can even allow two such centres, but they cannot be isotropic about a general point (because 
isotropy everywhere implies spatial homogeneity; see the discussion of FLRW models). 

The anisotropic cases (g = s < 2, r < 2) include solutions admitting an Abelian or non-Abelian 
group of isometrics G2, the G2 cosmologies, and spatially self-similar models (see e.g. lISTl ). 

Solutions with no symmetries at all have r = 0^s = 0, g^O. The real Universe, of course, belongs to 
this class; all the others are intended as approximations to this unique Universe. Remarkably, we know some 
exact solutions without symmetries, specifically (a) Szekeres' quasi-spherical models 11120111211 . that are 
in a sense non-linear FLRW perturbations Ill22i . with matter-comoving metric form 

ds^ = - dt^ + e^^ dx^ + e^^ {dy'^ + dz^) , A = A{t,x,y, z) , B = B{t,x,y, z) , (152) 

(b) Stephani's conformally flat models II123II124L and (c) Oleson's type N solutions (for a discussion of 
these and all the other inhomogeneous models, see Krasiriski fS] and Kramer et al [1251). One further inter- 
esting family without global symmetries are the Swiss-Cheese models made by cutting and pasting segments 
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of spherically symmetric models. These are discussed below. 

We now discuss the simplest useful anisotropic and inhomogeneous models, before turning to the Bianchi 
models in the next section. 

5.3 Bianchi Type I universes {s = 3) 

These are the simplest anisotropically expanding Universe models. The metric can be given in matter- 
comoving coordinates in the form Ml 261 

ds^ = -df + X^{t)dx^ + Y'\t)dy^ + Z^{t)dz^, u°=(5"o- (153) 

This is the simplest generalisation of the spatially flat FLRW models to allow for different expansion factors 
in three orthogonal directions; the corresponding average expansion scale factor is S{t) = {XYZY^^. 
They are spatially homogeneous, being invariant under an Abelian group of isometries Ga, simply transitive 
on spacelike surfaces {t = const}, so s — 3; in general g = =J> r = 3, but there are LRS and isotropic 
subcases (the latter being the Einstein-de Sitter universe). The space sections {t = const} are flat (when 
t = to, all the metric coefficients are constant), and all invariants depend only on the time coordinate t. The 
fluid flow (orthogonal to these homogeneous surfaces) is necessarily geodesic and irrotational. Thus these 
models obey the restrictions 

0^U-^CJ\ = Xa = Za=VaP, = ^i?ah ■ (154) 

The 1 + 3 covariant equations obeyed by these models follow from the 1 + 3 covariant equations in subsection 
l2.2l on making these restrictions. We can find a tetrad in the obvious way from the above coordinates (ei ' = 
X{t)^^ Si \ etc.); then the tetrad equations of the subsection l3.2l hold with 

= w"=w"=r2", = a"=rv/3, = e„(e) = e„(cr/37) , = e„(/^) = e„(p) . (155) 



It follows that the (Oa) -equation (|32] |. which is (Ci)" in the tetrad form, is identically satisfied, and also that 
Hab = and WbE"^'' = 0. From the Gauss embedding equation (|54] |. the shear obeys 



(5Va6)-=0 ^ aab = ^, (Sa6)-=0, (156) 

which implies 

^t' = |^, E2 = iSabE''\ (E2)-=0. (157) 

All of Einstein's field equations will then be satisfied if the conservation equation (l37l i. the Raychaudhuri 
equation ( |29] l. and the Friedmann-like equation (ISST i are satisfied. As in the FLRW case, the latter is the first 
integral of the other two. Assuming a 7-law equation of state, (|56] | will be satisfied and, using ( 1157b . equation 
i55[ becomes the generalised Friedmann equation, 

Y? M 

3^ = ^ + ^. (158) 

This shows that no matter how small the shear today, it will (for ordinary matter) dominate the very early 
evolution of the Universe model, which will then approximate Kasner's vacuum solution II125I . 

On writing out the tetrad components of the shear equation ( |156t , using the commutator relations (ISTT i to 
determine the shear components, one finds that the individual length scales are given by. 



X{t) = S{t) exp(Ei W{t)) , Y{t) = S{t) cxp(E2 W{t)) , Z{t) = S{t) exp(S3 W{t)) , 
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where 

^W-/^' ''''' 

and the constants Eq, satisfy 

El + S2 + S3 = , E2 + E^ + e2 . 
These relations can be satisfied by setting 

27r 47r 

Eq = (2E/3) sin aa , ai = a , a2 = a + — , = a + — , (160) 

o o 

and a is a constant. Thus, the solution is given by choosing a value for 7, and then integrating successively 
(fBSl l and ( fT59] l. 

Exercise: Show that, on using the obvious tetrad associated with the coordinates above, all the tetrad (and 
1 + 3 covariant) equations are then satisfied. 

For example, in the case of dust (7 = 1) we have: 



so 



Xit) = S{t) ^) , ni) = ^WU77vT ' Z{t) = Sit)l 



The generic case is anisotropic; LRS cases occur when a — tt/Q and a — 7r/2 in (11601 ). and isotropy when 
E = 0. 

At late times this isotropizes to give the Einstein-de Sitter model, and, hence, as mentioned above, can be 
a good model of the real Universe if E is chosen appropriately. However, at early times, the situation is quite 
different. As i ^ 0, provided E 7^ 0, then S{t) (VSE)!/^ ^1/3 ^^^j 

Plotting the function f{a) = f (5 + sin a), we see that the generic behaviour occurs for a 7^ tt/2; in this 
case two of the powers are positive but one is negative, so going backwards in time, the collapse along the 
preferred axis reverses and changes to a (divergent) expansion, while collapse continues (divergently) along 
the two orthogonal direction; the singularity is a cigar singularity. Going forward in time, a collapse along 
the preferred axis stops and reverses to become an expansion. However when a = 7r/2, one exponent is pos- 
itive but the other two are zero. Hence, going back in time, collapse continues divergently along the preferred 
direction in these LRS solutions back to the singularity, but in the orthogonal directions it slows down and 
halts; this is a pancake singularity. An important consequence in this special case is that particle horizons 
are broken in the preferred direction — communication is possible to arbitrary distance in a cylinder around 
this axis IH. 



One can work out detailed observational relations in these models. Because of the high symmetry, the 
null geodesies can be found explicitly; those along the three preferred axes are particularly simple. Redshift 
along each of these axes simply scales with the expansion ratio in that direction. Area distances can be found 
explicitly M127II128I . An interesting feature is that all observations will show an eight-fold discrete isotropy 
symmetry about the preferred axes. One can also work out helium production and CBR anisotropy in these 
models, following the pioneering paper by Thorne Iil29j . Because the shear can dominate the dynamics at 
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nucleosynthesis or baryosynthesis time, causing a speeding up of the expansion, one can get quite differ- 
ent resuhs than in the FLRW models. Consequently, one can use the nucleosynthesis observations to limit 
the shear constant E, but still allowing extra freedom at the time of baryosynthesis. The CBR quadrupole 
anisotropy will directly measure the difference in expansion along the three principal axes since last scatter- 
ing, and, hence, may also be used to limit the anisotropy parameter E. Nucleosynthesis gives stronger limits, 
because it probes to earlier times. These models have also been investigated in the case of viscous fluids and 
kinetic theory solutions (Misner 11301 ). with electromagnetic fields, and also the effects of 'reheating' on the 
CBR anisotropy and spectrum have been examined; see Rees 01311 . 

Thus, these models can have arbitrarily small shear at the present day, thus can be arbitrarily close to an 
Einstein-de Sitter universe since decoupling, but can be quite different early on. 

Exercise: Show how the solutions will be altered by (i) a fluid with simple viscosity: iTab — — iT^ab with 
constant viscocity coefficient 77, (ii) freely propagating neutrinos 1130|. 



5.4 Lemaitre-Tolman-Bondi family (s = 2) 

The simplest inhomogeneous models are those that are spherically symmetric. In general they are time- 
dependent, with 2-dimensional spherical surfaces of symmetry: s = 2, g = l=>r = 3. The geometry of 
this family (including the closely related models with flat and negatively curved 2-surfaces of symmetry), is 
examined in a 1 + 3 co variant way by van Elst and Ellis 01O7I . and a tetrad analysis is given by Ellis ll30ll 
(the pressure-free case) and Stewart and Ellis fSO"] (for perfect fluids). Here we only consider the dust case, 
because then a simple analytic solution is possible; the perfect fluid case includes spherical stellar models 
and collapse solutions (see, e.g., Misner, Thorne and Wheeler B132I ). 



The general spherically symmetric metric for an irrotational dust matter source in synchronous matter- 
comoving coordinates is the Lemaitre-Tolman-Bondi ('LTB') metric 111 171 - 111 191 

ds"^ = -dt"^ + X^{t,r)dr'^ +Y'^{t,r)dn'^ , = S^o ■ (161) 

The function Y ~ Y{t, r) is the areal radius, since the proper area of a sphere of coordinate radius r on a 
time slice of constant t is AttY'^ (upon re-establishing factors of tt). Solving Einstein's field equations [11 1911 
shows 

1 + 2,hj[r) 

wherey(t,r) = dY{t,r)/dr, md dfP = dO^ + sir? 9 d(jP',yN\\hY{t,r) obeying a generaUsed Friedmann 
equation, 

n^:0=±l/^^ + 2i?(r) , (163) 



and the energy density given by 



Equation ( |163l l can be solved in terms of a parameter -q — ri{t, r): 
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when 
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coshry— 1, ( sinhry — 77. 

£{r) = { 1, 00 = <( (l/2)?7^ C = < (l/6)r/^ when { E = , (166) 

1 — COST], [ — sin 77, 

for hyperboUc, parabolic and elliptic solutions, respectively. 





The LTB model is characterised by three arbitrary functions of the matter-comoving coordinate radius r. 
E = E{r) > — 1 has a geometrical role, determining the local 'embedding angle' of spatial slices, and also 
a dynamical role, determining the local energy per unit mass of the dust particles, and, hence, the type of 
evolution of y. M = M{r) is the effective gravitational mass within coordinate radius r. ts — tsir) is the 
local time at which y = 0, i.e., the local time of the big bang — we have a non-simultaneous bang surface. 
Specification of these three arbitrary functions — M(r), E{r) and tsir) — fully determines the model, and 
whilst all have some type of physical or geometrical interpretation, they admit a freedom to choose the radial 
coordinate, leaving two physically meaningful choices, e.g., r = r(M), E = E{M), ts — tsiM). For 
particular choices of this initial data, one obtains FLRW models, which, of course, are special cases of these 
spherical models with very specific initial data. In fact, the FLRW models are obtained if one sets 

2E[r) = -kr^, Y{t,r) = S{t) r , M{r) ^ — ^i{t)Y^ . (167) 
The LTB models have been used in a number of interesting ways in cosmology: 

* to give simple models of structure formation II133II1341 . e.g., by looking at evolution of a locally open 
region in a closed universe 1 135 1 and evolution of density contrast II136L 

* to give Universe models that are inhomogeneous on a cosmological scale Ill37lll38ll . 

* to examine inhomogeneous big bang structures II139I . 

* to examine CBR anisotropics 1 1 401- 111421 . 

* to investigate observational conditions for spatial homogeneity II143II - M1471 . 

* to trace the effect of averaging on spatial inhomogeneities |1481, and 

* to look at the relationship between cosmic evolution and closure of the Universe II 1491 . 

These aspects are discussed in Krasinski's book 1 8 j. Part 111. Here we will only summarize one interesting 
result: namely, regarding observational tests of whether the real Universe is more like a LTB inhomogeneous 
model, or a FLRW model. In Mustapha, Hellaby and Ellis II144L the following result is shown: 

Isotropic Observations Tlieorem (1): Any given isotropic set of source observations n{z) 
and m(z), together with any given source luminosity and number evolution functions L{z) and 
N{z), can be fitted by a spherically symmetric dust cosmology — a LTB model — in which ob- 
servations are spherically symmetric about us because we are located near the central worldline. 

This shows that any spherically symmetric observations we may eventually make can be accommodated 
by appropriate inhomogeneities in a LTB model — irrespective of what source evolution may occur. In par- 
ticular, one can find such a model that will fit the observations if there is no source evolution. The following 
result also holds: 



Isotropic Observations Theorem (2): Given any spherically symmetric geometry and any 
spherically symmetric set of observations, we can find evolution functions that will make the 
model compatible with the observations. This applies in particular if we want to fit observations 
to a FLRW model. 

Strictly speaking, the hyperbolic, parabolic and elliptic solutions obtain when YE/M > 0, = and < 0, respectively, since E = 
at a spherical origin in both hyperbolic and elliptic models. 
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The point of the first result is that it shows that these models — spherical inhomogeneous generalisations 
of the FLRW models — are viable models of the real Universe since decoupling, because they cannot be 
observationally disproved (at least not in any simple way). The usual response to this is: but the FLRW 
models are confirmed by observations. The second result clarifies this: yes they are, provided you allow an 
evolution function to be chosen specifically so that the initially discrepant number counts fit the FLRW model 
predictions. That is, we assume the FLRW geometry, and then determine what source evolution makes this 
assumption compatible with observations II150I . Without this freedom, the FLRW models are contradicted 
by, e.g., radio source observations, which without evolution are better fitted by a spatially flat (Euclidean) 
model. Thus, the FLRW model fit is obtained only because of this freedom, allowed because we do not 
understand source luminosity and number evolution. The inhomogeneous LTB models provide an alternative 
understanding of the data; the observations do not contradict them. 

Exercise: Suppose (a) observations are isotropic, and (b) we knew the source evolution function and se- 
lection function, and (c) we were to observationally show that after taking them into account, the observer 
area distance relation has precisely the FLRW form (I122l l for /3 = 1 and the number count relation implies 
the FLRW form (|129t . Assuming the space-time matter content is dust, (i) prove from this that the space-time 
is a FLRW space-time. Now (ii) explain the observational difficulties that prevent us using this exact result 
to prove spatial homogeneity in practice. 

An alternative approach to proving homogeneity is via the Postulate of Uniform Thermal Histories 
('PUTH') — i.e., the assumption that because we see similar kinds of objects at great distances and nearby, 
they must have had similar thermal histories. One might then hope that from this one could deduce spatial 
homogeneity of the space-time geometry (for otherwise the thermal histories would have been different). 
Unfortunately, the argument here is not watertight, as can be shown by a counterexample based on the LTB 
models (Bonnor and Ellis 1 143 1). Proving — rather than assuming — spatial homogeneity remains elusive. 
We cannot observationally disprove spatial inhomogeneity. However, we can give a solid argument for it via 
the EGS theorem discussed below. 

5.5 Swiss-Cheese models 

Finally, an interesting family of inhomogeneous models is the Swiss-Cheese family of models, obtained by 
repeatedly cutting out a spherical region from a FLRW model and filling it in with another spherical model: 
Schwarzschild or LTB, for example. This requires: 

(i) locating the 3-dimensional timelike junction surfaces Sj- in each of the two models; 

(ii) defining a proposed identification $ between and S_; 

(iii) determining the junction conditions that (a) the 3-dimensional metrics of and I]_ (the first 
fundamental forms of these surfaces) be isometric under this identification, so that there be no disconti- 
nuity when we glue them together — we arrive at the same metric from both sides — and (b) the second 
fundamental forms of these surfaces (i.e., the covariant first derivatives of the 3-dimensional metrics along 
the spacelike normal directions) must also be isometric when we make this identification, so that they too 
are continuous in the resultant space-time — equivalently, there is no discontinuity in the direction of the 
spacelike unit normal vector as we cross the junction surface S (this is the condition that there be no surface 
layer on E once we make the join; see Israel 1 151]). 

Satisfying these junction conditions involves deciding how the 3-dimensional junction surfaces I]± should 
be placed in the respective background space-times. It follows from them that 4 of the 10 components of Tab 
must be continuous: if n° denotes the spacelike (or, in some other matching problems, timelike) unit nor- 
mal to S and pj' the tensor projecting orthogonal to n°, then {Tabn""n^) and T^c n'^Pa'^ must be continuous, 
but the other 6 components TdPa'^Pb^ can be discontinuous (thus at the surface of a star, in which case 
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ri° will be spacelike, the pressure is continuous but the energy density can be discontinuous). Conserva- 
tion of the energy-momentum tensor across the junction surface S will then be satisfied by the constraints 

= (Gab - Tab) n'^nb and = (G,, - T^,) n^p,-. 

(iv) Having determined that these junction conditions can be satisfied for some particular identification 
of points, one can then proceed to identify these corresponding points in the two surfaces and thus 
gluing an interior Schwarzschild part to an exterior FLRW part, for example. Because of the reciprocal na- 
ture of the junction conditions, it is then clear we could have joined them the other way also, obtaining a 
well-matched FLRW interior and Schwarzschild exterior. 

(v) One can continue in this way, obtaining a family of holes of different sizes in a FLRW model with 
different interior fillings, with further FLRW model segments fitted into the interiors of some of these re- 
gions, obtaining a Swiss-Cheese model. One can even obtain a hierarchically structured family of spherically 
symmetric vacuum and non-vacuum regions in this way. 

It is important to note that one cannot match arbitrary masses. It follows from the junction conditions that 
the Schwarzschild mass in the interior of a combined FLRW-Schwarzschild solution must be the same as the 
mass that has been removed: Mschw = (47r/3) (/i ?'^)flrw- If the masses were wrongly matched, there 
would be an excess gravitational field from the mass in the interior that would not fit the exterior gravitational 
field, and the result would be to distort the FLRW geometry in the exterior region — which then would no 
longer be a FLRW model. Alternatively viewed, the reason this matching of masses is needed is that oth- 
erwise we will have fitted the wrong background geometry to the inhomogeneous Swiss-Cheese model — 
averaging the masses in that model will not give the correct background average 1 152] , and they could not 
have arisen from rearranging uniformly distributed masses in an inhomogeneous way (this is the content of 
Traschen's integral constraints fl53l ). Consequently, there can be no long-range effects of such matching: 
the Schwarzschild mass cannot cause large-scale motions of matter in the FLRW region. 

These models were originally developed by Einstein and Straus [1541 (see also Schiicking P1551) to ex- 
amine the effect of the expansion of the Universe on the solar system (can we measure the expansion of the 
Universe by laser ranging within the solar system?). Their matching of a Schwarzschild interior to a FLRW 
exterior showed that this expansion has no effect on the motion of planets in the Schwarzschild region. It 
does not, however, answer the question as to where the boundary between the regions should be placed — 
which determines which regions are affected by the universal expansion. Subsequent uses of these models 
have included: 

* examining Oppenheimer-Snyder collapse in an expanding universe II156I - II158I . 

* examining gravitational lensing effects on area distances II159I . 

* investigating CBR anisotropics 1 160I- I162J . 

* modelling voids in large-scale structure II163II1641 . perhaps using surface-layers B165I . 

* modelling the Universe as a patchwork of domains of different curvature k = 0, ± 1 [166]. 

Exercise: Show how appropriate choice of initial data in a LTB model can give an effective Swiss-Cheese 
model with one centre surrounded by a series of successive FLRW and non-FLRW spherical regions. Can 
you include (i) flat, (ii) vacuum (Schwarzschild) regions in this construction? 

One of the most intriguing questions is what non-spherically symmetric models can be joined regularly 
onto a FLRW model. Bonnor has shown that some Szekeres anisotropic and inhomogeneous models can be 
matched to a dust FLRW model across a matter-comoving spherical junction surface II167I . Dyer et al 01681 
have shown that one can match FLRW and LRS Kasner (anisotropic vacuum Bianchi Type I) models across 
a flat 3-dimensional timelike junction surface. Optical properties of these Cheese-slice models have been 
investigated in depth 111691 . 
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6 Bianchi models 

These are the models in which there is a group of isometries G3 acting simply transitively on spacelike sur- 
faces {t = const}, so they are spatially homogeneous. There is only one essential dynamical coordinate (the 
time t), and Einstein's field equations reduce to ordinary differential equations, because the inhomogeneous 
degrees of freedom have been 'frozen out'. They are thus quite special in geometrical terms; nevertheless, 
they form a rich set of models where one can study the exact dynamics of the full non-linear field equations. 
The solutions to the EFE will depend on the matter in the space-time. In the case of a fluid (with uniquely 
defined flow lines), we have two different kinds of models: 

Orthogonal models, with the fluid flow lines orthogonal to the surfaces of homogeneity (Ellis and Mac- 
Callum |49|, see also lISTl ): 

Tilted models, with the fluid flow lines not orthogonal to the surfaces of homogeneity; the components 
of the fluid's peculiar velocity enter as further variables (King and Ellis 1.1 16J . see also [170] ). 

Rotating models must be tilted ( cf. Eq. dZTl i ), and are much more complex than non-rotating models. 

6.1 Constructing Bianchi models 

There are essentially three direct ways of constructing the orthogonal models, all based on properties of a 
tetrad of vectors { Bq } that commute with the basis of KV's { }, and usually with the timelike basis vector 
chosen parallel to the unit normal Ua = ~ Vat to the surfaces of homogeneity, i.e., Cq = n. 

first approach (Taub II171I . Heckmann and Schiicking II126I ) puts all the time variation in the metric 
components: 

ds^ ^-dt^+ -icp{t) (e",(x^) dx') {e^j{x^) dx^) , (168) 

where e°'i{x'^) are 1-forms inverse to the spatial triad vectors ea'^{x''), which have the same commutators 
C"- fj^, a, (3,j, - ■ ■ = 1, 2, 3, as the structure constants of the group of isometries and commute with the unit 
normal vector Cq to the surfaces of homogeneity; i.e., Gq, = Eq,' d/dx^, Bq = d/dt obey the commutator 
relations 

[ea-, e(3]= ape^ ■, [eo, e^J^O, (169) 

where the af3 are the Lie algebra structure constants satisfying the Jacobi identities (I146l l. Einstein's field 
equations (H) become ordinary differential equations for the metric functions "fapit)- 

The second approach is based on use of the automorphism group of the isometry group with time- 
dependent parameters. We will not consider it further here (see Collins and Hawking il72i . Jantzen M173II1741 
and Wainwright and Ellis [51] for a discussion). 

The third approach (Ellis and MacCallum fA9\), which is in our view the preferable one, uses an or- 
thonormal tetrad based on the normals to the surfaces of homogeneity (i.e., Gq = n, the unit normal vector 
to these surfaces). The tetrad is chosen to be invariant under the group of isometries, i.e., the tetrad vec- 
tors commute with the KV's, and the metric components in the tetrad are space-time constants, cjab — Vat', 
now the dynamical time variation is in the commutation functions for the basis vectors, which then deter- 
mine the time-(and space-)dependence in the basis vectors themselves. Thus, we have an orthonormal basis 
{ Ga }a=o, 1,2,3, such that 

[ea,eb]=-/'ab{t)ec. (170) 

The commutation functions 'y"'bc{t), together with the matter variables, are then treated as the dynamical 
variables. Einstein's field equations ([U are first-order equations for these quantities, supplemented by the 
Jacobi identities for the -f'^bcit), which are also first-order equations. It is sometimes useful to introduce also 
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the Weyl curvature components as auxiHary variables, but this is not necessary in order to obtain solutions. 
Thus the equations needed are just the tetrad equations given in section[3] specialised to the case 

-0 = e„(7%,) . (171) 

The spatial commutation functions '^"^^[t) can be decomposed into a time-dependent matrix napit) and a 
vector a"(t) (see (IHOll). and are equivalent to the structure constants C"/}^ of the symmetry group at each 
point|3 In view of ( 1171b . the Jacobi identities dSSl l now take the simple form 

n'^l^af} = 0. (172) 

The tetrad basis can be chosen to diagonalise Uap to attain Uap — diag {ni, n2, n^) and to set a" — (a, 0, 0), 
so that the Jacobi identities are then simply nia — 0. Consequently, we define two major classes of structure 
constants (and so Lie algebras): 

Class A: a — , 
Class B: a ^ . 

Following Schiicking, the adaptation of the Bianchi classification of G3 isometry group types used is as in 
Figure 2. Given a specific group type at one instant, this type will be preserved by the evolution equations for 
the quantities na{t) and a{t). This is a consequence of a generic property of Einstein's field equations: they 
will always preserve symmetries in initial data (within the Cauchy development of that data); see Hawking 
andElhs [iSj. 

In some cases, the Bianchi groups allow higher symmetry subcases: isotropic (FLRW) or LRS models. 
Figure 3 gives the Bianchi isometry groups admitted by FLRW and LRS solutions li49J , i.e., these are the 
simply transitive 3-dimensional subgroups allowed by the full Gg of isometrics (in the FLRW case) and the 
G4 of isometrics (in the LRS case). The only LRS models not allowing a simply transitive subgroup G3 are 
the Kantowski-Sachs models for fc = 1. 

Tilted models can be constructed similarly, as discussed below, or by using non-orthogonal bases in 
various ways [11161 : those possibilities will not be pursued further here. 

6.2 Dynamics of Bianchi models 

The set of tetrad equations (section 3) with these restrictions will determine the evolution of all the com- 
mutation functions and matter variables, and, hence, determine the metric and also the evolution of the Weyl 
curvature (these are regarded as auxiliary variables). In the case of orthogonal models — the fluid 4-velocity 
u is parallel to the normal vectors n — the matter variables will be just the fluid density and pressure [49j; in 
the case of tilted models — the fluid 4-velocity u is not parallel to the normal vectors n — we also need the 
peculiar velocity of the fluid relative to the normal vectors fl 161, determining the fluid energy-momentum 
tensor decomposition relative to the normal vectors (a perfect fluid will appear as an imperfect fluid in that 
frame). Various papers relate these equations to variational principles and HamUtonian formalisms, thus 
expressing them in terms of a potential formalism that gives an intuitive feel for what the evolution will be 
like ll48l I175j . There have also been many numerical investigations of these dynamical equations and the 
resulting solutions. We will briefly consider three specific aspects here, then the relation to observations, and 
finally the related dynamical systems approach. 

6.2.1 Chaos in these universes? 

An ongoing issue since Misner's discovery of the 'Mixmaster' behaviour of the Type IX universes has been 
whether or not these solutions show chaotic behaviour as they approach the initial singularity (see II 1761 and 

-^That is, they can be brought to the canonical forms of the C" by a suitable change of group-invariant basis (the final normalisation 
to ±1 may require changing from normalised basis vectors); the transformation to do so is different at each point and at each time. 
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Figure 2: Canonical structure constants for different Bianchi types. The Class B parameter h is defined as 
h = a? /n2n3 (see, e.g., lISTI ). 



Isotropic Bianchi models 
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FLRW k = 0: Bianchi I, Bianchi VII_0 
FLRW k = -1: Bianchi V, Bianchi VII_h 



LRS Bianchi models 



Orthogonal 







c \neq 



Taub-NUT I 
Taub-NUT 3 
Taub-NUT 2 



Tilted 



[KS +1: no subgroup] 
Bianchi I, VII_0 



Bianchi III 



"KS 



Bianchi IX 
Bianchi II 
Bianchi VII_h, 

III 



Bianchi V, VII_ 

Farnsworth, 

Collins-Ellis 



Figure 3: The Bianchi models permitting higher symmetry subcases. The parameter c is zero iff the preferred 
spatial vector is hypersurf ace-orthogonal. 
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Hobill in lISTl '). The potential approach represents these solutions as bouncing in an expanding approximately 
triangular shaped potential well, with three deep troughs attached to the corners. The return map approxima- 
tion (a series of Kasner-like epochs, separated by collisions with the potential walls and consequent change 
of the Kasner parameters) suggests the motion is chaotic, but the question is whether this map represents 
the solutions of the differential equations well enough to reach this conclusion (for example, the potential 
walls are represented as flat in this approximation; and there are times when the solution moves up one of 
the troughs and then reflects back, but the return map does not represent this part of the motion). Part of the 
problem is that the usual definitions of chaos in terms of a Lyapunov parameter depend on the definition of 
time variable used, and there is a good case for changing to conformal Misner time in these investigations. 

The issue may have been solved now by an analysis of the motion in terms of the attractors in phase 
space given by Cornish and Levin 1 1771 , suggesting that the motion is indeed chaotic, independent of the 
definition of time used. There may also be chaos in Type XIII solutions. Moreover, chaotic behaviour near 
the initial singularity was observed in solutions when a source-free magnetic Maxwell field is coupled to fluid 
space-times of Type I II178I and Type VIq II 1791 . 

6.2.2 Horizons and whimper singularities 

In tilted Class B models, it is possible for there to be a dramatic change in the nature of the solution. This 
occurs where the surfaces of homogeneity change from being spacelike (at late times) to being timelike (at 
early times), these regions being separated by a null surface H, the horizon associated with this change of 
symmetry. At earlier times the solution is no longer spatially homogeneous — it is inhomogeneous and sta- 
tionary|3 Associated with the horizon is a singularity where all scalar quantities are finite but components 
of the matter energy-momentum tensor diverge when measured in a parallelly propagated frame as one ap- 
proaches the boundary of space-time (this happens because the parallelly propagated frame gets infinitely 
rescaled in a finite proper time relative to a family of KV's, which in the limit have this singularity as a fixed 
point). The matter itself originates at an anisotropic big bang singularity at the origin of the universe in the 
stationary inhomogeneous region. 

Details of how this happens are given in Ellis and King 111801 . and phase plane diagrams for the simplest 
models in which this occurs — tilted LRS Type V models — in Collins and Ellis |170|. These models 
isotropise at late times, and can be arbitraiily similar to a low density FLRW model at the present day. 

6.2.3 Isotropisation properties 

An issue of importance is whether these models tend to isotropy at early or late times. An important paper by 
Collins and Hawking [172 J shows that for ordinary matter many Bianchi models become anisotropic at very 
late times, even if they are very nearly isotropic at present. Thus, isotropy is unstable in this case. On the 
other hand, a paper by Wald [181] showed that Bianchi models will tend to isotropise at late times if there is 
a positive cosmological constant present, implying that an inflationary era can cause anisotropics to die away. 
The latter work, however, while applicable to models with non-zero tilt angle, did not show this angle dies 
away, and indeed it does not do so in general (Goliath and Ellis [182]). Inflation also only occurs in Bianchi 
models if there is not too much anisotropy to begin with (Rothman and Ellis [,183] ), and it is not clear that 
shear and spatial curvature are in fact removed in all cases II 1 841 . Hence, some Bianchi models isotropise due 
to inflation, but not all. 

To study these kinds of question properly needs the use of phase planes. These will be discussed after 
briefly considering observations. 

^^This kind of change happens also in the maximally extended Schwarzschild solution at the event horizon. 
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6.3 Observational relations 

Observational relations in these models have been examined in detail. 

(a) Redshift, observer area distance, and galaxy observations ((M, z) and {N, z) relations) are considered 
in MacCallum and Ellis II128I . Anisotropies can occur in all these relations, but many of the models will 
display discrete isotropics in the sky. 

(b) The effect of tilt is to make the universe look inhomogeneous, even though it is spatially homogeneous 
(King and Ellis |116|). This will be reflected in particular in a dipole anisotropy in number counts, which 
will thus occur in rotating universes |2O|0 

(c) Element formation will be altered primarily through possible changes in the expansion time scale at 
the time of nucleosynthesis ( ||129| 11861 11871 ). This enables us to put limits on anisotropy from measured 
element abundances in particular Bianchi types. This effect could in principle go either way, so a useful 
conjecture [ 188 1 is that in fact the effect of anisotropy will always — despite the possible presence of rotation 
— be to speed up the expansion time scale in Bianchi models. 

(d) CBR anisotropies will result in anisotropic Universe models, e.g., many Class B Bianchi models will 
show a hot-spot and associated spiral pattern in the CBR sky II189I - II191|1^ This enables us to put limits on 
anisotropy from observed CBR anisotropy limits (Collins and Hawking II189I . Bunn et al II192I ). If 'reheat- 
ing' takes place in an anisotropic universe, this will mix anisotropic temperatures from different directions, 
and hence distort the CBR spectrum 01311 . 

Limits on present-day anisotropy from the CBR and element abundance measurements are very stringent: 
Icol/Bo < 10^^ to 10~^^, depending on the model. However, because of the anisotropies that can build up 
in both directions in time, this does not imply that the very early Universe (before nucleosynthesis) or late 
Universe will also be isotropic. The conclusion applies back to last scattering (CBR measurements) and to 
nucleosynthesis (element abundances). In both cases the conclusion is quite model dependent. Although very 
strong limits apply to some Bianchi models, they are much weaker for other types. Hence, one should be a 
bit cautious in what one claims in this regard. 

6.4 Dynamical systems approach 

The most illuminating description of the evolution of families of Bianchi models is a dynamical systems 
approach based on the use of orthonormal tetrads, presented in detail in Wainwright and Ellis |51|. The 
main variables used are essentially the commutation functions mentioned above, but rescaled by a common 
time dependent factor 

6.4.1 Reduced differential equations 

The basic idea (Collins [193], Wainwright ri941) is to write the Einstein's field equations in a way that enables 
one to study the evolution of the various physical and geometrical quantities relative to the overall rate of 
expansion of the model, as described by the rate of expansion scalar Q, or, equivalently, the Hubble scalar 
iJ = i 0. The remaining freedom in the choice of orthonormal tetrad needs to be eliminated by specifying 
the variables fi" implicitly or explicitly (e.g., by specifying them as functions of the aap)- This also simplifies 
the other quantities (e.g., choice of a shear eigenframe will result in the tensor Cap being represented by two 

-■*They will also occur in FLRW modeLs seen from a reference frame that is not comoving; hence, they should occur in the real 
universe if the standard interpretation of the CBR anisotropy as due to our motion relative to a FLRW universe is correct; see Ellis and 
Baldwin 1 1851. 

-^This result is derived in a gauge-dependent way; it would be useful to have a gauge-invariant version. 
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diagonal terms). One so obtains a reduced set of variables, consisting of H and the remaining commutation 
functions, which we denote symbolically by 

X = (T^hclreducod) • (173) 

The physical state of the model is thus described by the vector {H, x). The details of this reduction differ 
for the Class A and B models, and in the latter case there is an algebraic constraint of the form 

gix) - , (174) 

where g is a homogeneous polynomial. 

The idea is now to normalise x with the Hubble scalar H. We denote the resulting variables by a vector 
y £ M", and write: 

y = |: • (175) 

These new variables are dimensionless, and will be referred to as Hubble-normalised variables. It is 

clear that each dynamical state y determines a 1-parameter family of physical states {H, x). The evolution 
equations for the 7°^^ lead to evolution equations for H and x and hence for y. In deriving the evolution 
equations for y from those for x, the deceleration parameter q plays an important role. The Hubble scalar 
H can be used to define a scale factor S according to dZSl ) 

H=^, (176) 

where • denotes differentiation with respect to t. The deceleration parameter, defined hy q ~ — S S/S"^ (see 
(|58]l). is related to H according to 

H=-{l + q)H^. (177) 

In order that the evolution equations define a flow, it is necessary, in conjunction with the rescaling (11751) . to 
introduce a dimensionless time variable t according to 

S = Soe^, (178) 

where So is the value of the scale factor at some arbitrary reference time. Since S assumes values < S < 
+ 00 in an ever-expanding model, r assumes all real values, with r — ^ — 00 at the initial singularity and 
T — > + 00 at late times. It follows from equations ( |176l l and ( |178t that 

and the evolution equation (I177l t for H can be written 

= -{l + q)H. (ISO) 

ctr 

Since the right-hand sides of the evolution equations for the 7''f,c are homogeneous of degree 2 in the j°'bc, 
the change ( |179t of the time variable results in H canceling out of the evolution equation for y, yielding an 
autonomous differential equation ('DE'): 

^ = /(?/), yeM"- (181) 

The constraint g{x) = translates into a constraint 

giy) = , (182) 

which is preserved by the DE. The functions / : R" ^ M" and g : M" R are polynomial functions in y. 
An essential feature of this process is that the evolution equation for H, namely ( II8OI 1. decouples from the 
remaining equations ( II8II 1 and (1182b . In other words, the DE ( 1181b describes the evolution of the non-tilted 
Bianchi cosmologies, the transformation (1175b essentially scaling away the effects of the overall expansion. 
An important consequence is that the new variables are bounded near the initial singularity. 
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6.4.2 Equations and orbits 

The first step in the analysis is to formulate Einstein's field equations, using Hubble-normalised variables, as 
a DE (1181b in M", possibly subject to a constraint ( 1182b . Thus one uses the tetrad equations presented above, 
now adapted to apply to the variables rescaled in this way. Since r assumes all real values (for models which 
expand indefinitely), the solutions of ( 1181b are defined for all r and hence define a flow {(pr} on R". The 
evolution of the cosmological models can thus be analyzed by studying the orbits of this flow in the physical 
region of the state space, which is a subset of R" defined by the requirement that the matter energy density fj, 
be non-negative, i.e., 

niy) = J^>0, (183) 
where the density parameter fl (see ( l58b ) is a dimensionless measure of /i. 

The vacuum boundary, defined hy fl{y) = 0, describes the evolution of vacuum Bianchi models, and 
is an invariant set which plays an important role in the qualitative analysis because vacuum models can be 
asymptotic states for perfect fluid models near the big-bang or at late times. There are other invariant sets 
which are also specified by simple restrictions on y which play a special role: the subsets representing each 
Bianchi type (Figure 2), and the subsets representing higher symmetry models, specifically the FLRW models 
and the LRS Bianchi models (according to Figure 3). 

It is desirable that the dimensionless state space D in is a compact set. In this case each orbit 
will have non-empty future and past limit sets, and hence there will exist a past attractor and a future 
attractor in state space. When using Hubble-normalised variables, compactness of the state space has a 
direct physical meaning for ever-expanding models: if the state space is compact, then at the big-bang no 
physical or geometrical quantity diverges more rapidly than the appropriate power of H, and at late times 
no such quantity tends to zero less rapidly than the appropriate power of H. This will happen for many 
models; however, the state space for Bianchi Type VIIq and Type VIII models is non-compact. This lack of 
compactness manifests itself in the behaviour of the Weyl curvature at late times. 

6.4.3 Equilibrium points and self-similar cosmologies 

Each ordinary orbit in the dimensionless state space corresponds to a one-parameter family of physical mod- 
els, which are conformally related by a constant rescaling of the metric. On the other hand, for an equilib- 
rium point y* of the DE ( 1181b (which satisfies f{y*) — 0), the deceleration parameter g is a constant, i.e., 
q{y*) = q*, and we find 

In this case, however, the parameter Hq is no longer essential, since it can be set to unity by a translation of 
T,T ^ T + const; then (|179b imphes that 

Ht=-^, (184) 

so that by ( 11731 ) and ( |175b the commutation functions are of the form (const) x t^^. It follows that the 
resulting cosmological model is self-similar. It then turns out that to each equilibrium point of the DE f |7iS7| ) 
there corresponds a unique self-similar cosmological model. In such a model the physical states at different 
times differ only by an overall change in the length scale. Such models are expanding, but in such a way 
that their dimensionless state does not change. They include the spatially flat FLRW model {Vl = 1) and 
the Milne model {Vl = 0). All vacuum and orthogonal perfect fluid self-similar Bianchi solutions have been 
given by Hsu and Wainwright 1,195,1 . The equilibrium points determine the asymptotic behaviour of other 
more general models. 
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6.4.4 Phase planes 

Many phase planes can be constructed explicitly. The reader is referred to Wainright and Ellis (5T\ for a 
comprehensive presentation and survey of results attained so far. Several interesting points emerge: 

* Relation to lower dimensional spaces: it seems that the lower dimensional spaces, delineating higher 
symmetry models, can be skeletons guiding the development of the higher dimensional spaces (the more 
generic models). This is one reason why study of the exact higher symmetry models is of significance. The 
way this occurs is the subject of ongoing investigation (the key issue being how the finite dimensional dynam- 
ical systems corresponding to models with symmetry are imbedded in and relate to the infinite dimensional 
dynamical system describing the evolution of models without symmetry). 

* Identification of models in state space: the analysis of the phase planes for Bianchi models shows that 
the procedure sometimes adopted of identifying all points in state space corresponding to the same model is 
not a good idea. For example the Kasner ring that serves as a framework for evolution of many other Bianchi 
models contains multiple realizations of the same Kasner model. To identify them as the same point in state 
space would make the evolution patterns very difficult to follow. It is better to keep them separate, but to 
learn to identify where multiple realizations of the same model occur (which is just the equivalence problem 
for cosmological models). 

* Isotropisation is a particular issue that can be studied by use of these planes Ill96lll82ll . It turns out that 
even in the classes of non-inflationary Bianchi models that contain FLRW models as special cases, not all 
models isotropise at some period of their evolution; and of those that do so, most become anisotropic again 
at late times. Only an inflationary equation of state will lead to such isotropisation for a fairly general class 
of models (but in the tilted case it is not clear that the tilt angle will die away 11821 ): once it has turned off, 
anisotropic modes will again occur. 

An important idea that arises out of this study is that of intermediate isotropisation: namely, models 
that become very like a FLRW model for a period of their evolution but start off and end up quite unlike these 
models. It turns out that many Bianchi types allow intermediate isotropisation, because the FLRW models 
are saddle points in the relevant phase planes. This leads to the following two interesting results: 

Bianchi Evolution Theorem (1): Consider a family of Bianchi models that allow intermedi- 
ate isotropisation. Define an e-neighbourhood of a FLRW model as a region in state space where 
all geometrical and physical quantities are closer than e to their values in a FLRW model. Choose 
a time scale L. Then no matter how small e and how large L, there is an open set of Bianchi 
models in the state space such that each model spends longer than L within the corresponding 
e-neighbourhood of the FLRW model. 

(This follows because the saddle point is a fixed point of the phase flow; consequently, the phase flow 
vector becomes arbitrarily close to zero at all points in a small enough open region around the FLRW point 
in state space.) 

Consequently, although these models are quite unlike FLRW models at very early and very late times, 
there is an open set of them that are observationally indistinguishable from a FLRW model (choose L long 
enough to encompass from today to last coupling or nucleosynthesis, and e to correspond to current observa- 
tional bounds). Thus, there exist many such models that are viable as models of the real Universe in terms of 
compatibility with astronomical observations. 

Bianchi Evolution Theorem (2): In each set of Bianchi models of a type admitting interme- 
diate isotropisation, there will be spatially homogeneous models that are linearisations of these 
Bianchi models about FLRW models. These perturbation modes will occur in any almost-FLRW 
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model that is generic rather than fine-tuned; however, the exact models approximated by these 
linearisations will be quite unlike FLRW models at very early and very late times. 

(Proof is by linearising the equations above (see the following section) to obtain the Bianchi equations 
linearised about the FLRW models that occur at the saddle point leading to the intermediate isotropisation. 
These modes will be the solutions in a small neighbourhood about the saddle point permitted by the linearised 
equations (given existence of solutions to the non-linear equations, linearisation will not prevent correspond- 
ing linearised solutions existing).) 

The point is that these modes can exist as linearisations of the FLRW model; if they do not occur, then 
initial data has been chosen to set these modes precisely to zero (rather than being made very small), which 
requires very special initial conditions. Thus, these modes will occur in almost all almost-FLRW models. 
Hence, if one believes in generality arguments, they will occur in the real Universe. When they occur, they 
will at early and late times grow until the model is very far from a FLRW geometry (while being arbitrarily 
close to an FLRW model for a very long time, as per the previous theorem). 

Exercise: Most studies of CBR anisotropics and nucleosynthesis are carried out for the Bianchi types 
that allow FLRW models as special cases (see Figure 3). Show that Bianchi models can approximate FLRW 
models for extended periods even if they do not belong to those types. What kinds of CBR anisotropies can 
occur in these models? (See, e.g., [51 J .) 



7 Almost-FLRW models 

The real Universe is not FLRW because of all the structure it contains, and (because of the non-linearity 
of Einstein's field equations) the other exact solutions we can attain have higher symmetry than the real 
Universe. Thus, in order to obtain realistic models we can compare with detailed observations, we need to 
approximate, aiming to obtain 'almost-FLRW' models representing a universe that is FLRW- like on a large 
scale but allowing for generic inhomogeneities on a small scale. 



7.1 Gauge problem 

The major problem in studying perturbed models is the gauge problem, due to the fact that there is no 
identifiable fixed background model in General Relativity. One can start with a unique FLRW Universe 
model with metric in some local coordinate system, and perturb it to obtain a more realistic model: 
5ab ~^ 9ab = Vab + <^9ab, but then the process has no unique inverse: the background model is not 
uniquely determined by the lumpy Universe model gab (no unique tensorial averaging process has been 
defined that will recover '§^1, from gab)- Many choices can be made. However, the usual variables describing 
perturbations depend on the way the (fictitious) background model g^f, is fitted to the metric of the real 
Universe, gab', these variables can be given any values one wants by changing this correspondence For 
example, the dimensionless energy density contrast S representing a density perturbation is usually defined 

Si.) = , (185) 

where is the actual value of the energy density at the point x, while //(x) is the (fictitious) background 
value there, determined by the chosen mapping of the background model into the realistic lumpy model (Ellis 
and Bruni L197J). This quantity can be given any value we desire by altering that map; we can, e.g., set 
it to zero by choosing the real surfaces of constant energy density (provided these are spacelike) to be the 



-*This is often represented implicitly rather than explicitly, by assuming that points with the same coordinate values in the background 
space and more realistic model map to each other; then the gauge freedom is contained in the coordinate freedom available in the realistic 
universe model. 



7 ALMOST-FLRW MODELS 



53 



background surfaces of constant time and, hence, of constant energy density. Consequently, perturbation 
equations written in terms of this variable have as solution both physical modes and gauge modes, the latter 
corresponding to variation of gauge choice rather than to physical variation. 

One way to solve this is by very carefully keeping track of the gauge choice used and the resulting gauge 
freedom; see Ma and Bertschinger |198| and Prof. Bertschinger's lectures here. The alternative is to use 
gauge-invariant variables. A widely used and fundamentally important set of such variables are those intro- 
duced by Bardeen 1 199], and used for example by Bardeen, Steinhardt and Turner [200|. Another possibility 
is use of gauge-invariant and 1 + 3 covariant ('GIC') variables, i.e., variables that are gauge-invariant and 
also 1 + 3 covariantly defined so that they have a clear geometrical meaning, and can be examined in any 
desired (spatial) coordinate system. That is what will be pursued here. 

In more detail: our aim is to examine perturbed models by using 1 + 3 covariant variables defined in the 
real space-time (not the background), deriving exact equations for these variables in that space-time, and then 
approximating by linearising about a RW geometry to get the linearised equations describing the evolution 
of energy density inhomogeneities in almost-FLRW universes. How do we handle gauge invariance in this 
approach? We rely on the 

Gauge Invariance Lemma (Stewart and Walker [2011): If a quantity T ... vanishes in the 
background space-time, then it is gauge-invariant (to all orders). 

[The proof is straightforward : If T' ... = 0, then ST -\,, = T- ... - T' ... = T which is manifestly 
independent of the mapping <!> from 5' to S" (it does not matter how we map T ... from S* to S* when T ... 
vanishes).] The application to almost-FLRW models follows (see Ellis and Bruni |197|), where we use an 
order-of-magnitude notation as follows: Given a smallness parameter e, 0[n] denotes ©(e"), and A ^ B 
means A — B = 0[2] (i.e., these variables are equivalent to 0[1]). When A w we shall regard A as 
vanishing when we linearise (for it is zero to the accuracy of relevant first-order calculations). Then, 

• Zeroth-order variables are p, p, 8, and their covariant time derivatives, fi, p, 6, 

• First-order variables are u'^, aab, <f\ T^ab, Eab, Hab, Xa, Za, and their covariant time and space 
derivatives. 

As these first-order variables all vanish in exact FLRW models, provided m° is uniquely defined in the real- 
istic (lumpy) almost-FLRW Universe model, they are all uniquely defined GIC variables. Thus, this set of 
variables provides what we wanted: 1 + 3 covariant variables characterising departures from a FLRW geom- 
etry (and, in particular, the spatial inhomogeneity of a universe) that are gauge-invariant when the universe 
is almost-FLRW. Because they are tensors defined in the real space-time, we can evaluate them in any local 
coordinate system we like in that space-time. 

7.1.1 Key variables 

Two simple gauge-invariant quantities give us the information we need to discuss the time evolution of energy 
density fluctuations. The basic quantities we start with are the orthogonal projections of the energy density 
gradient, i.e., the vector Xa = Va/i, and of the expansion gradient, i.e., the vector Za = VaS. The first 
can be determined (a) from virial theorem estimates and large-scale structure observations (as, e.g., in the 
POTENT programme), (b) by observing gradients in the numbers of observed sources and estimating the 
mass-to-light ratio ( Kristian and Sachs |98 1, Eq. (39) ), and (c) by gravitational lensing observations. How- 
ever, these do not directly correspond to the quantities usually calculated; but two closely related quantities 
do. The first is the matter-comoving fractional energy density gradient: 

Va = S—, (186) 
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which is gauge-invariant and dimensionless, and represents the spatial energy density variation over a fixed 
comoving scale. Note that S, and so Va, is defined only up to a constant by equation (l25l i: this allows it to 
represent the energy density variation between any neighbouring worldlines. The vector Va can be separated 
into a magnitude V and direction Ca 

Va=Vea, eae" = 1 , eau" = ^ {VaV^y/^ . (187) 

The magnitude V is the gauge-invariant variablj^ that most closely corresponds to the intention of the usual 
6 = Sfj,/]! given in ( 1185b . The crucial difference from the usual definition is that T) represents a (real) spatial 
fluctuation, rather than a (fictitious) time fluctuation, and does so in a GIC manner. An important auxiliary 
variable in what follows is the matter-comoving spatial expansion gradient: 

Za = SZa. (188) 

The issue now is, can we find a set of equations determining how these variables evolve? Yes we can; 
they follow from the exact 1 + 3 covariant equations of subsection 12. 2 1 

7.2 Dynamical equations 

We can determine exact propagation equations along the fluid flow lines for the quantities defined in the 
previous section, and then linearise these to the almost-FLRW case. The basic linearised equations are given 
by Hawking [22j (see his equations (13) to (19)); we add to them the linearised propagation equations for the 
gauge-invariant spatial gradients defined above M197I . 

7.2.1 Growth of inhomogeneity 

Taking the spatial gradient of the equation of energy conservation ( |37] | (for the case of a perfect fluid), we 

find ma 

Va(/i) + e Va(M + p) + (m + p) v^e = o , 

i.e., 

Using Leibniz' Rule and changing the order of differentiation in the second-derivative term (and noting 
that the pressure-gradient term cancels on using the momentum conservation equation (|38]|), we obtain the 
fundamental equation for the growth of inhomogeneity: 

+ f e X" = - (/i + p) - a\ + jf^" Xc , (189) 

with source term Z° . On taking the spatial gradient of the Raychaudhuri equation (l29b . we find the companion 
equation for that source term: 

^(a> + e z° = - i X" - (T% z'' + T]"'"' Zc + u^n-i v"(cr2 - i^^) + v''(VbM^ + ubv!') , (i90) 

where 

7^ = - i + Vau" + {iiau") - 2ct2 + 2cj2 + ^ + A. (191) 

These exact equations contain no information not implied by the others already given; nevertheless, they 
are useful in that they are exact equations directly giving the rate of growth of inhomogeneity in the generic 
(perfect fluid) case, the second, together with the evolution equations above, giving the rate of change of all 
the source terms in the first. 



Or, equivalently in the linear case, tlie spatial divergence of T>a ■ 
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The procedure now is to systematically approximate all the dynamical equations, and, in particular, the 
structure growth equations given above, by dropping all terms of second order or higher in the implicij^ 
expansion variable e. Thus, we obtain the linearised equations as approximations to the exact equations 
above by noting that in the almost-FLRW case, 

X<"> + feX" = -{p + p)Z'' + 0[2] (192) 

Z^^^+eZ" = -iX" + ^^°7^ + V"(V6w'') + 0[2] , (193) 

where 

7e = -ie^ + VaW^ + Af + A + Op] . (194) 

Then we linearise the equations by dropping the terms 0[2], so from now on in this section '=' means equal 
up to terms of order e^. 

7.3 Dust 

In the case of dust, p = = 0, and the equations (|192| i and ( |193l l for growth of inhomogeneity become 

S-^h''b{S'^X''y = -pZ"" (195) 
S-^h''b{S^Z''y ^ -^X", (196) 

This closes up to give a second-order equation (take the covariant time derivative of the first and substitute 
from the second and the energy conservation equation (IJTTi). To compare with the usual equations, change to 
the variables Va and ( see ( |186t and ( |188t ). Then the equations become 

V^"^ = -Z" (197) 
i^''^ = -^QZ" - ^pV . (198) 

These directly imply the second-order equation (take the covariant time derivative of the first equation!) 

= + I ei>(,)-iM2?a, (199) 

which is the usual equation for growth of energy density inhomogeneity in dust universes, and has the usual 
solutions: when fc = 0, therl^ S{t) cx t^/'^, and we obtain 

= d+,(x")i2/3 + d_ja;")i-i , d,a=0, (200) 

(where t is proper time along the flow lines)I3 This shows the growing mode that leads to structure formation 
and the decaying mode that dissipates previously existing inhomogeneities. It has been obtained in a GIC 
way: all the first-order variables, including in particular those in this equation, are gauge-invariant, and there 
are no gauge modes. Furthermore, we have available the fully non-linear equations, Eqs. ( I189l l and ( 11901) . 
and so can estimate the errors in the neglected terms, and set up a systematic higher-order approximation 
scheme for solutions of these equations. Solutions for other background models (with fc = ± 1 or A 7^ 0) can 
be obtained by substituting the appropriate values in ( |199l l for the background variables Q and p, possibly 
changing to conformal time to simplify the calculation. 

Exercise: What is the growth rate at late times in a low-density universe, when the expansion is curvature 
dominated and so is linear: S{t) — cq tl What if there is a cosmological constant, so the late time expansion 
is exponential: S{t) = cq exp{Ho t), where cq and Hq are constants? 

-''One could make this variable explicit, but there does not seem to be much gain in doing so. 

-'This is the background rather than the real value of this quantity, which is what should really be used here; it is determined in the 
real space-time by 5/5" = ^ ©■ However, as we are linearising, the difference makes only a second-order change in the coefficients in 
the equations, which we can neglect. 

^"We can also see that if = 0, there will be an exponential rather than power-law growth. 
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7.3.1 Other quantities 

We have concentrated here on the growth of inhomogeneities in the energy density and the expansion rate. 
However, all the 1 + 3 covariant equations of the previous sections apply and can be linearised in a straight- 
forward way to the almost-FLRW case (and one can find suitable coordinates and tetrads in order to employ 
the complete set of tetrad equations in this context, too). Doing so, one can in particular study vorticity per- 
turbations and gravitational wave perturbations; see the pioneering paper by Hawking |22|. We will not 
consider these further here; however, a series of interesting issues arise. The GIC approach to gravitational 
waves is examined in Hogan and Ellis 02O2I and Dunsby, Bassett and Ellis 1231 . 

A further important issue is the effect of perturbations on observations (apart from the CBR anisotropies, 
discussed below). It has been known for a long time that anisotropies can affect area distances as well as 
redshifts (see Bertotti |203|, Kantowski 1 159|); the Dyer-Roeder formula (P2041, see also |99|) can be used 
at any redshift for those many rays that propagate in the lower-density regions between inhomogeneities; 
however, this formula is not accurate for those ray bundles that pass very close to matter, where shearing 
becomes important. This is closely related to the averaging problem (see, e.g., Ellis |88 1 or Boersma 1205]): 
how can dynamics and observations of a Universe model which is basically empty almost everywhere average 
out correctly to give the same dynamics and observations as a Universe model which is exactly spatially 
homogeneous? What differences will there be from the FLRW case? We will not pursue this further here 
except to state that it is believed this does in fact work out OK: in the fully inhomogeneous case it is the 
Weyl curvature that causes distortions in the empty spaces between astrophysical objects (as in gravitational 
lensing), and, hence, causes convergence of both timelike and null geodesic congruences; these, however, 
average out to give zero average distortion and the same convergence effect as a FLRW space-time with zero 
Weyl curvature, but with the Ricci curvature causing focusing of these curves (see |206| for a discussion of 
the null case; however, some subtleties arise here in terms of the way areas are defined when strong lensing 
takes place M207lll45ll208l ). 

7.4 Perfect fluids 

A GIC analysis similar to the dust case has been given by Ellis, Hwang and Bruni determining FLRW pertur- 
bations for the perfect fluid case II209I [191 . This gives the single-fluid equation for growth of structure in an 
almost-FLRW Universe model (again, derived in a GIC manner), and includes as special cases a fully 1 + 3 
covariant derivation of the Jeans length and of the speed of sound for barotropic perfect fluids. To evaluate 
the last two terms in ( fT93] l when p ^ 0, using dSSl l we see that, to first order. 



Va(VbV'': 



'p) 



(201) 



But, for simplicity considering only the case of vanishing vorticity, to' 



00 we have 



(202) 



and, on using the Ricci identities for the V-derivatives and the zeroth-order relation "^Rab 
3-dimensional Ricci tensor, we obtain 



i hab for the 



Thus, on using ^R — JC — 6k /S^, we find 



(203) 




(204) 



When it is non-zero, (27) must be taken into account when commuting derivatives; see |19| . 
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introducing the notation V^V aP = b^^'^ aP- In performing this calculation, note that there will not be 
3-spaces orthogonal to the fluid flow if lo'^ ^ 0, but still we can calculate the 3-dimensional orthogonal 
derivatives as usual (by using the projection tensor /i^h); the difference from when = will be that the 
quantity we calculate as a curvature tensor, using the usual definition from commutation of second derivatives, 
will not have all the usual curvature tensor symmetries. Nevertheless, the zeroth-order equations, represent- 
ing the curvature of the 3-spaces orthogonal to the fluid flow in the background model, will agree with the 
Unearised equations up to the required accuracy. 

Now \fp — p{^, s), where s is the entropy per particle, we find 

\ / s=const V '^•^ / /j=const 

We assume we can ignore the second term (pressure variations caused by spatial entropy variations) relative to 
the first (pressure variations caused by energy density variations) and spatial variations in the scale function S 
(which would at most cause second-order variations in the propagation equations). Then (ignoring terms due 
to the spatial variation of dp/dfi, which will again cause second-order variations) we find in the zero-vorticity 
case, 

^[i/C.. + V.(V..^)] = -^^ (I) + ■ (206) 

This is the result that we need in proceeding with (|193l l. 



7.4.1 Second-order equations 

The equations for propagation can now be used to obtain second-order equations for PqEI For easy compar- 
ison, we follow Bardeen 1 199 1 by defining 

w=-, cl = ^ ^ (-\=w = -Q{l + w){cl-w). (207) 

Now covariant differentiation of ( |192| i. projection orthogonal to and linearisation gives a second-order 
equation for Va (we use Eqs. ( |29] ), ( |193l l, ( 12071 ) and ( 12061 ) in the process). We find 

-[2k 

= V^,^ + {j-2w + cl)eV^,)~[{^+iw~lw^~3cl)fi+{cl-w) — ]Va (208) 

This equation is the basic result of this subsection; the rest of the discussion examines its properties and 
special cases. It is a second-order equation determining the evolution of the GIC energy density variation 
variable Va along the fluid flow lines. It has the form of a wave equation with extra terms due to the expan- 
sion, gravitation and the spatial curvature of the Universe model|3 We bracket the last two terms together, 
because when we make a harmonic decomposition these terms together give the harmonic eigenvalues n^. 

This form of the equations allows for a variation of w — p/fi with time. However, if w = const, then 
from ( |207b — w, and the equation simplifies to 

2k 

= V^a) + {l~w)eV^a) ~^{l~w){l + 3w)^lVa + w[—Va'V'Va]. (209) 

The matter source term vanishes if w — 1 (the case of 'stiff matter' p — fi) or w = — ^ (the case 
p = — i /i, corresponding to matter with no active gravitational mass). Between these two limits ('ordinary 

^^And the variable <l?a = fiS^ "Da that corresponds more closely to Bardeen's variable; see 1191 . 
'^We have dropped A in these equations; it can be represented by setting w = — 1. 
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matter'), the matter term is positive and tends to cause the energy density gradient to increase ('gravitational 
aggregation'); outside these Hmits, the term is negative and tends to cause the energy density gradient to 
decrease ('gravitational smoothing'). Also, the sign of the damping term (giving the adiabatic decay of inho- 
mogeneities) is positive if | > w (that is, 2^ > 3p) but negative otherwise (they adiabatically grow rather 
than decay in this case). The equation reduces correctly to the corresponding dust equation in the case w = 0. 
Two other cases of importance are: 

* Speed of sound: when 8, fi, and k/S'^ can be neglected, we see directly from ( 12081 ) that Cs introduced 
above is the speed of sound (and that imaginary values of Cs, i.e., negative values of dp/dfx, then lead to 
exponential growth or decay rather than oscillations). 



* Radiation: In the case of pure radiation, 7 = | and w = i = c^. Then we find from ( |209t 

\ 1 r 

^) + 3 [ ^ 



3fc 2k 

= +iei?(a) + -^)Va + \[^Va~V^Va]. (210) 



7.4.2 Harmonic decomposition 



It is standard (see, e.g., 11221 and 01991 ) to decompose the variables harmonically, thus effectively separating 
out the time and space variations by turning the differential equations for time variation of the perturbations 
as a whole into separate time variation equations for each component of spatial variation, characterised by a 
matter-comoving wavenumber. This conveniently represents the idea of a matter-comoving wavelength for 
the matter inhomogeneities. In our case we do so by writing Va in terms of harmonic vectors q\^\ from 
which the background expansion has been factored out. 

We start with the defining equations of the 1 + 3 covariant scalar harmonics (3*^"^ 

2 

Q(")=0, V2g(") =-^Q(") , (211) 

corresponding to Bardeen's scalar Helmholtz equation (2.7) II199L but expressed 1 + 3 covariantly following 
Hawking ll22l . From these quantities we define the 1 + 3 covariant vector harmonics (cf. M199L equations 
(2.8) and (2.10); we do not divide by the wavenumber, however, so our equations are valid even if n = 0) 

Q(")^^V,Q(") ^ Oi-^^'^-O, qS«0, v^q(») = _ , (212) 

(the factor S ensuring these vector harmonics are approximately covariantly constant along the fluid flow 
lines in the almost-FLRW case). We write Va in terms of these harmonics: 

Va = Qi") , Val?(") « , (213) 

n 

where I?^"^ is the harmonic component of Va corresponding to the matter-comoving wavenumber n, con- 
taining the time variation of that component; to first order, I?'"' = (5/i'"-'//I. Putting this decomposition in 
the linearised equations ( 1208b and ( 12091 ). the harmonics decouple. Thus, e.g., we obtain from ( 12091 ) the rt-th 
harmonic equation 

„2 

= 2j(") + (|-w)ep(")-[i(l-u;)(l + 3w)Ai-w— , (214) 

(valid for each n > 0), where one can, if one wishes, substitute for /i from the zeroth-order Friedmann 
equation in terms of 8 and k/S'^. This equation shows how the growth of the inhomogeneity depends on the 
matter-comoving wavelength. For the case of radiation, w — ^, this is 

= 2?(") + iei>(")-[|Ai-3-^]2?^"^ . (215) 
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7.5 Implications 

To determine the solutions explicitly, we have to substitute for pi, 9 and S from the zeroth-order equations. 
The most important issue is which terms dominate. 

7.5.1 Jeans instability 

Jeans' criterion is that gravitational collapse will tend to occur if the combination of the matter term and the 
term containing the Laplace operator in ( 1208b or ( |209t is positive II210I : i.e., if 

2k 

^ {1 ~ W) {I + 3w) pVa > W[ ^Va - V^Va ] , (216) 
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when = w. Using the harmonic decomposition, this can be expressed in terms of an equivalent scale: from 
(1214b gravitational coUapse tends to occur for a mode 2?'"^ if 



n2 



l{\-w){l + iw)p>w—, (217) 



i.e., if 

1/2 

S{t) > n . (218) 



In terms of wavelengths the Jeans length is defined by 



nj -'''''^ Gpit) {l~w){l + 3w) ' ^^'^^ 

where we have re-established the fundamental constants c and G (w — (cs/c)^). Thus, gravitational col- 
lapse will occur for small n (wavelengths longer than Aj), but not for large n (wavelengths less than Aj), 
for the pressure gradients are then large enough to resist the collapse and lead to acoustic oscillations instead. 

For non-relativistic matter, \w\ ^ 1 and /i = Pm(?, where pm is the mass density of the matter, so 



Then the Jeans mass will be 



47r 47r „ /7r\3/2 , 

Mj = — p„ A3 = — cl (- j (221) 



For radiation, where w = ^ and /i = PrC^, collapse will occur if 



37r 



iGprit) 

The corresponding Jeans mass for matter coupled to the radiation will be 



(222) 



M,.-p,„A3, = -c^(^^j P™. (223) 
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7.5.2 Short-wavelength solutions 

For wavelengths much shorter than the Jeans length, equation (I214l l becomes the damped harmonic equation 

= 2j(») + (2_^)eX)(")+u; (224) 

giving oscillations. In the early Universe, during radiation dominated expansion before decoupling, the tight 
coupling of the dominant radiation and matter leads to a fluid with w = i; then the short-wavelength equation 
becomes 

2 

= 2?(") + iei>(") + i^I?("\ (225) 
giving the acoustic oscillations during that era for modes such that \ < Xj — CsC {Stt/AG^Y^^. 

7.5.3 Long-wavelength solutions 

For wavelengths much longer than the Jeans length, we can drop the Laplace operator terms in (I2141 l to obtain 

= + (| - w) e -i(l-u;)(l + 3«;)Ai2?^"^ . (226) 

Thus, the second-order propagation equations become ordinary differential equations along the fluid flow 
lines, easily solved for particular equations of state. In the case of radiation (w — ^) we find 

= 2?(") + i e i>(") - 1 j^i + , (227) 

when A > Aj = c,c{3n/4G^iy/'^. When fc = 0, then S{t) cx t^^^, and we obtain in this long-wavelength 
limit 

Va=d+a{x'')t + d^a{x^)t-^^\ d^a ^ , (228) 

(where t is proper time along the flow lines). The corresponding standard results in the synchronous and 
matter-comoving proper time gauges differ, being modes proportional to t and to t^^^; we obtain the same 
growth law as derived in the matter-comoving time orthogonal gauge and equivalent gauges. As our variables 
are GIC, we believe they show the latter gauges represent the physics more accurately than any other Note 
that, moreover, we obtain no fictitious modes (proportional to t^^) because we are using GIC variables. 

7.5.4 Change of behaviour with time 

Any particular inhomogeneity will have a constant matter-comoving size and, hence, constant matter-comoving 
wavelength A and constant matter-comoving wavenumber n as defined above. However, the Jeans length will 
vary with time. 

During the radiation era S cx t^/"^ and ^ = | i ^ (see JllSl l), so (dropping the dimensional constants) 
Pm oc t^^/"^, and the matter-comoving Jeans length 



Aj = i7r<2 (X t (229) 



will steadily grow to a value X^^^ at matter-radiation equality, while the Jeans mass of coupled matter will 

grow as 

Mj = ^ p,„ X^j cx = ^3/2 . (230) 



■^This is early enough that we can ignore the curvature term in the Friedmann equation. 
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Thereafter, until recombination, the Jeans mass stays constant: the matter and radiation are still tightly cou- 
pled but now the Universe model is matter-dominated and the speed of sound of the coupled fluid depends 
on the matter density: Cs{c/^/3) ( l + (3p„i/4pr) ) (seeRees [211]). After recombination the Jeans length and 
mass will rapidly die away because the matter and radiation decouple, leading to Cs ^ and so Xj 0. Each 
wavelength A longer than Xf^^ will have a growing mode as in (I228l l until the Jeans wavelength becomes 
greater than A; it will then stop growing and undergo acoustic oscillations which will last until decoupling 
when the Jeans length drops towards zero and matter-dominated growth starts according to ( I200l i. Growth of 
small perturbations eventually slows down when the Universe model becomes curvature-dominated at late 
times (when this happens depends on Oq; in a critical density universe, it never occurs). 

Thus, the key times for any wavelength after the initial perturbations have been seedec^l are (i) tj, when 
they become smaller than the Jeans length (if they are small enough that this occurs), (ii) t^qu, when the 
matter-dominated era starts (which will be before decoupling of matter and radiation, because fto > 0.1), 
(iii) tdcc, when decoupling takes place. The acoustic oscillations have constant amplitude in the radiation- 
dominated era from tj until tcqu, and then die away as t^^/^ in the matter-dominated era until they end at tdcc 
1121 II . Baryonic inhomogeneities 'freeze out' at that time; they then start growing by damped gravitational 
attraction. If they grow large enough, changing to non-linear collapse and ultimately star formation, then 
local energy generation starts. 

Exercise: Establish these behaviours from the equations given above. 

By contrast CDM freezes out at tcqu and starts growth at that time (Rees ||97l ). Thus, in a CDM- 
dominated Universe model, as is often supposed, the CDM fluctuations that govern structure formation 
start gravitational growth earlier than the baryons. They then govern the growth of inhomogeneities, attract- 
ing the baryons into their potential wells; a 2-fluid description representing the separate average velocities 
and their relative motion (see below) is needed to examine this. 

This picture has to be modified, however, by allowing for diffusion effects. Kinetic theory is the best way 
to tackle this. The result is damping of perturbations below diffusion scales which depend on whether or not 
Hot Dark Matter ('HDM') is present; baryonic fluctuations on small scales are attenuated by photon viscosity 
and free-streaming of neutrinos M211l[97l . 

7.6 Other matter 

Many other cases have been examined in this GIC formalism. We list them with major references. 

7.6.1 Scalar fields 

The case of scalar fields is dealt with in a GIC way by Bruni, Ellis and Dunsby 12121 . This analysis leads to 
the usual conserved quantities and theory of growth of inhomogeneities in an inflationary era. A key element 
here is choice of 4-velocity; for small perturbations there is a unique obvious choice, namely, choosing 
orthogonal to the surfaces on which the scalar field </> is constant (assuming these are spacelike). The energy- 
momentum tensor then has the form of a perfect fluid, but with energy density and pressure depending on 
both kinetic and potential energy terms for 0. 

7.6.2 Multi-fiuids and imperfect fiuids 

The physically important case of multi-fluids is dealt with by Dunsby, Bruni and Ellis 12131 ; e.g., enabling 
modelling of perturbations that include a matter-radiation interaction. The key element again is choice of 

^'The usual assumption is tliat perturbations are essentially unaffected by all the strong interactions in the early universe after the end 
of inflation, including the ending of pair production (when matter ceases to be relativistic), decoupHng of neutrinos, and the irreversible 
interactions during baryosynthesis and nucleosynthesis, and decoupling with Hot Dark Matter ('HDM') or CDM. 
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4-velocity. Each component has a separate 4-velocity u^-.y and there are various options now for the ref- 
erence 4-velocity u". When Hnear changes are made in this choice of 4-velocity, the essential effect is to 
alter the measured momentum in the ii"-frame (Maartens et al II214I ). The equations are simplified most by 
choosing as the centre of mass 4-velocity for the sum of all components, and the 4-velocity of the 

1- th component as its centre of mass 4-velocity. One must then carefully check the separate momentum and 
energy equations for each component, as well as for the matter as a whole. These determine the evolution of 
the 1 + 3 covariantly defined relative velocities: — u'^-^ — u^^y and the separate matter densities 

Exercise: Establish the equations for the relative velocity and the energy density inhomogeneities in the 

2- fluid case. 

The case of imperfect fluids is closely related, and the same issue of choice of 4-velocity arises. As 
pointed out earlier, it is essential to use realistic equations of state in studying perturbations of an imperfect 
fluid, such as described by the Miiller-Israel-Stewart theory 1121 51 12161 (see Maartens and Triginier [217] for 
a detailed GIC analysis of such imperfect fluids). 

7.6.3 Magnetic fields 

These have been examined in a GIC way by Barrow and Tsagas |218|, using the 1 + 3 covariant spUtting of 
the electromagnetic field and Maxwell's field equations ITJ. 

7.6.4 Newtonian version 

A Newtonian version of the analysis can be developed fully in parallel to the relativistic version 02191 . based 
on the Newtonian analogue of the approach to cosmology presented in |6 1, and including derivations of the 
Newtonian Jeans length and Newtonian formulae for the growth of inhomogeneities. 

Exercise: Establish these equations, and hence determine the main differences between the Newtonian 
and relativistic versions of structure formation. 

7.6.5 Alternative gravity 

The same GIC approach can be used to analyze higher-derivative gravitational theories; details are in 022OI . 

7.7 Relation to other formalisms 

The relation between the GIC approach to perturbations and the very influential gauge-invariant formaUsm 
by Bardeen Ml 991 has been examined in depth 02211 . The essential points are that 

* as might be expected, the implications of both approaches for structure formation are the same, 

* the implications of the GIC formalism can be worked out in any desired local coordinate system, in- 
cluding Bardeen's coordinates (which are incorporated into that approach in an essential way from the start), 

* Bardeen's approach is essentially based on the linearised equations, while the GIC starts with the fuU 
non-linear equations and linearises them, as explained above. This enables an estimate of the errors involved, 
and a systematic n-th order approximation scheme, 

* the GIC formalism does not use a non-local splitting 02221 into scalar, vector and tensor modes, and 
only uses a harmonic splitting (into wavelengths) at a late stage of the analysis; these are both built into 
Bardeen's approach ab initio. 

Both approaches have the advantage over gauge approaches that they do not involve gauge modes, and the 
differential equations are of minimal order needed to characterise the physics of the problem. Many papers 
on the use of various gauges are rather confused; however, the major paper by Ma and Bertschinger [|1981 
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clarifies the relations between important gauges in a clear way. As shown there, the answers obtained for 
large-scale growth of inhomogeneities is indeed gauge-dependent, and this becomes significant particularly 
on very large scales. The GIC formalism obviates this problem. However, whatever formalism is used, the 
issue is how the results relate to observations. The perturbation theory predicts structure growth, gravitational 
wave emission, gravitational lensing, and background radiation anisotropics. The latter are one of the most 
important tests of the geometry and physics of perturbed models, and are the topic of the final section. 



8 CBR anisotropies 

Central to present day cosmology is the study of the information obtainable from measurements of CBR 
anisotropies. A GIC version of the pioneering Sachs-Wolfe paper II223I . based on photon path integration 
and calculation of the redshift along these paths (cf. the integration in terms of Bardeen's variables by Panek 
| 224| ). is given by Dunsby [2251 and Challinor and Lasenby II226I . However, a kinetic theory approach 
enables a more in-depth study of the photons' evolution and interactions with the matter inhomogeneities, 
and so is the dominant way of analyzing CBR anisotropies. 



8.1 Covariant relativistic kinetic theory 

Relativistic kinetic theory (see, e.g., Il68l and 0227I - II229I ) provides a self-consistent microscopically based 
treatment where there is a natural unifying framework in which to deal with a gas of particles in circum- 
stances ranging from hydrodynamical to free-streaming behaviour. The photon gas undergoes a transition 
from hydrodynamical tight coupling with matter, through the process of decoupling from matter, to non- 
hydrodynamical free-streaming. This transition is characterised by the evolution of the photon mean free 
path length from effectively zero to effectively infinity. This range of behaviour can appropriately be de- 
scribed by kinetic theory with non-relativistic classical Thomson scattering (see, e.g., Jackson [230] or 
Feynman 1 112311 ). and the baryonic matter with which radiation interacts can reasonably be described hydro- 
dynamically during these times. 

In this approach, the single-particle photon distribution function f{x'^,p°') over a 7-dimensional phase 
space 1681 represents the number of photons measured in the 3-volume element dV at the event that have 
momenta in the momentum space volume element vr about the momentum through the equation 

dN = f{x\p'')i~Pau'')7rdV , (231) 

where m° is the observer's 4-velocity and the redshift factor {—pau"") makes / into a (observer-independent) 
scalar. The rate of change of / in photon phase space is determined by the relativistic generalisation of 
Boltzmann's equation 

L{f) = C[f] , (232) 

where the Liouville operator 

Lif) ^dl^ + K^ (233) 
dx^ dv dp°- dv 

gives the change of / in parameter distance dv along the geodesies that characterise the particle motions. 
The collision term C[f] determines the rate of change of / due to emission, absorption and scattering pro- 
cesses; it can represent Thomson scattering, binary collisions, etc. Over the period of importance for CBR 
anisotropies, i.e., considerably after electron-positron annihilation, the average photon energy is much less 
than the electron rest mass and the electron thermal energy may be neglected so that the quantised Compton 
interaction between photons and electrons (the dominant interaction between radiation and matter) may rea- 
sonably be described in the non-quantised Thomson limit0 After decoupling, there is very little interaction 



The 1 + 3 covariant treatment described in tliis section also neglects polarisation effects. 
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between matter and the CBR so we can use Liouville's equation: 

C[f]=0 => L{f) = 0. (234) 
The energy-momentum tensor of the photons is 

r;^= [ p-p'fn, (235) 

where tt is the volume element in the momentum tangent space T^^. This satisfies the conservation equations 
( [35] l and ( [36l ). and is part of the total energy-momentum tensor T"^^ determining the space-time curvature by 
Einstein's field equations 

Exercise: The same theory applies to particles with non-zero rest-mass, (i) What is the form of Boltz- 
mann's collision term for binary particle collisions? (ii) Show that energy-momentum conservation in in- 
dividual collisions will lead to conservation of the particle energy-momentum tensor T"''. (iii) Using the 
appropriate integral definition of entropy, determine an H-theorem for this form of collision, (iv) Under what 
conditions can equilibrium exist for such a gas of particles, and what is the equilibrium form of the particle 
distribution function? (See ISj.) 

8.2 Angular harmonic decomposition 

In the 1 + 3 covariant approach of 1 2321 12331 Fl the photon 4-momentump° (where paP° = 0)is split relative 
to an observer moving with 4-velocity as 

= S {u" + e") , eae"" = 1 , Cau" = , (236) 

where E — {—paU°') is the photon energy and e° = E^^ h°-bP^ is the photon's spatial propagation 
direction, as measured by a matter-comoving (fundamental) observer. Then the photon distribution function 
is decomposed into 1 + 3 covariant Iiarmonics via the expansion 1232112371 

f{x,p) = f{x, E,e)^J2 ^) e^'=F + F^e'' + Fa^e'^e'' + • • • = ^ F^.e^-^^^ , (237) 

where e^^ = e^^e"^ • • • e°*, and e*-"^*^ is the symmetric trace-free part of e^'. The 1 + 3 covariant dis- 
tribution function anisotropy multipoles Fai are irreducible since they are Projected, Symmetric, and 
Trace-Free ('PSTF'), i.e., 

Fa-.-b = F(^a--b) ^ Fa-.-b ~ -F'(a---6) : Fa-.-bU^ = = Fa-.-bch!"^ Fa, — Fi^At) ■ 

They encode the anisotropy structure of the photon distribution function in the same way as the usual splier- 
ical Iiarmonic expansion 

but with two major advantages: (a) the Fa, are 1 + 3 covariant, and thus independent of any choice of 
coordinates in momentum space, unlike the /™; (b) Fa, is a rank-^ tensor field on space-time for eachfixed 
E, and directly determines the £-multipole of radiation anisotropy after integration over E. The multipoles 
can be recovered from the photon distribution function via 

r 2^f£')^ 
Fa,^ AJ^ J f{x,e)e^A,)dn, with A, = 4^ (2/+ I)! ' ^^^^^ 

'^Based on the Ph.D. Thesis of R Treciokas, Cambridge University, 1972, combined with the covariant formalism of F A E Pirani 
1234 1 (see also K S Thome 12351) by G F R Ellis when in Hamburg (1st Institute of Theoretical Physics) in 1972; see also 1691 1701 . A 
similar formahsm has been developed by M L Wilson t236j . 
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where dJl = d e is a soUd angle in momentum space. A further useful identity is 02321 

47r f ^odd, 
^' dn = I (239) 

^ [ ^(Qia2 ;ja3a4 . . .^jQ^-ia^) £ gVeU . 

The first three multipoles determine the radiation energy-momentum tensor, which is, from ( I235l l and 
(O, 

T:\x) = Jp-p'fix,p) d'p = M„ + i /i"") + 2 q^^ u") + , (240) 

where d^p = i? dE dil is the covariant momentum space volume element on the future null cone at the event 
X. It follows from ( 1237b and (1240b that the dynamical quantities of the radiation (in the w°-frame) are: 

= 47r / E^FdE , ^ — E^ dE , T:f = — E^ F"^ dE . (241) 

We extend these dynamical quantities to all multipole orders by defining the 1 + 3 covariant brightness 
anisotropy multipoleS 12371 



Ua, = / E"" Fa, dE = U^a,) , (242) 

so that n = Mfl/47r, n° = 3 q^J-ii: and = 15 tt^VStt. 
Writing Boltzmann's equation ( 1232b in the form 

|^pV„(/)-r%,pV^^-C[/], (243) 

the collision term is also decomposed into 1 + 3 covariant harmonics: 

C[f]^Y. ^-4. E) e^'-=b + bac^ + bate'^e" + ■■■ , (244) 

where the 1 + 3 covariant scattering multipoles = ^(A^) encode irreducible properties of the particle 
interactions. Then Boltzmann's equation is equivalent to an infinite hierarchy q/ 1 + 3 covariant multipole 
equations 

LA,{x,E) = bAAFA,Ax,E)] , 

where La, = ^(a^) ^re the anisotropy multipoles of df /dv, and will be given in the next subsection. These 
multipole equations are tensor field equations on space-time for each value of the photon energy E (but note 
that energy changes along each photon path). Given the solutions FAe{x, E) of the equations, the relation 
(1237b then determines the full photon distribution function f{x,E, e) as a scalar field over phase space. 



8.3 Non-linear 1 + 3 covariant multipole equations 

The full Boltzmann equation in photon phase space contains more information than necessary to analyze 
radiation anisotropies in an inhomogeneous Universe model. For that purpose, when the radiation is close to 
black body, we do not require the full spectral behaviour of the photon distribution function multipoles, but 
only the energy-integrated multipoles. The monopole leads to the average (all-sky) temperature, while the 
higher-order multipoles determine the temperature anisotropies. The GIC definition of the average temper- 
ature T is given according to the Stefan-Boltzmann law by 

fi„ {x)^iTT I E^ F{x,E)dE ^rT^{x) , (245) 



^^'^Because photons are massless, we do not need the complexity of the moment definitions used in 12321 . From now on, all energy 
integrals will be understood to be over the range < E < cx). 
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where r is the radiation constant. If / is close to a Planckian distribution, then T is the thermal black body 
average temperature. But note that no notion of background temperature is involved in this definition. There 
is an all-sky average implied in ( 12451 ). Fluctuations across the sky are measured by energy-integrating the 
higher-order multipoles (a precise definition is given below), i.e., the fluctuations are determined by the 11^^ 
(£ > 1) defined in (l242l) . 



The form of C[f] in ( I244l l shows that 1 + 3 covariant equations for the temperature fluctuations arise from 
decomposing the energy-integrated Boltzmann equation 



dv 



E^ C[f] dE 



(246) 



into 1 + 3 covariant multipoles. We begin with the right-hand side, which requires the 1 + 3 covariant form 
of the Thomson scattering term. Defining the 1 + 3 covariant energy-integrated scattering multipoles 



E^ bA, dE = K 



{At) 



we find that EH 



K 



- n,a^ [ n" - 4 nv^ - I U-\b ] + 0[3] 



10 2 B 7 B B 

jjahc _ 3 n^''^^) - 4 U'''""^V^d 1 + C[3] , 



and, for i > 3, 



n 



1 



2^ + 3 



0[3], 



(247) 

(248) 
(249) 

(250) 
(251) 

(252) 



where the expansion is in terms of the peculiar velocity Vg of the baryons relative to the reference frame 
Parameters are the free electron number density and the Thomson scattering cross section ct^, the 

latter being proportional to the square of the classical electron radius 12301 12311 . The first three multipoles 
are affected by Thomson scattering differently than the higher-order multipoles. 

Equations ( I248I )-( I252I ). derived in [214], are a non-linear generalisation of the results given by Challinor 
and Lasenby B238L They show the coupling of baryonic bulk velocity to the radiation multipoles, arising 
from 1 + 3 covariant non-linear effects in Thomson scattering. If we linearise fully, i.e., neglect all terms 
containing Vg except the term in the dipole K'^\ which is first-order, then our equations reduce to those 
in 0238L The generalised non-linear equations apply to the analysis of 1 + 3 covariant second-order effects 
on an FLRW background, to first-order effects on a spatially homogeneous but anisotropic background, and 
more generally, to any situation where the baryonic frame is in non-relativistic motion relative to the funda- 
mental u" -frame. 



Next we require the anisotropy multipoles Lai of df /dv. These can be read off for photons directly 
from the general expressions in 12321 . which are exact, 1 + 3 covariant, and also include the case of massive 
particles. For clarity and completeness, we outline an alternative, 1 + 3 covariant derivation (the derivation 
in 1232 1 uses tetrads). We require the rate of change of the photon energy along null geodesies, given by 
16911701 

^ = ~ [ i e + (u.e'') + {aabe'^e") ] E' , (253) 
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which follows directly from p^'Vbp"' — with (I236I I and ( fTTT l. Then 



where a prime denotes d/dE. The first term is readily put into irreducible PSTF form using identities in 
02321 . p 470. In the second term, when the round brackets are expanded, only those terms with at most one 
yflr survive, and 

U^'Fa... - - U^'Fa... , u''V''Fb... = - ( i e /l"^^ + (t"^ + 7?"''= LO, ) Ffo... . 

Thus, the 1 + 3 covariant anisotropy multipoles La^ of df /dv are 

- [E'^' FA.b]' u" - E' [E'-' F.A.J u^, 



f,P^ P c (^+l)(^ + 2) ^-(£+2) [77^+3 F ]' 6 

2^ 



(2^ + 3) 
E-^ hA, ■ 



E^'^F^A,., ^a,)' - E'-^ [E^-' F^A,.,]' <Ja,.,a,) (254) 



This regains the result of |I232| ( equation (4. 12) ) in the massless case. The form given here benefits from 
the streamlined version of the 1 + 3 covariant formalism. We re-iterate that this result is exact and holds for 
any photon or (massless) neutrino distribution in any space-time. 

We now multiply (1254b by E'^ and integrate over all photon energies, using integration by parts and the 
fact that E" Fa--- as E ^ oo for any positive n. We obtain the evolution equations that determine the 
brightness anisotropics multipoles 11^^ : 

n(A,> + 1 e n^, + V(,,n^,_,) + v^n^.b (255) 

(f -d- i)(p — 2) 

~ (21 + 3) + + 3) ^(a. n^.-.> + e^'Vbcia. ^A,-.)' 

(£-l){£+l)ii + 2) , U , 

" (2£ + 3)(2£+5) " n^.bc + ^^j:^^ (a.n^,_,).-(^ + 2)a(,,,,_,n^,_,) = Ka,. 

Once again, this is an exact result, and it holds also for any collision term, i.e., any Kai- For decoupled 
neutrinos, we have K^'^ = in this equation. For photons undergoing Thomson scattering, the right-hand 
side of ( |255] | is given by (I252l i. which is exact in the kinematical and dynamical quantities, but first-order in 
the relative baryonic velocity. The equations (I252l l and (1255b thus constitute a non-linear generalisation of 
the FLRW-linearised case given by Challinor and Lasenby [238] . 

The monopole and dipole of equation (1255b give the evolution equations for the energy and momentum 
densities: 

ri + |en + iVan'' + |?ian'' + ^(Tafcn'^'' = (256) 
ri<''> + 1 e n" + v"n + § VbH'^'' 

+ f U6n"'' + 4nu'' + ?7°'"=a;fcnc + cr"''nfc = K\ (257) 



8 CBR ANISOTROPIES 



68 



(these are thus the equations giving the divergence of T"'' in (I235l l). For photons, the equations are: 

^i, + |e^I,+\7aq^ + 2uaq^ + aabnf = n^a^ I ^^nV^ - 1>Ba) + 0[3] , (258) 

+ a\q'^+ri-''-LJtq„c + UbK'' = ^e'^t H - € + ^nb) + 0[3] , (259) 

(the present versions of the fluid energy and momentum conservation equations ( |35] | and (|36]|). 

The non-Hnear dynamical equations are completed by the energy-integrated Boltzmann multipole equa- 
tions. For photons, the quadrupole evolution equation is 

^U^b) ^iQ^ab^^ ^ab ^ I y (a^b) ^ 8^ ^ ^^abc 

= - n,a^ ( A _ 1 yb) _ 8^ n-'.c _ 2 „(a (^[3] ^ (260) 

(a fluid description gives no analogue of this and the following equations). The higher-order multipoles 
(£ > 3) evolve according to 

- l)(^ + 1)(£ + 2) ^ 5e - + 2) a^'^^'^-i n-4-^> 



(2£ + 3)(2^ + 5) (2^ + 3) 



2^ + 3 



jApa , 



0[3] . (261) 



For ^ = 3, the J\^^'^-^v'^J^ term on the right-hand side of equation ( 12611 ) must be multiplied by |. For 
neutrinos, the equations are the same except without the Thomson scattering terms. Note that these equations 
link angular multipoles of order £ — 2, £ — 1, ^, ^ + 1, ^ + 2, i.e., they link^ve successive harmonic terms. 
This is the source of the harmonic mixing that occurs as the radiation propagates. 

These equations for the radiation (and neutrino) multipoles generahse the equations given by Challinor 
and Lasenby |238l, to which they reduce when we remove all terms Oiev^) and ©(e^). In this case, on 
introducing a FLRW-Iinearisation, there is major simplification of the equations: 

/i„ + |e/i„+V„g„" « 0, (262) 
9i'^^ + |09." + |A^>" + iVV«+V<'' « n,aA\^n^l~ql). (263) 



and, for ^ > 3, 



ri^-^^^ + f e n^^ + v^^^n-^'-i^ + \ ^^^^^^ « - n^a^ W 



(265) 



Note that these equations now link only angular multipoles of order £ ~ 1, £, £ + 1, i.e., they link three 
successive terms. This is a major qualitative difference from the full non-linear equations. 
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8.4 Temperature anisotropy multipoles 

Finally, we return to the definition of temperature anisotropies. As noted above, these are determined by the 
■ We define the temperature fluctuation t{x, e) via the directional temperature which is determined by 
the directional bolometric brightness: 



T{x,e)^T{x)[l + T(x,e) 



^ J E^f{x,E,e)dE 



1/4 



(266) 



This is a GIC definition which is also exact. We can rewrite it explicitly in terms of the 11^^ : 

1/4 



r(x, e) 



^« e>i 



- 1 = Tae^ + r^fcc^e" + • • • = ^ r^, (x) e^' 



(267) 



e>i 



In principle, we can extract the 1 + 3 covariant ureducible PSTF temperature anisotropy multipoles ta^ 

T(^Ae) by using the inversion (I238l l: 



TAeix) ^ I T{x,e)e(Ai)d^l 



(268) 



In the almost-FLRW case, when r is ©[1], we regain from (1267b the linearised definition given in II239I : 

Ha, , (269) 



where £ > 1. In particular, the dipole and quadrupole are 

in" 



and t"'' w 2_ 

4m« '^fJ-n 



(270) 



We can normalise the dynamical brightness anisotropy multipoles 11^^ of the radiation to define the 
dimensionless 1 + 3 covariant brightness temperature anisotropy multipoles (i > 1) 



(271) 



Thus, the — T(^Ai) equal to the temperature anisotropy multipoles plus non-hnear corrections. In 
terms of these quantities, the hierarchy of radiation multipoles becomes: 



T 
T 



_4 l + ie 



T 



5 -' 



I 2 _ _ i-j-ab a rj-b 2 • n~o,b ^abc ^ 



T 



(272) 



0[3] , (273) 



T 



bed 



(274) 
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and, for £ > 3: 



1 



2^ + 3 



(2^ + 3) ' 



(^-l)(£+l)(^ + 2) 
(2<' + 3) (2£+5) 



(2^ + 3) 



T 



„ rfAlhc 



■0[3]. 



(2^ + 3) 



(275) 



For ^ = 3, the Thomson scattering term 'j'i^e-iy'^i) must be muhiphed by a factor |. 

The 1 + 3 covariant non-Hnear muhipole equations given in this form show more clearly the evolution 
of temperature anisotropies (including the monopole, i.e., the average temperature T), in general linking five 
successive harmonics. Although the only determine the actual temperature fluctuations t^, to linear 
order, they are a useful dimensionless measure of anisotropy. Furthermore, equations ( I272I )-( I275I ) apply as 
the evolution equations for the temperature anisotropy multipoles when the radiation anisotropy is small (i.e., 
~ TAc), but the space-time inhomogeneity and anisotropy are not restricted. This includes the particular 
case of small CBR anisotropies in general Bianchi models, or in perturbed Bianchi models. 

FLRW-linearisation, i.e., the case when only first-order effects relative to the FLRW Umit are consid- 
ered, reduces the above equations to the linearised form, generically linking three successive harmonics: 



T 
T 



T 



b 



(276) 

(277) 
(278) 



and, for £ > 3: 



■AM) 



(2^+3) 



A,b 



(279) 



These are GIC multipole equations leading to the Fourier mode formulation of the energy-integrated Boltz- 
mann equations used in the standard literature (see, e.g., Il240j and the references therein) when they are 
decomposed into spatial harmonics and associated wavelengths]^ This then allows examination of diffusion 
effects, which are wavelength-dependent. 

These linearised equations, together with the linearised equations governing the kinematical and grav- 
itational quantities, may be 1 + 3 covariantly split into scalar, vector and tensor modes, as described in 
02211 12381 . The modes can then be expanded in 1 + 3 covariant eigentensors of the matter-comoving Lapla- 
cian, and the Fourier coefficients obey ordinary differential equations. 



Exercise: Determine the resulting hierarchy of mode equations. Show from these equations that there is 
a wavelength As such that for shorter wavelengths the perturbations are heavily damped (the physical reason 
is photon diffusion). 

Numerical integrations of these equations are performed for scalar modes by Challinor and Lasenby 112381 . 
with further analytic results given in |241 1. These are the Sachs- Wolfe family of integrations, of fundamental 



'They contain the free-streaming subcase for = 0. 
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importance in determining CBR anisotropies, as discussed in many places; see, e.g., Hu and Sugiyama 12401 . 
Hu and White II242II (the GIC version will be given in the series of papers by Gebbie, Maartens, Dunsby and 
Ellis). However, they assume an almost-FLRW geometry. We turn now to justifying that assumption. 

8.5 Almost-EGS-Theorem and its applications 

One of our most important understandings of the nature of the Universe is that at recent times it is well- 
represented by the standard spatially homogeneous and isotropic FLRW models. The basic reason for this 
belief is the observed high degree of isotropy of the CBR, together with a fundamental result of Ehlers, Geren 
and Sachs |69| (hereafter 'EGS'), taken nowadays (see, e.g., |3|) to establish that the Universe is almost- 
FLRW at recent times (i.e., since decoupling of matter and radiation). 

The EGS programme can be summarized as follows: Using (a) the measured high isotropy of the CBR 
at our space-time position, and (b) the Copernican assumption that we are not at a privileged position in the 
Universe, the aim is to deduce that the Universe is accurately FLRW. EGS gave an exact theorem of this kind: 
if a family of freely-falling observers measure self-gravitating background radiation to be everywhere exactly 
isotropic in the case of non-interacting matter and radiation, then the Universe is exactly FLRW. This is taken 
to establish the desired conclusion, in view of the measured near-isotropy of the radiation at our space-time 
location. However, of course the CBR is not exactly isotropic. Generally, we want to show stability of 
arguments we use [243 J ; in this case, we wish to do so by showing the EGS result remains nearly true if 
the radiation is nearly isotropic, thus providing the foundation on which further analyses, such as that in the 
famous Sachs-Wolfe paper II223I . are based. More precisely, we aim to prove the following theorem II244I . 

Almost-EGS-Theorem: If the Einstein-Liouville equations are satisfied in an expanding 
Universe model, where there is present pressure-free matter with 4-velocity vector field 
(Mqu" = — 1) such that (freely-propagating) background radiation is everywhere almost-isotropic 
relative to u"" in some domain U, then the space-time is almost-FLRW in U . 

This description is intended to represent the situation in the Universe since decoupling to the present 
day. The pressure-free matter represents the galaxies on which fundamental observers live, who measure the 
radiation to be almost isotropic. 

8.5.1 Assumptions 

In detail, we consider matter and radiation in a space-time region U . In our application, we consider U to 
be the region within and near our past light cone from decoupling to the present day (this is the observable 
space-time region where we would like to prove the Universe is almost-FLRW; before decoupling a different 
analysis is needed, for collisions dominate there, and we do not have sufficient data to comment on the situa- 
tion far from our past light cone). Our assumptions are that, in the region considered, 

(1) Einstein's field equations are satisfied, with the total energy-momentum tensor Tab composed of non- 
interacting matter and radiation components: 

yab ^ j-ab _^ j,ab ^ ^ ^j,ab ^ q ^ ^ ^j,ab ^ q (280) 

The independent conservation equations express the decoupling of matter from radiation; the matter energy- 
momentum tensor is T^^ — pu°' u^, (fij^j — p is the matter energy density) and, at each point, the radiation 
energy-momentum tensor is T^^ = J f 

where tt is the momentum space volume element. The only 
non-zero energy-momentum tensor contribution from the matter, relative to the 4-velocity u'^, is the energy 
density; so the total energy density is 

p = fi^+p = O[0], (281) 
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while all other energy-momentum tensor components are simply equal to the radiation contributions. Without 
confusion we will omit a subscript 'i?' on these terms. 

(2) The matter 4-velocity field is geodesic and expanding: 

e = O[0]> 0, (282) 

(the first requirement in fact follows from momentum conservation for the pressure-free matter; see (|38]|). 
The assumption of an expanding Universe model is essential to what follows, for otherwise there are counter- 
examples to the result f2451 . 

(3) The radiation obeys Liouville's equation ( 1234b : L{f) = 0, where L is the Liouville operator This 
means that there is no entropy production, so that q°- and iTab are not dissipative quantities, but measure the 
extent to which / deviates from isotropy. 

(4) Relative to — i.e., for all matter-comoving observers — the photon distribution function / is 
almost isotropic everywhere in the region U. Formally, in this region, F — F{x, E) and its time derivatives 
are zeroth-order, while Fa,, = Fa,, {x, E) plus their time and space derivatives are at most first-order for 
f > 0. In brief: 

F,F^O[Q], Fa,. F^a,), VaFA, = , (283) 
(and we assume all higher derivatives of Fa, are also 0[1] for I > 0). 

It immediately follows, under reasonable assumptions about the phase space integrals of the harmonic 
components, that the scalar moments are zeroth-order but the tensor moments and their derivatives are first- 
order We assume the conditions required to ensure this are true, i.e., we additionally suppose 

fi^,p^=O[0], g(a), Vag6 = 0[l], TT.fe, 7T(,b) , VaTT,,, = O [1] (284) 

are satisfied as a consequence of (I283l l, and that the same holds for the higher time and space derivatives of 
TTab and for all higher-order moments (specifically, those defined in (I292l i below). The^rif aim now is to 
show that 

(a) the kinematical quantities and the Weyl curvature are almost-FLRW, i.e., aab and uj°' are first-order, 
which then implies that Eab and Hab are also first-order. The second aim is to show that then, 

(b) there are coordinates such that the metric tensor takes a perturbed RW form. 

8.5.2 Proving almost-FLRW kinematics 

Through an appropriate integration over momentum space, the zeroth harmonic of Liouville's equation (I234l l 
gives the energy conservation equation for the radiation ( cf. (125 8l l ), 

/i„ + I e + V„<z" + aabTT'"'' - , (285) 

while the first harmonic gives the momentum conservation equation ( cf. ( 12591 ) ), 

q{a) + 4 ^ 1 ya^^ ^ y^^ab _^ ^ ^b ^ ,^abc = Q , (286) 

which implies by ( I284l i. 

V,M„ = 0[l] . (287) 

Taking spatial derivatives of (I286l l. the same result follows for the higher spatial derivatives of ixb.\ in partic- 
ular, its second derivatives are at most first-order. 
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Now the definition of the V-derivative leads to the identity (|27] |. giving 

(VaVb - Vf,Va)AiH = 2TJabcUj'fl^ . (288) 

In an expanding Universe model, by ( I284l i. the energy conservation equation ( I285l l for the radiation shows 

Un + lQf^n=0[l]. (289) 



Thus, because the model is expanding, fiR is zeroth-order However, the left-hand side of ( 12881 ) is first-order; 
consequently, by the higher-derivative version of ( 12871 ). 

w° = 0[l] . (290) 

The second harmonic of Liouville's equation ( 1234b leads, after an appropriate integration over momentum 
space, to an evolution equation for the anisotropic stress tensor iTab (which is the present version of ( 12601 )): 



r^'^''^ + I e vr'^'' + ^ cr'^'' + f V^'^g''^ + + 2 tj" 77^^^" vr^^'' + f ^c^" tt^^" ^crf H'^'"'' = , (291) 

315 



where (cf. ( |242] |) 



35 

Consequently, by ( |290] i. (l283T i and ( |284] |. 



J^fc^—/ V ^F"'"' dE , iiabcd ^ E^F'''""^dE. (292) 



^Tafc - , (293) 



and (on taking derivatives of the above equation) the same is true for its time and space derivatives. Then, to 
first order, the evolution equations become 

fl, + |e^l, + Vaq'' ~ 0, (294) 

gM + le^a^ 1 V^/i^+Vfevr"'' « 0, (295) 
-k^"''^ + f e tt'''' + ^ ^„ a'''' + I V^'^q''^ + J"'' « , (296) 

showing the equations can only close at first order ;/ we can argue that J°'' = 0[2]. Any such approximation 
must be done with utmost care because of the radiation multipole truncation theorem (see below). 

It now follows from the shear propagation equation (ISTT i and the iJQ(,-equation ( l34b that all the Weyl 
curvature components are also at most first-order: 

Eab - , Hab =0[1]. (297) 
Consequently, the (div i?) -equation ( |45T l shows that 

Xa = Vaf^ = 0[l] ^ VaP = 0[l], (298) 

by ( 12811 ). Finally, the (OQ!)-equation ([32]), or spatial derivatives of ( |289t , show that 

Za = VaO = 0[1] . (299) 

This establishes that the kinematical quantities for the matter flow and the Weyl curvature are almost-FLRW: 
all the quantities that vanish in the FLRW case are at most first-order here. Thus, the only zeroth-order 1 + 3 
covariantly defined quantities are those that are non-vanishing in FLRW models. 
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8.5.3 Proving almost-FLRW dynamics 

It follows that the zeroth-order equations governing the dynamics are just those of a FLRW model. Because 
the kinematical quantities and the Weyl curvature are precisely those we expect in a perturbed FLRW model, 
we can linearise the 1 + 3 covariant equations about the FLRW values in the usual way, in which the back- 
ground model will — as we have just seen — obey the usual FLRW equations. Hence, the usual linearised 
1 + 3 covariant FLRW perturbation analyses can be applied II 1971 12211 . leading to the usual results for growth 
of inhomogeneities in almost-FLRW Universe models. 

8.5.4 Finding an almost-RW metric 

Given that the kinematical quantities and the Weyl curvature take an almost-FLRW form, the key issue in 
proving existence of an almost-RW metric is choice of a time function to use as a cosmic time (in the realistic, 
inhomogeneous Universe model). The problem is that although we have shown the vorticity will be small, it 
will in general not be zero; hence, there will be no time surfaces orthogonal to the matter flow lines fSHS). 
The problem can equivalently be viewed as the need to find a vorticity-free (i.e., hypersurface-orthogonal) 
congruence of curves to use as a kinematical reference frame, which, if possible, one would like to also be 
geodesic. As we cannot assume the matter 4-velocity fulfills this condition, we need to introduce another 
congruence of curves, say u", that is hypersurface-orthogonal and that does not differ too greatly from 
so the matter is moving slowly (non-relativistically) relative to that frame. Then we have to change to that 
frame, using a 'hat' to denote its kinematical quantities. (See II244I .') 

8.5.5 Result 

This gives the result we want; we have shown that freely propagating almost-isotropic background radiation 
everywhere in a region U implies the Universe model is almost-FLRW in that region. Thus, we have proved 
the stability of the Ehlers, Geren and Sachs 1 69 1 result. This result is the foundation for the important analysis 
of Sachs and Wolfe II223I and all related analyses, determining the effect of inhomogeneities on the CBR by 
integration from last scattering till today in an almost-FLRW Universe model, for these papers start off with 
the assumption that the Universe is almost-FLRW since decoupling, and build on that basis. 

It should also be noted that the assumption that radiation is isotropic about distant observers can be par- 
tially checked by testing how close the CBR spectrum is to black body in those directions where we detect the 
Sunyaev-ZeI'dovich effect HI 041 (see also Goodman [246]). Since that effect mixes and scatters incoming 
radiation, substantial radiation anisotropics relative to those clusters of galaxies inducing the effect would 
result in a significant distortion of the outgoing spectrum. 

Exercise: How good are the limits we might obtain on the CBR anisotropy at distance points on our past 
light cone by this method? Are there any other ways one can obtain such limits? 

8.6 Other CBR calculations 

One can extend the above analysis to obtain model-independent limits on the CBR anisotropy (Maartens et 
al II239II ). However, there are also a whole series of model-dependent analyses available. 

8.6.1 Sachs- Wolfe and related effects 

We will not pursue here the issue of Sachs- Wolfe effect II223I type calculations, integrating up the redshift 
from the surface of last scattering to the observer and so determining the CBR anisotropy, and, hence, cos- 
mological parameters |247|, because they are so extensively covered in the literature. The GIC version is 
covered by Challinor and Lasenby l,238j and developed in depth in the series of papers by Gebbie, Dunsby, 
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Maartens and Ellis II214II2371I241I . However, a few comments are in order. 

One can approach this topic 1 + 3 covariantly by a direct generalisation of the photon type of calculation 
(see 112251 1226II ). or by use of the relativistic kinetic theory formalism developed above, as in the papers just 
cited. Two things need to be very carefully considered. These are, 

(i) The issue of putting correct limits on the Sachs-Wolfe integral. This integral needs to be taken to the 
real surface of last scattering, not the background surface of last scattering (which is fictitious). If this is not 
done, the results may be gauge-dependent. One should place the limits on this integral properly, which can 
be done quite easily by calculating when the optical depth due to Thomson scattering is unity II2241 12481 . 

(ii) In order to obtain solutions, one has somehow to cut off the harmonic series, for otherwise there is an 
infinite regress whereby higher-order harmonics determine the evolution of lower-order ones, as is clear from 
the multipole equations above. However, truncation of harmonics in a kinetic theory description needs to be 
approached with great caution, because of the following theorem for collision-free radiation II233II : 

Radiation Multipole Truncation Theorem: In the exact Einstein-Liouville theory, trun- 
cation of the angular harmonics seen by a family of geodesic observers at any order whatever 
leads to the vanishing of the shear of the family of observers; hence, this can only occur in highly 
restricted spaces. 

The proof is given in 1,233 1. This is an exact result of the full theory. However, it remains true in the linear 
theory if the linearisation is done carefully, even though the term responsible for the conclusion is second- 
order and so is omitted in most linearised calculations. 

The point is that in that equation, at the relevant order, there are only second-order terms in the equation, 
so one cannot drop the terms responsible for this conclusion (one can only drop them in an equation where 
there are also linear terms, so the second-order terms are negligible compared with the linear terms; that is 
not the case here). Thus, a proper justification for the acceptability of some effective truncation procedure 
relies on showing how a specific approximation procedure gives an acceptable approximation to the results 
of the exact theory, despite this disjuncture. 

Exercise: Give such a justification. 
8.6.2 Other models 

As mentioned previously in these lectures, there is an extensive literature on the CBR anisotropy: 

• in Bianchi (exact spatially homogeneous) models; see, e.g., II 1901 - 11 19211 . 

• in Swiss-Cheese (exact inhomogeneous) models; see, e.g., P16H and II142I , 

• in small universes (FLRW models with compact spatial topology that closes up on a small scale, so 
that we have seen round the universe already since decoupling ll89l ); see, e.g., 124911 and ll92l . 

These provide important parametrized sets of models that one can use to test and exploit the restrictions the 
CBR observations place on alternative Universe models that are not necessarily close to the standard models 
at early or late times, but are still potentially compatible with observations, and deserve full exploration. 

In each case, whether carrying out Sachs-Wolfe type or more specific model calculations, we get tighter 
limits than in the almost-EGS case, but they are also more model-dependent. In each case they can be carried 
out using the 1 + 3 covariantly defined harmonic variables and can be used to put limits on the 1 + 3 covari- 
antly defined quantities that are the theme of this article. 
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The overall conclusion is that — given the Copernican assumption underlying use of the almost-FLRW 
models — the high degree of observed CBR isotropy confirms use of these models for the observable Universe 
and puts limits on the amplitude of anisotropic and inhomogeneous modes in this region (i.e., within the 
horizon and since decoupling). Other models remain viable as representations of the Universe at earlier and 
later times and outside the particle horizon. 

9 Conclusion and open issues 

The sections above have presented the 1 + 3 covariant variables and equations, together with the tetrad equa- 
tions, a series of interesting exact cosmological models, and a systematic procedure for obtaining approxi- 
mate almost-FLRW models and examining observations in these models. It has been shown how a number 
of anisotropic and inhomogeneous cosmological models are useful in studying observational limits on the 
geometry of the real Universe, and how interesting dynamical and observational issues arise in considering 
these models; indeed issues such as whether inflation in fact succeeds in making the Universe isotropic or not 
cannot be tackled without examining such models II250II . 

9.1 Conclusion 

The standard model of cosmology is vindicated in the following sense (cf. 1941 14l1 ): 

• The Universe is expanding and evolving, as evidenced by the {magnitude, redshift)-relation for galaxies 
and other sources, together with number count observations and a broad compatibility of age estimates; 

• It started from a hot big bang early stage which evolved to the presently observed state, this early era 
being evidenced by the CBR spectrum and concordance with nucleosynthesis observations; 

• Given ^ Copernican assumption, the high degree of isotropy of the CBR and other observations sup- 
port an almost-FLRW (nearly homogeneous and isotropic) model within the observable region (inside 
our past light cone) since decoupling and probably back to nucleosynthesis times; 

• There are globally inhomogeneous spherically symmetric models that are also compatible with the 
observations if we do not introduce the Copernican assumption II144L and there are inhomogeneous 
and anisotropic modes that suggest the Universe is not almost-FLRW at very early and very late times 
ll5nil96l . and on very large scales (outside the particle horizon) |66|. 

The issue of cosmological parameters may be resolved in the next decade due to the flood of new data com- 
ing in — from deep space number counts and redshift measurements, observations of supernovae in distant 
galaxies, measurement of strong and weak gravitational lensing, and CBR anisotropy measurements (see 
Coles and Ellis 1.41 J for a discussion). 

The further extensions of this model proposed from the particle physics side are at present mainly com- 
patible but not yet as compelling from an observational viewpoint, namely: 

• An early inflationary era helps resolve some philosophical puzzles about the structure of the Universe 
related specifically to its homogeneity and isotropy Il64l - ll66l . but the link to particle physics will not 
be compelling until a specific inflaton candidate is identified; 

• Structure formation initiated by quantum fluctuations at very early stages in an inflationary model are a 
very attractive idea, with the proposal of a CDM-dominated late Universe and associated predictions for 
the CBR anisotropy giving a strong link to observations that may be confirmed in the coming decade; 
however, theoretical details of this scheme, and particularly the associated issues of biasing and the 
normalisation of matter to radiation, are still to be resolved; 

'^''Necessarily philosophically based 11501 . 
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• The density of matter in the Universe is probably below the critical value predicted by the majority of 
inflationary modelsJ4rj; compatibility with inflation may be preserved by either moving to low density 
inflationary modelo or by confirmation of a cosmological constant that is dynamically dominant at 
the present time@ the latter proposal needs testing relative to number counts, lensing observations, 
and structure formation scenarios; 

• The vibrant variety of pre-inflationary proposals involve a huge extrapolation of presently known 
physics to way beyond the testable domain and, given their variety, do not as yet give a compelling 
unique view of that era. 

Hence, it is suggested BTI that these proposals should not be regarded as part of the standard model but rather 
as interesting avenues under investigation. 

9.2 Open issues 

Considered from a broader viewpoint, substantial issues remain unresolved. Among them are, 

(1) The Newtonian theory of cosmology is not yet adequately resolved. Newtonian theory is only a good 
theory of gravitation when it is a good approximation to General Relativity; obtaining this limit in non-linear 
cosmological situations raises a series of questions and issues that still need clarifying ll24l . particularly re- 
lating to boundary conditions in realistic Newtonian cosmological models. 

(2) We have some understanding of how the evolution of families of inhomogeneous models relates to 
that of families of higher symmetry models. It has been indicated that a skeleton of higher symmetry mod- 
els seems to guide the evolution of lower symmetry models in the state space (the space of cosmological 
space-times) |51 1. This relation needs further elucidation. Also, anisotropic and inhomogeneous inflationary 
models are relatively little explored and problems remain II 1841 12531 ; 

(3) We need to find a suitable measure of probabiUty in the full space of cosmological space-times, and 
in its involutive subspaces. The requirement is a natural measure that is plausible. Some progress has been 
made in the FLRW subcase 12541 - 112571 . but even here it is not definitive. Closely related to this is the issue 
of the stability of the results we derive from cosmological modelling II243I ; 

(4) We need to be able to relate descriptions of the same space-time on different scales of description. 
This leads to the issue of averaging and the resulting effective (polarization) contributions to the energy- 
momentum tensor, arising because averaging does not commute with calculating the field equations for a 
given metric I.88J , and we do not have a good procedure for fitting a FLRW model to a lumpy realistic model 
of the Universe II258I . It includes the issue (discussed briefly above) of how the abnost-every where empty 
real Universe can have dynamical and observational relations that average out to high precision to the FLRW 
relations on a large scale. 

(5) Related to this is the question of definition of entropy for gravitating systems in general^ and cosmo- 
logical models in particular This may be expected to imply a coarse-graining in general, and so is strongly 
related to the averaging question. It is an important issue in terms of its relation to the spontaneous formation 
of structure in the early Universe, and also relates to the still unresolved arrow of time problem [259], which 
in turn relates to a series of further issues concerned with the effect on local physics of boundary conditions 
at the beginning of the Universe (see, e.g., Il28l ). 



■^'For an early proposal, see 12511 and (73). 

'^^But way below that predicted by present field theory 12521 . 

'^^In the case of black holes, there is a highly developed theory; but there is no definition for a general gravitational field; see, e.g., 
12591 . 
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(6) One can approach relating a model cosmology to astronomical observations by a strictly observational 
approach ||98l[T0 ri. as opposed to the more usual model-based approach as envisaged in the main sections 
above and indeed in most texts on cosmology. Intermediate between them is a best-fitting approach II258L 
Use of such a fitting approach is probably the best way to tackle modelUng the real Universe II260I . but it is 
not yet well-developed. 

(7) Finally, underlying all these issues is the series of problems arising because of the uniqueness of the 
Universe, which is what particularly gives cosmology its unique character and underlies the special problems 
in cosmological modelling and application of probability theory to cosmology |261 1. Proposals to deal with 
this by considering an ensemble of Universe models realized in one or other of a number of possible ways are 
in fact untestable and, hence, of a metaphysical rather than physical nature; but this needs further exploration. 

There is interesting work still to be done in all these areas. It will be important in tackling these issues 
to, as far as possible, use gauge-invariant and 1 + 3 covariant methods, because coordinate methods can be 
misleading. 
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